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Executive  Summary 

Motivation  The  efficiency  of  modem  Hall  thrusters  is  impacted  by  the  ability  to  maintain  a  high 
ratio  of  the  ion  to  electron  current  density,  Ji/Je,  at  the  exit  of  the  acceleration  channel.  In  principal, 
there  is  no  limit  to  this  ratio,  since  the  electrons  that  are  necessary  to  maintain  the  current  balance  at 
the  anode  can  be  generated  within  the  ionization  zone  upstream  of  the  channel  exit.  However,  in 
practice,  this  metric  of  performance  in  a  typical  Hall  thruster  can  range  from  ~  1-10.  This  is  in  part 
due  to  the  loss  of  electrons  and  ions  through  recombination  reactions  on  the  channel  wall,  and  also 
due  to  the  anomalous  migration  of  electrons  from  the  exit  plane  through  the  region  of  strong 
magnetic  fields.  It  is  believed  that  this  poor  confinement  of  electrons  is  in  part  due  to  turbulent 
fluctuations  in  plasma  properties,  which  give  rise  Bohm-like  electron  transport  in  a  direction 
orthogonal  to  the  applied  magnetic  field.  The  Stanford  Research  was  aimed  at  understanding  this 
transport,  and  any  affects  that  wall  scattering  may  have  on  also  enhancing  electron  current  flow 
within  the  discharge  channel.  The  broader  long-term  objective  is  to  develop  passive  and/or  active 
control  strategies  to  enhance  or  suppress  electron  flow  in  regions  of  the  channel  where  such 
actuation  is  desirable.  Feedback  -  controlled  actuation  methods  to  reduce  plasma  fluctuations  where 
only  marginally  successful  back  in  the  mid  70 ’s  in  early  versions  of  Hall  thrusters  [1,2].  Using 
modem  diagnostic  and  characterization  strategies  to  obtain  an  improved  understanding  of  the  role 
that  fluctuations  and  wall  scattering  plays  in  electron  transport,  it  is  anticipated  that  we  can  design 
more  effective  control  and  actuation  strategies.  Such  strategies  will  require,  however,  the 
development  of  robust  models  for  electron  flow,  and  the  incorporation  of  such  models  into 
simulators  that  are  able  to  accurately  capture  thruster  response. 

Accomplishments  The  research  funded  under  this  grant  has  been  extremely  successful,  and  over  the 
course  of  the  two  years  and  four  months  of  funding,  has  led  to  12  archival  publications,  24 
conference  publications  or  abstracts,  and  2  invention  disclosures.  Here  we  provide  a  brief  description 
of  the  results,  and  direct  the  reader  to  the  extensive  publication  list,  where  the  findings  are  described 
in  considerable  detail.  In  some  cases,  papers  are  provided  as  attachments  in  the  Appendix. 

Using  extensive  measurements  of  the  time-averaged  plasma  properties  in  the  discharge  channel, 
including  laser-induced  fluorescence  measurements  of  neutral  and  ionized  xenon  velocities  [3,  see 
also  Publications  8,  12,  18,  and  21,  and  Appendix  I]  and  various  electrostatic  probe  diagnostics  [4, 
see  also  Publications  7,  12,  15,  16,  23,  26,  31,  32,  and  Appendix  II]  we  were  able  to  determine  the 
effective  Hall  parameter,  (<a«r)eff,  the  inverse  of  which  is  a  measure  of  the  crossed-field  electron 
mobility,  when  the  electron  cyclotron  frequency  is  much  greater  than  the  electron  collision 
frequency.  These  values  are  then  compared  to  the  classical  (collision-driven)  Hall  parameters 
expected  for  a  quiescent,  magnetized  plasma  (see  Fig.  1).  The  results  indicate  that  in  the  vicinity  of 
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FIG.  1.  Comparison  of  additive  collision  frequency  model  to 
experimental  and  classical  inverse  Hall  parameter  profiles  at 
various  axial  locations  within  the  discharge  channel.  The  exit 
plane  is  at  an  axial  location  of  x  =  0  mm,  and  the  anode  is  at  x  = 
-80  mm 


the  anode,  where  there  are  fewer  plasma 
instabilities,  the  electron  transport  mechanism  is 
found  to  be  in  agreement  with  a  “classical” 
model  based  on  elastic  electron  collisions  with 
the  background  neutral  xenon.  However,  we  find 
that  in  the  vicinity  of  the  discharge  channel  exit, 
where  the  magnetic  field  is  the  strongest  and 
where  there  are  intense  fluctuations  in  the 
plasma  properties,  the  inferred  Hall  parameter 
departs  from  the  classical  value.  Instead,  it  is 
found  to  be  closer  to  the  Bohm  value  of 
((ocer)cft  =16  [5],  These  results  provided  the  first 
direct  support  for  a  simple  model  for  the  Hall 
parameter  (and  hence  electron  mobility  or 
conductivity)  that  is  based  on  the  scalar  addition 
of  the  electron  collision  frequencies;  a  model 
that  was  proposed  by  Boeuf  and  Garrigues  [6]: 


/total 


GO  T 

V  ce  J  oscillations 


(1) 


The  solid  line  in  Fig.  1 .  is  the  result  of  this  model,  where  the  effects  of  wall  scattering  is  neglected, 
and  where  the  contribution  to  the  total  inverse  Hall  parameter  due  to  fluctuations  is  modeled  simply 
as  1/16  in  accordance  with  the  Bohm  theory  [5].  While  the  exact  mechanism  for  this  anomalous 
electron  mobility  is  still  the  subject  of  intense  research,  the  results  that  we  have  obtained  recently, 
draw  attention  to  the  possible  role  that  fluctuations  have  on  enhancing  electron  transport  in  regions 
where  the  electrons  are  highly  magnetized.  The  compelling  evidence  for  the  relationship  between 
plasma  fluctuations  and  anomalous  electron  transport  is  the  correlation  between  the  axial  location  of 
this  enhanced  transport,  and  the  location  where  intense  plasma  density  fluctuations  are  seen  to  occur. 
Using  Langmuir  probes  (operating  in  saturated  ion-collection  mode)  inserted  into  the  plasma  channel 
to  map  out  time  dependent  plasma  density  fluctuations  at  frequencies  f  <  400  kHz,  we  have 

extensively  characterized  the  dynamical 
behavior  of  plasma  density  within  the  discharge 
channel  [7,  see  also  Publications  3,  11,  13,  14, 
25,  29,  33,  and  Appendix  III].  An  example  of  a 
spectral  map  (frequency-position)  depicting  the 
local  fluctuations  in  plasma  density  within  the 
discharge  channel  is  shown  in  Fig.  2. 

A  comparison  (conditional  cross-spectral 
analysis)  of  the  signal  from  two  Langmuir 
probes,  located  at  the  same  axial  position,  but 
separated  by  some  angle  on  the  azimuth,  was 
used  to  measure  the  azimuthal  component  of  the 
phase  velocity  (wave  dispersion).  Fig.  3  shows 
a  typical  dispersion  plot,  rich  in  features,  clearly 
showing  the  presence  of  four  characteristic 


-40  -30  -20  -10  0  10 


Probe  Position  Along  Thruster  Channel 

Fig.  2.  Spectral  map  of  plasma  density  fluctuations  in  the 
channel  of  a  Hall  thruster.  The  maximum  color  scale  represents 
fluctuations  in  density  that  are  comparable  to  the  mean  value. 
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modes,  the  most  prominent  of  which  is 
that  of  a  pure  azimuthally  propagating 
disturbance  (m=l)  with  frequencies  in  the 
5-20  kHz  range,  and  of  phase  velocity  Vp  ~ 
1000  -  4000  m/s.  The  dotted  line  depicts 
the  expected  dispersion  for  a  pure 
azimuthal  mode  and  a  15°  tilt  angle  with 
respect  to  the  azimuthal  plane,  consistent 
with  the  observation  of  “ionization 
spokes”  seen  in  earlier  Hall  discharge 
devices  [8].  A  second  mode  is  that 
associated  with  axial  disturbances  at 
higher  frequencies,  the  strongest  of  which 
are  identified  on  the  figure  by 
superimposing  dashed-dot  vertical  lines. 
The  third  mode  that  is  identifiable  in  Fig. 
2.  is  that  of  an  m  =  12  azimuthal  wave 
(guided  by  the  theoretical  dashed  line), 
which  has  a  wavelength  equal  to  that  of 
the  separation  between  the  two  probes,  and 
is  believed  to  be  activated  by  the 
placement  of  the  probes  within  the  discharge  annulus  (30°  separation  results  in  a  wavelength  of 
1/1 2th  the  channel  circumference).  Finally,  the  fourth  mode  identified  by  the  solid  (m=4)  line  appears 
actually  as  a  lower  velocity  limit  to  the  high  frequency  disturbances.  The  appearance  of  this 
demarcation  is  attributed  to  an  interaction  between  the  axial  disturbances  and  the  four-fold  azimuthal 
asymmetry  imposed  by  the  four  equally  spaced  magnetic  solenoids  connected  to  the  outer  pole-piece 
[7],  It  is  noteworthy  that  the  dispersion  characteristics  of  the  disturbances  are  seen  to  vary 
dramatically  with  changes  in  operating  regimes  (e.g.,  applied  voltage),  or  changes  in  position  within 
the  discharge  channel  [7]. 

Under  this  grant,  we  have  complemented  these  experimental  studies  with  theoretical  work  on  the 
excitation  and  propagation  of  theses  disturbances,  including  the  determination  of  linear  growth  rates, 
so  as  to  better  understand  the  physical  nature  of  these  azimuthal  and  axial  waves,  and  the 
mechanisms  responsible  for  their  growth  [8,  see  also  Publications  10,  35,  and  Appendix  IV],  Using 
the  extensive  experimental  data  on  the  time-averaged,  spatially  varying  plasma  properties  within  our 
discharge,  we  have  performed  a  theoretical  study  of  the  response  of  this  plasma  to  small  (linear) 
perturbations  in  its  properties.  As  a  starting  point  for  this  analysis,  we  assume  a  two-dimensional 
fluid  description  that  includes  a  simplified  equation  for  the  electron  energy,  and  constrain  the 
azimuthal  wave-vector  such  that  we  excite  only  the  dominant  (m=l)  azimuthal  modes.  The  growth 
rate  and  frequencies  of  predominantly  axial  and  azimuthally  propagating  plasma  disturbances  are 
obtained  by  numerical  solution  of  the  resulting  eigenvalue  problem  under  a  quasi-uniform  plasma 
condition,  along  the  entire  discharge  channel.  The  results  identify  the  persistence  of  a  relatively  low 
frequency  instability  that  is  concentrated  largely  in  the  vicinity  of  the  exit  plane,  where  the  magnetic 
field  is  at  its  maximum  value,  consistent  with  experimental  observations  for  relatively  low  operating 
voltages  (~100V). 

The  results  of  a  typical  calculation  are  shown  in  Fig.  4.  Here  we  display  the  determined 
frequencies  of  the  unstable  wave  versus  axial  position,  for  various  propagation  angles.  It  is  apparent 


FIG.  3.  Wave  dispersion  plots  obtained  from  two  azimuthal  probes 
separated  by  30°,  for  low  discharge  voltage  conditions  (100  V),  and  at  x 
=  0  mm.  The  lines  shown  are  the  expected  theoretical  dispersion  lines  for 
azimuthal  waves  with  mode  numbers:  an  apparent  m  =  0  mode  (vertical 
dot-dashed  lines),  m  =  1  (dotted  line),  m  =  4  (sold  line),  and  m  =  12 
(dashed  line)  modes. 
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from  the  figure  that  the  disturbances  at  higher 
frequencies  favor  longitudinal  (axial)  propagation, 
whereas  azimuthal  propagation  is  largely 
associated  with  frequencies  below  20  kHz,  in 
qualitative  agreement  with  what  is  seen 
experimentally  (Fig.  3).  An  examination  of  the 
amplitudes  also  indicate  that  the  disturbances  are 
concentrated  between  the  axial  locations  x  =  -20 
to  10  mm,  also  in  agreement  with  the 
experimentally  measured  fluctuations  as  seen  in 
Fig.  2.  The  lower  frequency  azimuthal  waves 
extend  throughout  a  greater  region  of  the  channel 
whereas  the  higher  frequency  waves  are 
concentrated  mostly  in  the  region  between  x  =  -10 
mm  and  x  =  0.  The  sensitivity  of  the  behavior  of 
these  instabilities  to  various  plasma  properties 
indicate  that  the  disturbances  at  low  discharge  voltages  are  ionization  waves,  influenced  in  part  by 
the  spatial  variation  in  other  plasma  properties  (e.g.,  magnetic  field,  ion  velocity).  The  development 
of  an  analytical  model  for  these  low  frequency  disturbances,  and  an  extension  of  the  numerical 
simulations  to  better  understand  the  dynamical  behavior  of  these  discharges  at  high  operating 
voltages  is  presently  underway. 


Propagation  Angl*  [Dag] 


Axial  Position  [mm] 


FIG.  4.  Predicted  frequencies  for  unstable  (exponentially 
growing)  disturbances  along  the  discharge  channel,  for  various 
propagation  angles.  The  background  plasma  conditions  are 
those  measured  for  100V  discharge  operation. 


Finally,  it  has  been  an  objective  of  this  program,  to  understand  the  physics  of  these  devices,  so  as  to 
be  able  to  scale  the  thrusters  down  to  very  low  power  levels.  The  operation  of  E  X  B  discharge 
accelerators  at  these  low  power  ranges  introduce  many  engineering  challenges.  These  include 
increased  heat  loads  imposed  on  the  structure  due  to  the  nature  of  the  scaling  of  these  plasma 
devices,  and  the  increased  frequency  (possibly  into  the  MHz  regime,  also  due  to  the  nature  of  the 


FIG.  5.  Photograph  of  linear  thruster  operating  at  about 
75W  of  discharge  power. 


scaling)  of  the  strong  “loop”  instabilities  and  its 
impact  on  the  plasma,  and  power  processing  unit, 
amongst  others.  At  Stanford,  we  have  developed 
and  tested  a  linear-geometry,  open  electron  drift 
Hall  discharge.  The  linear  geometry  makes  it 
feasible  to  use  permanent  magnets  instead  of 
electromagnets  (this  is  important  if  micro¬ 
fabrication  methods  are  to  be  introduced).  The 
impetus  for  a  linear  geometry  was  largely  the 
result  of  the  experiments  described  above  on  co¬ 
axial  geometries,  which  indicate  that  the  drifting 
electrons  may  not  make  a  complete  revolution 
around  the  thruster  before  traversing  the  region  of 
strong  magnetic  fields. 

A  photograph  of  the  linear-geometry  thruster  that 
we  have  fabricated  and  tested  while  operating  at 
about  75W  of  power  is  shown  in  Fig.  5.  Details 
of  the  geometry  and  operating  characteristics  have 
been  presented  in  papers  (see  Publications  6,  17, 
and  Appendix  V],  although  we  have  not  yet 
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published  any  details  on  the  performance,  which 
is  quite  encouraging.  Using  facilities  down  at 
Edwards  Air  Force  Base,  we  were  able  to 
characterize  the  performance  of  this  linear 
thruster  over  a  wide  range  of  conditions.  An 
example  of  the  measured  thrust  efficiency  is 
shown  in  Fig.  6.  It  is  apparent  that  the  thrust 
efficiency  shows  a  marked  increase  with 
decreasing  power,  unlike  the  behavior  seen  in 
coaxial  devices.  At  the  low  power  operating 
point,  where  the  efficiency  exceeded  14%,  the 
measured  specific  impulse  approached  1100 
seconds.  As  far  as  we  know,  no  other  ultra-low 
power  Hall  discharge  has  performed  so  well  at 
these  conditions,  despite  the  lack  of  design 
optimization.  In  the  continuation  grant,  we 
propose  to  further  examine  the  performance  of 
Hall  thrusters  of  radically  altered  geometries, 
including  co-axial  designs  without  inner  magnetic 

cores,  such  as  the  designs  being  explored  at  the  Princeton  Plasma  Physics  Laboratory  by  Fisch, 

Raitses,  and  colleagues. 
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Laser-Induced  Fluorescence 


Measurements  of  Velocity  within  a  Hall  Discharge 

W.A.  Hargus,  Jr*  and  M.A.  Cappelli 
Mechanical  Engineering  Department 
Thermosciences  Division 
Stanford  University 
Stanford,  CA  94305 

Abstract 

The  results  of  a  study  of  laser  induced  fluorescence  velocimetry  of  neutral  and  singly  ionized  xenon 
in  the  plume  and  interior  portions  of  the  acceleration  channel  of  a  Hall  thruster  plasma  discharge  operating  at 
powers  ranging  from  250  to  725  W  are  described.  Axial  ion  and  neutral  velocity  profiles  for  four  discharge 
voltage  conditions  (100  V,  160  V,  200  V,  250  V)  are  measured  as  are  radial  ion  velocity  profiles  in  the  near 
field  plume.  Ion  velocity  measurements  of  axial  velocity  both  inside  and  outside  the  thruster  as  well  as  radial 
velocity  measurements  outside  the  thruster  are  performed  using  laser  induced  fluorescence  with  nonresonant 
signal  detection  on  the  xenon  ion  5d[4]7/2-6p[3]5/2  excitation  transition  while  monitoring  the  signal  from 
the  6s[2]3/2-6p[3]5/2  transition.  Neutral  axial  velocity  measurements  are  similarly  performed  in  the  inte¬ 
rior  of  the  Hall  thruster  using  the  6s[3/2]2  -6p[3/2]2  transition  with  resonance  fluorescence  collection. 
Optical  access  to  the  interior  of  the  Hall  thruster  is  provided  by  a  1  mm  wide  axial  slot  in  the  insulator  outer 
wall.  While  the  majority  of  the  ion  velocity  measurements  used  partially  saturated  fluorescence  to  improve 
the  signal  to  noise  ratio,  one  radial  trace  of  the  ion  transition  was  taken  in  the  linear  fluorescence  region  and 
yields  an  xenon  ion  translational  temperature  between  400  and  800  K  at  a  location  13  mm  into  the  plume. 

Introduction 

Due  to  their  high  specific  impulse  and  high  thrust  efficiencies,  Hall  thrusters  are  now  being  consid¬ 
ered  for  use  on  commercial,  research,  and  military  spacecraft.  This  technology  provides  economic  advan¬ 
tages  for  a  number  of  missions  and  its  use  can  be  translated  into  lower  launch  mass,  longer  time  on  station, 
or  larger  payloads  [1].  In  order  to  extend  the  performance  and  operating  envelope  of  Hall  thrusters,  there  is  a 
need  for  increased  understanding  of  the  complex  phenomena  that  control  propellant  ionization  and  accelera¬ 
tion  within  the  discharge.  In  order  to  more  fully  understand  the  physics  in  these  discharges,  several  labora¬ 
tory  model  Hall  thrusters  have  been  constructed  at  Stanford  University.  These  thrusters  serve  as  test  articles 
for  model  development  and  advanced  plasma  diagnostics  including  laser  induced  fluorescence  (LIF),  probes 
of  various  types,  and  thrust  measurements  [2-5]. 
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Laser  based  techniques  have  been  developed  to  nonintrusively  probe  neutral  and  ionized  xenon 
[2,5].  Such  measurements  in  the  plumes  of  Hall  and  other  types  of  ion  thrusters  provide  important  informa¬ 
tion  on  the  expansion  of  the  plasma  plume  and  its  potential  effect  on  satellite  propulsion  design  and  integra¬ 
tion.  Laser  based  diagnostic  measurements  have  been  previously  employed  to  examine  plume  plasma 
properties  in  other  electric  propulsion  devices.  For  example,  the  hydrogen  arcjet  has  been  extensively  stud¬ 
ied  using  lasers  to  measure  velocity,  temperature,  and  electron  number  density  [6].  However,  few  measure¬ 
ments  (especially  optical  measurements)  have  been  made  within  the  interior  discharge  of  these  devices.  The 
high  spatial  resolution  of  single  point  LIF  is  essential  in  probing  nonuniform  plasma  environments  such  as 
those  inside  Hall  thrusters  and  other  electric  propulsion  devices. 

Ion  velocity  measurements  of  axial  velocity,  both  inside  and  outside  the  thruster  as  well  as  radial 
velocity  measurements  outside  the  thruster,  are  performed  using  LIF  with  nonresonant  signal  detection 
exciting  the  xenon  ion  5d[4]7/2-6p[3]5/2  electronic  transition  while  monitoring  signal  from  the 
6s[2]3/2  -  6p[3]5/2  transition.  Neutral  axial  velocity  measurements  are  similarly  performed  in  the  interior  of 
the  Hall  thruster  using  the  6s[3/2]°  -6p[3/2]2  transition  with  resonance  fluorescence  collection.  The 
majority  of  the  velocity  measurements  use  partially  saturated  fluorescence  to  improve  the  signal  to  noise 
ratio. 


Theory 

Laser  Induced  Fluorescence 

Laser  induced  flourescence  is  used  to  detect  velocity-induced  shifts  in  the  spectral  absortpion  of 
xenon  atoms  and  ions.  The  fluorescence  is  monitored  as  a  continuous-wave  laser  is  tuned  in  frequency  over 
the  transition  of  interest,  of  energy  hv12 .  Measurements  can  be  made  with  high  spatial  resolution,  deter¬ 
mined  by  the  intersection  of  the  probe  laser  beam  with  the  fluorescence  optical  collection  volume. 

In  laser-induced  fluorescence  velocimetry,  an  atom  or  ion  moving  with  a  velocity  component  u  rel¬ 
ative  to  the  direction  of  the  incoming  laser  will  absorb  the  light  at  a  frequency  shifted  from  that  of  stationary 
absorbers  due  to  the  Doppler  effect.  The  magnitude  of  this  frequency  shift  8v  12  is: 


8v12  =  v 


u 

12- 


(1) 


The  measured  fluorescence  signal  is  given  by  [7-9]: 

Sf  =  Tldachv12A21N2  (2) 

where  pd  is  the  efficiency  of  the  detection  system,  ac  accounts  for  geometric  factors  involving  the  collec¬ 
tion  system,  and  A-,,  is  the  Einstein  coefficient  for  spontaneous  emission  of  the  relevant  transition.  For  low 
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laser  intensities,  rate  equation  analysis  indicates  that  the  upper  level  population  N2 ,  and  therefore  the  fluo¬ 
rescence  signal,  is  linearly  dependent  on  laser  intensity  at  steady  state,  i.e., 

N2~IvB,A  (3) 

where  lv  is  the  spectral  irradiance  at  frequency  v  ,  B12  is  the  Einstein  stimulated  absorption  coefficient,  and 
<|)v  is  the  transition’s  spectral  line  shape  which  accounts  for  the  variation  of  the  absorption  or  laser  excitation 
with  frequency.  The  line  shape  is  determined  by  the  environment  of  the  absorbing  atoms,  so  an  accurate 
measurement  of  the  line  shape  function  can  lead  to  the  determination  of  various  plasma  parameters.  How¬ 
ever,  for  velocity  measurements,  partially  saturated  fluorescence  with  a  distorted  line  shape  can  still  provide 
an  adequate  measure  of  the  mean  velocity.  This  was  experimentally  verified  with  several  saturation  studies 
which  also  examine  the  Doppler  shift.  The  variation  of  the  resultant  velocities  was  found,  in  this  study,  to  be 
less  than  the  experimental  uncertainty  for  the  ions  (±500  m/s),  or  for  the  neutrals  (±60  m/s). 

Several  factors  affect  the  line  shape  and  give  rise  to  broadening  and/or  shift  of  the  spectral  line.  In 
high  temperature  plasmas,  the  most  significant  is  Doppler  broadening  due  to  the  absorber’s  random  thermal 
motion,  characterized  by  the  molecular  kinetic  temperature,  Tkin .  When  the  absorbing  species  velocity  dis¬ 
tribution  is  Maxwellian  in  shape,  the  Doppler  broadening  results  in  a  Gaussian  line  shape.  Collisional  inter¬ 
actions  between  the  absorbers  and  other  species  in  the  plasma  give  rise  to  spectral  line  shapes  that  are  often 
Lorentzian.  This  includes  interactions  with  charged  particles  (Stark  broadening)  and  uncharged  particles 
(van  der  Waals  broadening).  If  both  Doppler  broadening  and  collisional  broadening  are  important  and  inde¬ 
pendent,  the  resulting  line  shape  is  a  convolution  of  the  Gaussian  and  Lorentzian  line  shape  into  a  Voigt  line 
shape  [9]. 

The  absorption  line  shape  is  an  intrinsic  property  of  the  absorbers,  whereas  the  fluorescence  excita¬ 
tion  line  shape  is  the  variation  in  the  detected  fluorescence  signal  with  frequency  as  the  laser  is  tuned  across 
the  absorption  line  feature.  If  the  laser  excitation  significantly  perturbs  the  populations  of  the  coupled  levels, 
it  is  said  to  be  saturating  the  transition  and  the  fluorescence  signal  is  then  a  nonlinear  function  of  laser  inten¬ 
sity.  In  cases  where  the  laser  intensity  is  significantly  below  the  saturation  level  and  the  laser  linewidth  is 
small  compared  to  the  measured  linewidth,  the  fluorescence  excitation  line  shape  reflects  the  spectral 
absorption  line  shape  as  given  by  Eqns.  2  and  3.  When  the  laser  intensity  is  sufficiently  high  to  saturate  the 
transition,  the  fluorescence  excitation  line  shape  is  broader  than  the  spectral  line  shape  and  the  fluorescence 
intensity  is  less  than  it  would  be  if  it  were  linear  with  the  laser  intensity  lv .  The  saturation  intensity,  defined 
as  that  intensity  which  produces  a  fluorescence  signal  half  of  what  it  would  be  if  the  fluorescence  was  linear 
with  Iv ,  depends  inversely  on  the  line  strength  of  the  particular  transition.  Stronger  transitions  have  a 
smaller  saturation  intensity  and  thus  a  larger  saturation  effect  for  a  given  laser  intensity. 

Hyperfine  Structure 
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The  nine  isotopes  of  xenon,  the  propellant  most  commonly  used  in  Hall  plasma  thrusters,  each  have 
a  slight  difference  in  their  electron  transition  energies  due  to  their  differences  in  mass  [10].  The  odd  mass 
isotopes  are  further  spin  split  due  to  nuclear  magnetic  dipole  and  electric  quadrapole  moments.  Nuclei  which 
have  an  odd  number  of  protons  and/or  an  odd  number  of  neutrons  posses  an  intrinsic  nuclear  spin  Ih/2tt , 
where  I  is  integral  or  half-integral  depending  on  if  the  atomic  mass  is  even  or  odd,  respectively  [11].  For 
nuclei  with  non-zero  nuclear  spin  (angular  momentum),  the  interaction  of  the  nucleus  with  the  bound  elec¬ 
trons  lead  to  the  splitting  of  levels  with  angular  momentum  J  into  a  number  of  components,  each  corre¬ 
sponding  to  a  specific  value  of  the  total  angular  momentum  F  =  I  +  J  [12].  As  a  result  of  this  interaction,  F 
is  a  conserved  quantity  while  I  and  J  individually  are  not.  The  interaction  is  weak  allowing  the  hyperfine 
splitting  of  each  level  to  be  taken  independently  of  the  other  levels.  The  number  of  nuclear  spin  split  hyper- 
fine  components  is  21  +  1  if  J  >  I  and  2J  +  1  if  J  <  I ,  with  F  taking  on  the  values  F  =  J  + 1,  J  + 1  - 1, ....  |J  -I| 
[12,13]  while  satisfying  the  selection  rules  imposed  on  F,  i.e.,  AF  =  0,  ±1 ,  unless  F  =  0,  in  which  case 
AF*0. 


With  these  selection  rules  on  the  quantum  numbers  for  a  particular  electronic  transition,  and  with 
knowledge  of  the  hyperfine  structure  constants  which  characterize  the  magnetic  dipole  and  electric  quadru¬ 
ple  moments  of  the  nucleas  [11],  the  hyperfine  energy  shifts  from  the  position  of  the  energy  for  the 
unshifted  level  with  angular  momentum  J  can  be  easily  calculated  [13].  The  relative  intesities  of  transitions 
between  these  levels  is  derived  assuming  Russel-Saunders  coupling  [13],  allowing  the  complete  construc¬ 
tion  of  the  fluorescence  lineshape.  Of  course,  the  intensities  of  the  isotope  shifted  transitions  are  propor¬ 
tional  to  each  isotope’s  relative  abundance  [11,14-15].  The  relative  intensities  of  the  nuclear  spin  split 
hyperfine  splitting  are  governed  by  two  summation  rules  [12].  First,  the  sum  of  the  intensities  of  all  the  lines 
of  the  hyperfine  structure  of  a  transition  J  — » J'  (the  prime  refers  to  the  upper  level  involved  in  the  transition) 
originating  from  a  component  F  of  the  level  J  is  proportional  to  the  statistical  weight  of  this  component, 
2F  +  1  .  Second,  the  sum  of  the  intensities  of  all  the  lines  of  the  hyperfine  structure  the  transition  J  — >  J'  end¬ 
ing  on  the  component  F'  of  the  level  J'  is  proportional  to  the  statistical  weight  of  this  component,  (2F  +  1) . 
With  these  two  sum  rules,  a  system  of  linear  equations  are  solved  for  the  relative  intensities  of  the  nuclear 
spin  split  components  of  each  isotope. 

Xenon  Spectroscopy 

For  the  results  reported  here,  we  probed  the  823.2  nm  and  834.7  nm  electronic  transitions  of  neutral 
and  singly  ionized  xenon,  respectively.  The  isotopic  and  nuclear-spin  effects  contributing  to  the  hyperfine 
structure  of  the  6s[3/2]°  -6p[3/2]2  neutral  xenon  transition  at  823.2  nm  produce  21  spectral  lines,  15  of 
which  are  shown  in  the  transition  schematic  in  Fig.  la  (the  other  six  are  of  the  same  transition  shown  for  the 
even  isotopes).  Similarly,  the  5d[4]7/2  -  6p[3]5/2  xenon  ion  transition  at  834.7  nm  has  a  total  of  19  isotopic 
and  spin  split  components,  13  of  which  are  shown  in  Fig.  lb  (again,  the  other  six  are  of  the  same  transition 
shown  for  the  even  isotopes).  The  hyperfine  splitting  constants  are  only  known  for  a  limited  set  of  energy 
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levels  [16-20].  While  isotope  shifts  and  all  of  the  hyperfine  splitting  constants  are  available  for  the  upper 
and  lower  levels  of  the  neutral  823.2  nm  transition,  the  834.7  nm  xenon  ion  transition  only  has  data  on  the 
nuclear  spin  splitting  constants  of  the  upper  state,  and  no  information  is  available  on  the  transition  dependent 
isotope  shifts.  For  LIF  measurements  primarily  aimed  at  determining  velocities  within  the  plasma  flow,  it  is 
often  convenient  to  probe  more  accessible  transitions  for  which  there  is  incomplete  knowledge  of  the  isoto¬ 
pic  and  nuclear  spin  splitting  constants.  Manzella  has  shown  that  the  5d[4]7/2-6p[3]5/2  xenon  ion  transi¬ 
tion  at  834.7  nm  may  be  used  to  make  velocity  measurements  in  a  Hall  thruster  plume  [21].  A  convenient 
feature  of  this  transition  is  the  presence  of  a  relatively  strong  line  originating  from  the  same  upper  state 
(6s[2]3/2- 6p[3]5/2  transition  at  541.9  nm)  which  allows  for  nonresonant  fluorescence  collection  [22].  A 
nonresonant  fluorescence  scheme  is  preferred  where  there  is  the  possibility  of  laser  scattering  from  internal 
surfaces  of  the  discharge. 


Experiment 

Test  Facility 

The  test  facility  consists  of  a  non-magnetic  cylindrical  stainless  steel  vacuum  vessel  1  m  in  diame¬ 
ter  and  1.5  m  in  length,  pumped  by  two  50  cm  diffusion  pumps  (backed  by  a  425  1/s  mechanical  pump) 
mounted  at  each  end.  Figure  2  shows  a  schematic  of  the  facility.  The  base  pressure  with  no  flow  is  approxi¬ 
mately  10"6  Torr  as  measured  by  an  ionization  gauge  calibrated  for  nitrogen.  Chamber  pressures  during  dis¬ 
charge  operation  result  in  pressures  of  approximately  10'4  Torr  (3xl0"5  Torr  on  xenon,  ±50%).  Gas  flow  to 
the  Hall  thruster  anode  and  cathode  is  throttled  by  two  Unit  Instruments  1200  series  mass  flow  controllers 
factory  calibrated  for  xenon.  The  gas  (propellant)  used  in  this  study  was  research  grade  (99.995%)  xenon. 

The  thruster  is  mounted  on  a  two  axis  translation  system  having  a  vertical  range  of  travel  of  30  cm 
and  axial  (horizontal)  range  of  8  cm,  each  with  a  resolution  of  about  10  /tm.  Although  the  repeatability  is 
considerably  courser,  it  is  still  significantly  less  than  the  dimensions  of  the  laser  probe  volume  in  the  axial 
dimension.  For  internal  LIF  measurements,  the  thruster  is  mounted  on  a  platform  that  provides  optical 
access  through  a  slot  machined  into  the  acceleration  channel  wall  (described  below).  The  view  of  the 
thruster  mounting  scheme  is  also  shown  in  Fig.  2. 

Hall  Thruster 

The  thruster  used  in  this  study  has  been  described  in  detail  extensively  elsewhere  [4].  It  has  four 
outer  magnetic  windings  and  a  single  inner  winding  generating  the  magnetic  field,  which  is  a  maximum  near 
the  discharge  exit  plane.  The  insulator  is  constructed  from  2  sections  of  cast  alumina  tubing  84  mm  in  length 
cemented  to  a  machinable  alumina  back  plate.  The  anode,  which  also  serves  as  the  gas  distributer,  is  posi¬ 
tioned  at  the  back  of  the  discharge  channel  constructed  from  the  coaxial  assembly  of  the  alumina  tubes.  The 
channel  has  an  outer  diameter  of  95  mm  and  a  12  mm  width.  A  cross-sectional  view  of  the  assembled  dis¬ 
charge  thruster  is  shown  in  Fig.  3.  Not  shown  in  Fig.  3  is  the  external  hollow  cathode  (ion  Tech  HCN-252), 
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positioned  with  its  orifice  20  mm  beyond  the  thruster  exit  plane,  used  to  neutralize  the  ion  beam  and  sustain 
the  anode  discharge. 

For  the  LIF  studies  within  the  discharge  chamber,  a  slot  parallel  to  the  axis  1  mm  wide  was 
machined  into  the  outer  insulator.  The  slot  width  was  selected  to  be  less  than  the  local  electron  Larmor 
radius  to  minimize  its  influence  on  the  discharge  properties.  The  slot  provides  optical  access  to  the  interior 
of  the  Hall  thruster,  and  operation  with  the  slot  did  not  differ  significantly  from  the  operation  with  the 
unslotted  insulator  (identical  current-voltage  characteristics).  Despite  the  presence  of  the  slot,  optical  access 
to  the  interior  is  still  limited  near  the  exit  plane  by  the  front  plate  of  the  magnetic  circuit,  precluding  optical 
measurements  very  near  the  exit  plane.  The  front  iron  plate  is  not  cut  since  this  would  modify  the  local  mag¬ 
netic  field. 

Laser  Induced  Fluorescence 

The  experimental  apparatus  used  for  the  LIF  measurements  includes  a  tunable  Coherent  899-21 
single  frequency  titanium  sapphire  laser,  actively  stabilized  to  provide  line  widths  on  the  order  of  1  MHz 
with  near  zero  frequency  drift.  The  titanium  sapphire  laser  is  pumped  by  a  Coherent  Nd:YAG  solid  state 
Verdi  pump  laser,  which  provides  5  W  of  single  mode  pump  power  at  532  nm.  The  tunable  laser  wavelength 
is  monitored  by  a  Burleigh  Instruments  WA-1000  scanning  Michaelson  interferometer  wavemeter  with  a 
resolution  of  0.01  cm'1. 

The  laser  beam  is  chopped  to  allow  for  phase  sensitive  detection  and  directed  into  the  Hall  thruster 
plume  by  a  series  of  mirrors.  For  axial  velocimetry  measurements,  the  slightly  divergent  beam  (1.7  milliradi- 
ans  full  angle)  is  focused  to  a  submillimeter  beam  waist  by  a  50  mm  diameter,  1.5  m  focal  length  lens  as 
shown  in  Fig.  2.  For  radial  velocimetry  measurements,  the  beam  enters  through  a  side  window  and  is 
focused  by  a  50  mm  diameter,  50  cm  focal  length  lens  as  shown  in  Fig.  5.  The  collection  optics  for  both 
radial  and  axial  velocity  measurements  consist  of  a  75  mm  diameter,  60  cm  focal  length,  collimating  lens. 
The  collected  light  is  then  focused  on  to  the  entrance  slit  of  a  0.5  m  Ebert-Fastie  monochromator  with  a  50 
mm  diameter,  30  cm  focal  length,  lens.  An  optical  field  stop  is  placed  between  the  two  lenses  to  match  the  F/ 
#  of  the  optical  train  with  that  of  the  monochromator.  The  monochromator  is  used  in  the  collection  optical 
train  as  a  narrow  band  optical  filter  so  that  only  light  from  the  transition  of  interest  is  collected.  With 
entrance  and  exit  slits  fully  open  (425  /tm),  the  600  groove/mm  plane  grating  (blazed  for  600  nm)  within  the 
monochromator  allows  the  exit  slit  mounted  Hamamatsu  R928  photomultiplier  tube  (PMT)  to  sample  a 
wavelength  interval  of  approximately  1  nm.  The  orientation  of  the  monochromator  allows  the  slit  aperture  to 
define  the  length  of  the  probe  beam  along  which  the  fluorescence  is  collected.  The  sample  volume  for  the 
axial  data  presented  in  this  work  is  approximately  100  ^m  in  diameter  and  2  mm  in  length.  Similarly,  the 
sample  volume  of  the  radial  data  is  approximately  100  jam  in  diameter  and  1  mm  in  length.  For  neutral 
xenon  LEF  velocimetry  measurements,  a  portion  of  the  probe  beam  is  split  from  the  main  beam,  passed 
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through  a  xenon  glow  discharge  tube,  and  used  as  an  stationary  absorption  reference.  A  silicon  photodiode 
monitors  the  absorption  signal.  This  use  of  the  glow  discharge  tube  is  only  possible  for  neutral  xenon.  The 
glow  discharge  does  not  support  a  sufficient  population  of  excited  state  ions.  However,  the  unshifted  ion 
transition  line  center  was  confirmed  to  correspond  to  that  taken  from  the  literature  by  the  observation  of  very 
slow  moving  ions  near  the  anode,  and  by  radial  measurements  near  the  acceleration  channel  center. 

The  probe  laser  beam  is  chopped,  and  the  LIF  signal  is  collected  using  a  Stanford  Research  Sys¬ 
tems  SRS-850  digital  lockin  amplifier.  Concurrently,  the  absorption  signal  from  the  stationary  reference  is 
collected  using  an  SRS-530  lockin  amplifier.  Data  from  the  absorption  signal,  laser  power  output,  and  the 
wavemeter  are  stored  on  the  digital  lockin  amplifier  using  3  available  analog  inputs  along  with  the  LIF  sig¬ 
nal.  Typical  tests  consist  of  a  12-20  GHz  scan  of  the  probe  laser  frequency  over  a  3  minute  period.  The  beam 
is  chopped  at  a  frequency  of  1.5  kHz,  and  both  lockin  amplifiers  use  1  s  time  constants.  Data  is  sampled  at  8 
Hz,  producing  four  traces  of  approximately  2,000  points  for  each  velocity  data  point.  Several  unsaturated 
traces  using  lower  laser  intensities,  10  s  time  constants,  and  10  min  scans  were  also  collected. 

Results  and  Analysis 

The  Hall  thruster  is  operated  at  four  conditions  for  the  work  reported  here.  At  each  condition,  the 
peak  magnetic  field  is  125  G,  the  mass  flow  to  the  anode  is  2  mg/s,  and  the  mass  flow  to  the  external  hollow 
cathode  is  0.3  mg/s.  The  test  conditions  correspond  to  discharge  voltages  of  100,  160,  200,  and  250  V.  The 
anode  currents  for  these  conditions  are  2.1,  2.4,  2.6,  and  2.9  A,  respectively.  The  total  power  consumed  by 
the  cathode  and  magnet  circuit  is  approximately  30  W.  It  should  be  noted  that  the  power  dissipated  in  the 
ballast  resistors  on  the  anode  and  cathode  keeper  lines  (~10  W)  is  not  included  in  these  calculations. 

All  spatially  resolved  measurements  are  referenced  to  a  two  coordinate  system  shown  in  Fig.  5.  The 
position  in  the  radial  coordinate  is  referenced  to  the  radial  location  corresponding  to  the  midpoint  of  the 
acceleration  channel  using  the  variable  D ,  defined  as  positive  toward  the  thruster  centerline.  The  axial  coor¬ 
dinate  is  given  by  Z  which  is  the  distance  from  the  thruster  exit  plane  and  is  defined  as  positive  along  the 
thrust  vector. 

Saturation  Study 

The  5d[4]7/2  —  6p[3]5/2  xenon  ion  transition  at  834.7  nm  was  probed  to  extract  local  velocity  data 
from  the  Doppler  shift  of  the  measured  fluorescence.  Figure  6  shows  the  saturation  curve  of  the  transition  at 
a  location  13  mm  from  the  exit  plane  at  the  center  of  the  acceleration  channel.  Here,  the  laser  beam  is  propa¬ 
gating  normal  to  the  thrust  vector  and  the  measurement  is  of  the  radial  velocity  component.  A  typical  satu¬ 
rated  trace  (laser  power  of  20  mW)  used  to  determine  the  velocity  in  the  probed  volume  is  compared  to  an 
unsaturated  trace  (laser  power  of  0.2  mW)  from  the  same  location  in  Fig.  7.  Both  traces  are  normalized  in 
the  figure  to  unity  peak  signal. 
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Velocity  measurements  of  the  saturation  curve  data  making  up  Fig.  6  are  shown  in  Fig.  8.  The  mean 
velocity  for  these  data  points  yields  a  value  of  -62  m/s  with  a  standard  deviation  of  65  m/s  and  a  range  of  127 
m/s.  The  absolute  accuracy  of  the  measurements  is  limited  by  the  specified  uncertainty  of  the  wavemeter 
reading,  which  is  ±500  m/s  at  this  wavelength.  The  important  conclusion  that  may  be  drawn  from  this  data  is 
that  the  5d[4]7/2-6p[3]5/2  xenon  ion  transition  of  xenon  provides  useful  LIF  velocimetry  data  even  when 
partially  saturated.  This  allows  the  collection  of  partially  saturated  fluorescence  signals,  maximizing  the  sig¬ 
nal  to  noise  ratio  and/or  allowing  for  faster  scans  while  still  being  able  to  extract  velocity  data  from  the  fluo¬ 
rescence  signal.  These  measurements  also  serve  as  a  confirmation  of  an  improved  accuracy  of  the  unshifted 
wavelength  since  a  finite  radial  velocity  is  not  expected  at  this  location.  Similar  results  were  also  obtained  in 
a  neutral  xenon  LIF  saturation  study. 

Xenon  Ion  LIF  Velocimetry 

The  axial  LIF  velocimetry  measurements  consist  of  two  data  sets  for  each  operating  condition.  The 
first  data  set  includes  ion  velocity  measurements  taken  externally  extending  from  the  exit  plane  to  approxi¬ 
mately  Z  =  35  mm.  The  second  data  set  contains  internal  axial  velocity  measurements  taken  from  the  exit 
plane  to  approximately  Z  =  -75  mm.  These  ranges  are  imposed  by  the  limited  travel  of  the  translation  stage 
that  provides  axial  motion  to  the  Hall  thruster.  Axial  velocimetry  measurements  were  taken  every  2.5  mm 
with  a  sample  probe  volume  of  100  /im  in  diameter  and  2  mm  long.  For  several  of  the  test  conditions  exam¬ 
ined,  profiles  of  the  axial  velocity  across  the  coordinate  D  were  also  examined.  Measurements  of  near  plume 
radial  ion  velocities  were  also  performed.  These  measurements  are  limited  to  the  plume  due  to  the  limited 
optical  access  to  the  interior  of  the  Hall  thruster. 

The  complete  axial  velocity  profiles  for  the  four  cases  examined  are  shown  in  Fig.  9.  The  error  bars 
correspond  to  the  ±500  m/s  uncertainty  associated  with  the  determination  of  the  magnitude  of  the  Doppler 
shift  relative  to  the  unshifted  line  center.  The  axial  velocity  profiles  exhibit  the  expected  behavior.  The 
velocity  is  near  zero  near  the  anode  (Z  =  -78  mm),  and  begins  to  rise  near  Z  =  -10  mm  at  the  edge  of  the 
acceleration  zone.  The  ions  are  rapidly  accelerated  in  the  region  of  the  exit  plane  and  reach  their  full  velocity 
in  the  neighborhood  of  Z  =  20  mm.  This  latter  position  corresponds  to  the  location  of  the  hollow  cathode 
neutralizer  relative  to  the  body  of  the  thruster  and  is  sometimes  referred  to  as  the  cathode  plane  in  the  litera¬ 
ture  [21]. 

The  length  of  the  acceleration  region  in  all  cases  is  invariant  at  30  mm.  Therefore,  increases  in  the 
anode  potential  result  in  a  concomitant  increase  of  the  axial  electric  field  component  within  the  thruster 
acceleration  channel.  The  initial  acceleration  is  seen  to  begin  at  10  mm  within  the  thruster  where  the  mag¬ 
netic  field  has  a  value  of  approximately  85%  of  the  centerline  maximum.  The  propellant  acceleration  is 
completed  20  mm  beyond  the  exit  plane  when  the  magnetic  field  has  a  value  of  approximately  25%  of  the 
centerline  maximum. 
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Significant  acceleration  of  the  ionized  propellant  occurs  outside  the  Hall  thruster.  This  is  consistent 
with  the  results  first  identified  in  a  Hall  discharge  with  a  much  shorter  acceleration  channel  [2].  It  is  note¬ 
worthy  that  the  velocity  increment  imparted  to  the  propellant  outside  the  of  Hall  thruster  is  essentially  con¬ 
stant,  invariant  to  discharge  voltage,  with  an  average  value  of  5,000  m/s.  Only  for  the  100  V  case,  does  the 
majority  of  the  acceleration  occur  outside  the  thruster.  Higher  discharge  voltages  appear  to  have  a  constant 
percentage  of  the  acceleration  occurring  externally.  It  is  also  equally  interesting  to  cast  this  result  in  terms  of 
the  kinetic  energy  of  the  propellant.  In  the  case  of  the  100  V  discharge  operation,  approximately  90%  of  the 
energy  is  deposited  into  the  propellant  between  the  exit  and  cathode  planes.  For  the  250  V  case,  the  fraction 
of  energy  deposition  beyond  the  exit  plane  is  nearer  to  60%.  It  appears  then  that  the  majority  of  the  energy 
deposition  into  the  Hall  thruster  propellant  occurs  outside  the  thruster.  However,  since  the  thrust  is  equiva¬ 
lent  to  the  momentum  flux,  the  majority  (65%)  of  the  thrust  is  still  generated  within  the  thruster  in  all  cases 
above  100  V.  It  is  further  noteworthy  that  in  all  the  cases  approximately  60  eV  is  unaccounted.  This  value  is 
constant  to  within  the  uncertainties  of  the  velocity  measurement  and  implies  that  the  mechanism  responsible 
for  this  loss  is  invariant  with  the  applied  anode  potential.  This  energy  loss  may  be  a  product  of  the  anode  and 
cathode  potential  falls,  or  other  mechanisms  common  to  these  discharges  (e.g.  distribution  of  in  the  ion  pro¬ 
duction,  or  resistive  losses  in  regions  of  low  ionization  fractions).  It  is  noted  that  these  measurements  are 
limited  to  axial  velocities  and  do  not  account  for  losses  due  to  plume  divergence. 

Several  radial  profiles  of  the  axial  ion  velocity  in  Fig.  10  illustrate  the  radial  variation  of  the  mea¬ 
sured  axial  velocities  for  a  discharge  voltage  of  160  V  at  two  locations  in  the  plume  and  one  within  the 
thruster  for  a  discharge  voltage  of  200  V.  The  width  of  the  acceleration  channel  is  12  mm  (-6  mm  <  D  <  6 
mm).  The  axial  velocity  profile  is  nearly  flat  which  strongly  implies  lines  of  constant  potential  in  the  radial 
direction.  The  slightly  concave  distribution  inside  the  channel  (200  V  case)  is  consistent  with  the  expected 
equipotential  surfaces  associated  with  the  radial  magnetic  field. 

Radial  velocity  measurements  were  also  performed  in  the  plume.  Combined  with  the  above  axial 
velocity  measurements  and  the  observation  that  the  axial  velocity  appears  to  be  independent  of  D,  vector 
plots  of  the  near  plume  may  be  constructed.  Figure  11  shows  a  vector  plot  constructed  for  a  200  V  discharge 
voltage.  In  this  case,  radial  velocities  vary  linearly  with  D  with  near  zero  velocity  at  D  =  0  mm  and  peak  at 
values  above  6,000  m/s  as  close  as  Z  =  13  mm  and  D  =  8  mm.  A  prominent  feature  of  a  Hall  thruster  plume 
is  the  luminous  central  cone  originating  at  the  thruster  axis  adjacent  to  the  center  pole  piece.  Although  the 
propellant  stream  exits  the  thruster  in  an  annulus,  an  intense,  optically  emitting  conical  feature  extends  a  sig¬ 
nificant  distance  into  the  vacuum  chamber  (10  to  30  cm  for  our  background  pressures)  and  is  especially  evi¬ 
dent  at  higher  discharge  voltages.  The  inward  focus  of  the  divergent  propellant  flow  is  likely  a  contributing 
mechanism  in  the  formation  of  this  intriguing  feature  which  is  not  precisely  understood. 

Ion  Line  Shape  Analysis 
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The  hyperfine  splitting  constants  for  the  xenon  ion  5d[4]7/2-6p[3]5/2  transition  at  834.7  nm  are 
only  known  for  the  upper  6p[3]5/2  state.  The  isotope  shifts  for  this  transition  are  also  unknown.  While  this 
does  not  greatly  impact  the  accuracy  to  which  velocities  may  be  determined,  it  does  impact  the  ability  to  use 
this  transition  for  accurate  measurements  of  the  ion  velocity  distribution,  or  temperature.  We  point  out  here 
that  the  notion  of  a  temperature  for  the  ions  in  this  flow  is  highly  tenuous,  as  the  ions  may  be  treated  kineti- 
cally,  do  not  suffer  many  collisions,  and  are  born  (created)  at  distributed  axial  locations  within  the  discharge. 
Furthermore,  a  time  averaged  broadening  of  the  fluorescence  excitation  spectra  can  arise  due  to  oscillations 
in  the  zone  of  ionization.  King  has  shown  by  mass  and  resolved  energy  analysis  that  the  axial  velocity  of  the 
ions  has  an  energy  distribution  of  approximately  10  eV  (in  the  far  field)  due  in  part  to  the  plasma  oscillations 
within  the  Hall  thruster  [23].  The  issue  of  the  distribution  of  ion  velocities  due  to  the  axial  location  where 
they  were  born  is  minimized  by  examining  the  fluorescence  spectra  in  the  radial  direction.  The  spectra  is 
taken  from  the  position  with  the  minimum  measured  velocity,  approximately  100  m/s  (±500  m/s)  at  a  loca¬ 
tion  of  D  =  0  mm  and  Z  =  13  mm.  In  order  for  a  temperature  to  be  extracted  from  the  measured  lineshape, 
the  ion  population  is  assumed  to  be  Maxwellian,  or  at  least  frozen  in  a  close  facsimile. 

An  estimate  of  the  ion  kinetic  temperature  with  uncertainties  of  40-70%  is  possible.  A  radial  excita¬ 
tion  unsaturated  fluorescence  trace  is  shown  in  Fig.  12  and  compared  to  a  line  shape  model  developed  by 
Cedolin  [24].  The  model  uses  the  5d[3]7/2  lower  level  hyperfine  spin  splitting  and  isotopic  shift  data  from 
the  605.1  nm  transition  and  the  measured  splitting  data  for  the  upper  6p[3]5/2  level.  Lorentzian  broadening 
is  neglected  and  only  Doppler  broadening  is  considered.  The  best  fit  of  this  model  predicts  a  kinetic  temper¬ 
ature  of  approximately  450  K.  The  model  does  not  completely  predict  the  outlying  features,  but  this  is 
expected  since  the  spin  splitting  constants  for  the  5d[3]7/2  state  are  used  for  the  lower  level  rather  than 
those  for  the  5d[4]7/2  state  probed  here.  If  hyperfine  splitting  is  ignored  and  only  the  isotope  shifts  corre¬ 
sponding  to  the  values  for  the  605.1  nm  transition  are  used,  the  model  predicts  a  kinetic  temperature  of 
approximately  800  K.  Neglecting  all  hyperfine  splitting  mechanisms,  including  isotopic  and  nuclear  spin 
splitting,  an  upper  bound  on  the  temperature  of  1,700  K  is  determined  from  the  fluorescence  Doppler  half 
width.  The  uncertainty  of  this  measurement  is  in  large  part  due  to  the  uncertainties  of  the  spectral  data  as 
well  as  due  to  the  noise  in  the  fluorescence  signal.  A  similar  measurement  in  the  plume  of  a  SPT-100  by 
Manzella  yielded  a  kinetic  temperature  of  approximately  800  K  [21].  It  should  be  noted  that  Manzella  used 
an  incorrect  value  of  J  for  the  lower  state  which  was  first  misidentified  by  Humphreys  [25]  and  propagated 
by  Moore  [26]  before  finally  being  corrected  by  Hansen  and  Persson  [22], 

Neutral  LIF  Velocimetry 

Figure  13  shows  the  results  of  axial  neutral  xenon  velocity  measurements  within  the  acceleration 
channel.  The  four  cases  examined  show  very  similar  behavior.  The  initial  velocity  near  the  anode  is  very 
low.  The  neutral  velocity  slowly  rises  until  a  position  of  approximately  20  mm  within  the  thruster.  At  this 
point,  where  the  ion  acceleration  is  also  seen  to  begin,  the  neutrals  are  accelerated  at  a  higher  rate  until  near 
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the  exit  plane  where  the  acceleration  appears  to  slow  and  even  possibly  reverse  when  the  thruster  is  operated 
at  higher  voltages.  The  decrease  in  neutral  xenon  velocity  is  likely  due  to  thruster  ingestion  of  background 
xenon.  Since  the  effect  appears  to  grow  with  increased  discharge  voltage,  it  is  possible  that  a  portion  of  the 
propellant  flow  reflected  from  nearby  vacuum  facility  walls  (now  entirely  consisting  of  neutrals)  is  ingested 
by  the  thruster. 

Flow  from  the  cathode  may  be  eliminated  as  the  source  of  the  apparent  deceleration  of  the  neutrals 
for  a  variety  of  reasons.  First,  the  flow  from  the  cathode,  although  15%  of  the  anode  flow,  is  exiting  from  a  2 
mm  orifice  approximately  12  cm  above  the  sample  volume.  Flow  from  the  cathode  should  be  sufficiently 
diffuse  and  not  affect  the  neutral  velocity  measurements  in  the  sample  volume.  Second,  the  cathode  is  2  cm 
downstream  of  the  exit  plane  and  angled  30°  from  the  front  plate  pointed  along  the  thrust  axis.  It  is  difficult 
to  envision  the  cathode  affecting  the  neutral  velocity  measurements  within  the  acceleration  channel.  It  is 
therefore  almost  certain  that  random  neutral  flux  from  the  chamber  background  is  responsible  for  the  appar¬ 
ent  drop  in  neutral  velocity  seen  near  the  exit  plane  in  Fig.  13. 

Due  to  the  highly  nonequilibrium  nature  of  the  Hall  thruster,  it  is  important  to  understand  the  appar¬ 
ent  acceleration  of  the  neutrals  beginning  40  mm  upstream  of  the  exit  plane.  The  plasma  within  the  Hall 
thruster  is  required  to  be  of  low  collisionality  by  the  constraint  that  the  magnetic  field  limit  the  electron  flux 
to  the  anode.  The  disparate  velocities  of  the  ions  and  neutrals  strongly  suggest  that  the  neutral  and  ion  popu¬ 
lations  are  not  coupled.  As  such,  the  apparent  acceleration  of  the  neutrals  may  actually  be  an  artifact  of  the 
time  of  flight  of  the  neutrals  through  the  volumetric  zone  of  ionization.  Slower  neutrals,  or  neutrals  that 
travel  a  longer  effective  path  length  due  to  collisions  with  the  walls  of  the  acceleration  channel,  have  a 
greater  probability  of  being  ionized  than  do  neutrals  in  the  high  energy  portion  of  the  velocity  distribution. 
Therefore,  neutrals  from  the  high  energy  range  of  the  velocity  distribution  are  more  likely  to  reach  the 
upstream  portion  of  the  acceleration  channel.  In  this  case,  there  is  no  actual  acceleration  of  the  neutrals,  but 
rather  a  depletion  of  the  slower  moving  neutrals  by  ionization. 

The  depletion  of  the  slower  neutral  velocity  classes  accounting  for  the  apparent  acceleration  of  the 
neutrals  as  seen  in  Fig.  13  may  be  qualitatively  explained  by  considering  the  one  dimensional  Boltzmann 
equation  [27]. 

uA[nf(u)]  =  -[nf(u)]neSj  (18) 

Where  n  is  the  neutral  number  density,  u  is  the  neutral  velocity  class,  z  is  the  spatial  coordinate,  f(u)  is  the 
velocity  distribution  function,  ne  is  the  local  plasma  electron  number  density,  and  Sj  is  the  ionization  rate 
coefficient.  The  generalized  one  dimensional  Boltzmann  equation  is  simplified  here  by  assuming  the  process 
is  steady  with  no  external  forces  and  that  the  sole  depletion  mechanism  for  the  neutral  velocity  classes  is 
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electron  collisional  ionization.  Implicit  in  these  assumptions  is  that  the  ions  and  neutrals  do  not  significantly 
interact. 

Equation  18  may  be  easily  integrated  with  respect  to  nf(u)  to  produce  an  analytic  solution  if  the  ioniza¬ 
tion  rate  coefficient  Sj  and  the  electron  number  density  ne  are  assumed  to  be  constant.  In  this  case,  the  pop¬ 
ulation  of  the  u  velocity  class  nf(u)  exponentially  decays  with  the  spatial  variable  z  moderated  by  the 
value  of  u . 

nf(u)  -  e-z/u  (2°) 

The  qualitative  conclusion  drawn  from  this  straightforward  analysis  of  the  one  dimensional  Boltz¬ 
mann  equation  is  that  the  neutral  density  in  a  model  plasma  is  depleted  along  the  Z  axis  due  to  neutral-elec¬ 
tron  collisional  ionization.  This  relative  depletion  depends  on  the  velocity  class.  Neutrals  in  the  lower 
velocity  classes  have  a  greater  probability  of  being  ionized  when  passing  through  the  zone  of  ionization  than 
neutrals  of  higher  velocity  classes.  Therefore,  the  apparent  acceleration  of  the  neutrals  in  Fig.  13  is  most 
likely  the  result  of  the  depletion  of  the  slower  moving  neutrals  rather  than  an  acceleration  process. 

Conclusions 

Measurements  of  xenon  ion  and  neutral  velocities  were  performed  in  the  plume  and  into  the  interior 
of  the  thruster  through  a  1  mm  wide  slot  in  the  outer  insulator  wall.  From  these  measurements,  information 
on  propellant  energy  deposition,  electric  field  strength,  and  flow  divergence  were  extracted.  Xenon  ion 
velocity  measurements  of  axial  velocity  both  inside  and  outside  the  thruster  as  well  as  radial  velocity  mea¬ 
surements  outside  the  thruster,  were  performed  using  LIF  with  nonresonant  signal  detection  using  the  xenon 
ion  5d[4]7/2-6p[3]5/2  electronic  transition  while  monitoring  signal  from  the  6s[2]3/2- 6p[3]5/2  transi¬ 
tion.  Neutral  velocity  measurements  were  similarly  performed  in  the  interior  of  the  Hall  thruster  using  the 
6s[3/2]°  -6p[3/2]2  transition  with  resonance  fluorescence  collection.  Most  velocity  measurements  used 
partially  saturated  fluorescence  to  improve  the  signal  to  noise  ratio.  One  radial  trace  of  the  xenon  ion  transi¬ 
tion  was  taken  in  the  linear  fluorescence  region  and  yielded  a  plume  ionic  translational  temperature  between 
400  and  800  K.  However,  since  the  hyperfine  structure  constants  are  not  known  for  the  5d[4]7/2  level,  the 
constants  for  the  5d[3]7/2  level  were  used  instead.  This  result  should  therefore  be  viewed  with  caution.  An 
upper  limit  on  the  kinetic  temperature  using  only  the  full  width  at  half  maximum  assuming  no  hyperfine 
splitting  yields  a  temperature  of  1,700  K.  From  the  neutral  velocity  measurements,  the  neutrals  appear  to  be 
accelerated  within  the  thruster.  Conclusions  drawn  from  analysis  of  the  one  dimensional  Boltzmann  equa¬ 
tion  imply  that  neutrals  are  depleted  along  the  Z  axis  due  to  neutral-electron  collisional  ionization,  and  the 
relative  depletion  depends  on  the  velocity  class  of  the  neutral  atom.  This  preferentially  depletes  the  lower 
velocity  classes  producing  the  apparent  acceleration  of  the  neutrals. 
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Fig.  6.  Saturation  curve  of  radial  xenon  ion  measurements  for  an  anode  discharge  voltage  of  200  V  at  Z  =  13 
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V  discharge  voltage. 
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Fig.  13.  Axial  neutral  velocity  measurements  at  D  =  0  mm  for  discharge  voltages  of  100  V,  160  V,  200  V,  and 
250V. 
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Anomalous  electron  mobility  in  a  coaxial  Hall  discharge  plasma 

Nathan  B.  Meezan,  William  A.  Hargus,  Jr.,  and  Mark  A.  Cappelli 
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(Received  12  June  2000;  published  24  January  2001) 

A  comprehensive  analysis  of  measurements  supporting  the  presence  of  anomalous  cross-field  electron  mo¬ 
bility  in  Hall  plasma  accelerators  is  presented.  Nonintrusive  laser-induced  fluorescence  measurements  of  neu¬ 
tral  xenon  and  ionized  xenon  velocities,  and  various  electrostatic  probe  diagnostic  measurements  are  used  to 
locally  determine  the  effective  electron  Hall  parameter  inside  the  accelerator  channel.  These  values  are  then 
compared  to  the  classical  (collision-driven)  Hall  parameters  expected  for  a  quiescent  magnetized  plasma.  The 
results  indicate  that  in  the  vicinity  of  the  anode,  where  there  are  fewer  plasma  instabilities,  the  electron- 
transport  mechanism  is  likely  elastic  collisions  with  the  background  neutral  xenon.  However,  we  find  that  in 
the  vicinity  of  the  discharge  channel  exit,  where  the  magnetic  field  is  the  strongest  and  where  there  are  intense 
fluctuations  in  the  plasma  properties,  the  inferred  Hall  parameter  departs  from  the  classical  value,  and  is  close 
to  the  Bohm  value  of  (w(.er)effa«16.  These  results  are  used  to  support  a  simple  model  for  the  Hall  parameter 
that  is  based  on  the  scalar  addition  of  the  electron  collision  frequencies  (elastic  collision  induced  plus  fluctua¬ 
tion  induced),  as  proposed  by  Boeuf  and  Garrigues  [J.  Appl.  Phys.  84,  3541  (1998)].  The  results  also  draw 
attention  to  the  possible  role  of  fluctuations  in  enhancing  electron  transport  in  regions  where  the  electrons  are 
highly  magnetized. 
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I.  INTRODUCTION 

The  use  of  coaxial  Hall  discharges  as  high  specific  im¬ 
pulse  plasma  accelerators  for  satellite  propulsion  date  back  to 
the  early  1960’s.  In  a  Hall  discharge  (Fig.  1),  the  plasma  is 
sustained  in  imposed  orthogonal  electric  and  magnetic  fields. 
The  discharge  electrons,  a  large  fraction  of  which  are  emitted 
by  an  external  cathode,  are  magnetized,  whereas  the  more 
massive  propellant  ions,  usually  xenon,  are  not.  Conse¬ 
quently,  the  electrostatic  fields  established  by  the  retarded 
electron  flow  accelerate  the  ions  to  high  velocities,  typically 
50%-60%  of  the  discharge  voltage.  In  a  coaxial  geometry, 
the  electrons  are  constrained  to  move  in  the  closed,  azi¬ 
muthal  EXB  drift,  with  cross-field  diffusion  providing  the 
necessary  current  to  sustain  the  discharge. 

The  first  detailed  investigation  of  the  properties  of  these 
coaxial  Hall  discharges  suggested  that  the  axial  electron  cur¬ 
rent  density  is  greater  than  that  expected  by  classical  trans¬ 
port  [1].  This  led  to  the  first  speculation  of  the  presence  of  an 
“anomalous”  cross-field  drift  of  electrons,  due  in  part  to 
fluctuations  in  the  electric  field  and  plasma  density.  Fluctua¬ 
tions  in  the  external  circuit  and  plasma  properties  in  these 
devices  are  now  well  documented  [2-5],  and  recent  numeri¬ 
cal  simulations  of  Hall  discharges  generally  require  the  ad 
hoc  introduction  of  an  anomalous  electron  mobility  to  pro¬ 
duce  reasonably  accurate  results  [6-9].  Despite  many  experi¬ 
mental  and  theoretical  studies,  the  precise  nature  of  the  in¬ 
stabilities  responsible  for  these  fluctuations  and  the 
mechanism  for  their  excitation  are  not  yet  well  understood. 
Furthermore,  other  mechanisms,  including  collisions  be¬ 
tween  energetic  electrons  and  the  channel  wall,  and  experi¬ 
mental  facility  effects  such  as  background  gas  ingestion  may 
also  play  important  roles  in  electron  transport. 

In  this  paper,  we  examine  the  electron-transport  mecha¬ 
nism  in  a  coaxial  Hall  discharge  through  a  collection  of  ap¬ 
plied  diagnostics  that  serve  to  map  the  electron  mobility  in 


the  discharge  channel.  The  two-dimensional  electron  mo¬ 
mentum  equation  for  a  weakly  ionized  plasma  gives  the  elec¬ 
tron  current  density. 


-enM^  =  enP 


\  ?n£  =  en  /E*\  1 

BrJo)ce  \BrJ  o)ceT 


(1) 


Here,  ne  is  the  electron  density,  uez  is  the  electron  drift  ve¬ 
locity,  vne  is  the  momentum-transfer  collision  frequency  for 
neutral-electron  collisions,  r  is  the  mean  time  between  colli¬ 
sions,  and  the  electron  cyclotron  frequency  a)ce»vne .  It  is 
assumed  that  the  applied  electric  field  Ez  and  magnetic  field 
strength  Br  are  in  the  axial  and  radial  directions,  respec- 


FIG.  1.  Schematic  of  a  typical  Hall  discharge  accelerator. 
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FIG.  2.  Cutaway  schematic  of  the  Hall  discharge  studied,  (a) 
Front  magnetic  plate;  (b)  central  magnetic  plate;  (c)  anode  gas  in¬ 
put;  (d)  alumina  insulator  wall;  (e)  magnet  core;  and  (f)  magnetic 
windings. 


tively.  From  Eq.  (1),  we  see  that  the  cross-field  electron  mo¬ 
bility  is  inversely  proportional  to  the  electron  Hall  parameter, 
o>r,r, 


Il-L 

Br)<*>ceT 


(2) 


Equation  (1)  can  be  rearranged  to  give 

1  _tBr 
(a>ceT)  \EZ 

The  Hall  parameter  (or  electron  mobility)  can  therefore  be 
determined  at  any  point  in  the  discharge  channel  where  the 
plasma  properties  on  the  right-hand  side  of  Eq.  (3)  are 
known.  Furthermore,  the  resulting  Hall  parameter  can  be 
compared  to  that  calculated  on  the  basis  of  known  properties 
that  independently  determine  the  electron  cyclotron  fre¬ 
quency  and  the  mean  time  between  electron  collisions.  In 
this  paper,  we  review  our  recent  experimental  measurements 
of  plasma  properties  within  the  discharge  channel,  which 
permit  this  direct  comparison.  The  results  support  the  pres¬ 
ence  of  an  anomalous  transport  mechanism  in  regions  of  the 
plasma  known  to  support  fluctuations,  and  provide  the  first 
direct  quantitative  measure  of  the  axial  variation  in  the  elec¬ 
tron  mobility. 


II.  EXPERIMENTS 
A.  Prototype  Hall  discharge 

The  Hall  plasma  source  studied  here  was  constructed  at 
Stanford  as  a  test  bed  for  studying  the  discharge  physics  and 
is  not  intended  to  serve  as  an  operational  prototype  plasma 
accelerator.  A  cutaway  schematic  of  the  plasma  source  is 
shown  in  Fig.  2.  The  source  consists  of  an  annular  alumina 
channel  90  mm  in  diameter,  11  mm  in  width,  and  80  mm  in 
length.  A  magnetic  circuit  consisting  of  four  outer  coils,  one 


FIG.  3.  Axial  profile  of  radial  magnetic-field  strength  for  the 
Hall  discharge  studied. 

inner  coil,  and  three  iron  plates  provides  a  magnetic  field 
(mostly  radial  in  direction)  peaked  near  the  exit  of  the  dis¬ 
charge  channel,  as  shown  in  Fig.  3.  A  hollow  stainless-steel 
ring  with  32  holes  of  0.5-mm  diameter  serves  both  as  the 
anode  and  the  propellant  (gas)  input  of  the  discharge. 

A  commercial  hollow  cathode  (Ion  Tech  HCN-252)  is 
used  to  neutralize  the  resulting  ion  beam  and  support  the 
necessary  electric  field.  The  cathode  is  mounted  approxi¬ 
mately  2  cm  downstream  of  the  plasma  source  exit.  The 
cathode  body  was  kept  at  the  vacuum  chamber  ground  po¬ 
tential.  To  allow  optical  and  probe  access  inside  the  dis¬ 
charge  channel,  a  1  mm  wide  slot  was  cut  along  the  outer 
insulator  wall,  parallel  to  the  discharge  axis.  The  perfor¬ 
mance  of  the  discharge  is  not  noticeably  altered  by  the  pres¬ 
ence  of  the  slot. 

Two  independent  mass  flow  controllers  control  the  gas 
flow  to  the  anode  and  cathode.  Direct  current  power  supplies 
are  used  to  independently  power  the  discharge,  cathode 
heater,  cathode  keeper,  and  magnetic  windings.  The  dis¬ 
charge  (anode)  power  supply  is  always  operated  in  the 
voltage-regulated  limit.  The  experimental  facility  consists  of 
a  nonmagnetic  stainless-steel  tank  approximately  1  m  in  di¬ 
ameter  and  1.5  m  in  length  pumped  by  two  50-cm  diameter 
diffusion  pumps  and  backed  by  a  425  L/s  mechanical  pump. 
The  facility  is  equipped  with  several  viewports  and 
feedthroughs  for  optical  and  probe  access  to  the  discharge. 
The  base  pressure  of  the  facility  with  no  gas  flow  is  approxi¬ 
mately  10“ 6  Torr  as  measured  by  an  ionization  gauge  uncor¬ 
rected  for  gas  species.  The  background  pressure  during  dis¬ 
charge  operation  at  a  xenon  flow  rate  of  2.3  mg/s  is  typically 
10“4  Torr.  Although  this  pressure  is  an  order  of  magnitude 
lower  than  that  of  Janes  and  Lowder  [1],  it  is  still  consider¬ 
ably  higher  than  chamber  pressures  that  are  generally  accept¬ 
able  for  the  collection  of  accelerator  performance  data 
[10,1 1].  The  ingestion  of  background  gas  near  the  exit  of  the 
discharge  channel  may  influence  the  discharge  characteris¬ 
tics. 


B.  Plasma  diagnostic  techniques 

A  variety  of  optical  and  probe  diagnostics  were  utilized  to 
measure  the  plasma  properties  needed  for  the  determination 
of  the  effective  Hall  parameter  [using  Eq.  (3)].  The  potential 
distribution  in  the  discharge  was  measured  using  an  emissive 
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probe.  This  technique  is  described  in  more  detail  in  another 
paper  [12],  and  is  only  briefly  summarized  here.  A  tungsten 
filament  inserted  into  the  discharge  channel  was  heated  using 
a  variable-current  alternating  current  (ac)  power  supply  and  a 
second,  floating  step-down  transformer.  The  filament  current 
was  increased  until  the  potential  of  the  floating  transformer 
with  respect  to  ground  reached  a  nearly  saturated  value.  This 
value  was  taken  as  the  plasma  potential  [13].  This  method  of 
finding  the  plasma  potential  is  less  accurate  than  the  more 
common  “inflection  point”  method  [14];  however,  this  al¬ 
ternative  method  involves  the  collection  of  electron  current 
by  the  probe.  Due  to  the  small  size  of  the  discharge  channel, 
it  has  been  our  experience  that  collecting  electrons  from  the 
plasma  can  severely  perturb  the  operation  of  the  discharge.  A 
probe  collecting  electron  current  causes  a  much  greater  dis¬ 
turbance  on  discharge  performance  than  a  floating  probe  or 
one  collecting  ion  current.  The  floating  potential  of  the 
plasma  was  also  recorded  at  each  point  using  the  same  probe 
without  current  passing  through  the  filament.  The  magnetic- 
field  distribution  inside  the  channel  was  measured  with  the 
discharge  off  using  a  Hall  effect  sensor  probe. 

The  electron  current  density  is  determined  from  the  total 
discharge  current  density  (discharge  current  divided  by  the 
channel  cross  sectional  area)  less  the  ion  current  density: 

J  ez  HA  channel  '  (^) 

The  ion  current  density  is  determined  either  directly  using  a 
guarded  ion  flux  probe,  or  from  a  measurement  of  the  elec¬ 
tron  (ion)  density  using  a  cylindrical  Langmuir  probe  (dis¬ 
cussed  below)  and  from  laser-induced  fluorescence  (LIF) 
measurements  of  the  xenon  ion  velocity.  These  LIF  measure¬ 
ments  are  described  in  detail  elsewhere  [12,15].  They  in¬ 
volve  propagating  a  narrow-linewidth,  continuous-wave  la¬ 
ser  beam  through  the  plasma  along  the  discharge  axis  and 
measuring  the  fluorescence  signal  as  a  function  of  laser  fre¬ 
quency.  Due  to  the  Doppler  effect,  the  absorption  line  center 
is  shifted  from  that  of  a  stationary  plasma  reference.  The 
bulk  velocity  of  the  absorbing  xenon  ions  can  then  be  deter¬ 
mined  from  this  frequency  shift.  In  the  studies  reported  here, 
a  Coherent  899-21  tunable  titanium  sapphire  laser  was  used 
to  excite  the  5^/[4]7/2— 6/7[3]5/2  transition  in  Xell  at  834.7 
nm,  with  nonresonance  fluorescence  detected  at  541.9  nm.  A 
similar  method  was  also  used  to  measure  the  velocity  of 
excited  neutral  xenon  atoms  using  the  65*[ |-]2— 6/?[ f  ]2  tran¬ 
sition  at  823.2  nm,  only  in  this  case,  resonance  fluorescence 
is  detected.  The  neutral  atom  velocity  together  with  the  over¬ 
all  xenon  mass  flow  rate  served  to  provide  a  measure  of  the 
neutral  xenon  number  density,  from  which  the  classical  Hall 
parameter  can  be  determined. 

The  electron  density  in  the  discharge  was  measured  using 
two  different  techniques.  In  the  first  set  of  experiments,  a 
cylindrical  Langmuir  probe  was  inserted  into  the  discharge 
channel  through  the  optical  access  slot.  The  probe  was  a 
0.5-mm  diameter  tungsten  wire  bent  at  a  right  angle  so  that  a 
5-mm  length  of  the  probe  was  parallel  to  the  ion  flow  in  the 
channel.  The  ion-saturation  current  was  collected  at  a  probe 
bias  of  —46  V  and  the  electron  (ion)  density  was  determined 
using  a  thin  (“free-fall”)  sheath  approximation  [16]  for  ion 
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collection.  This  required  an  estimate  of  the  electron  tempera¬ 
ture,  as  described  below.  Inserting  the  probe  through  the  ac¬ 
cess  slot  minimized  its  disturbance  of  the  plasma,  but  also 
prevented  measurements  near  the  exit  of  the  discharge. 

Shortcomings  of  this  probe  characterization  method  in¬ 
cluded  spatial  averaging  along  the  length  of  the  probe, 
“shadowing”  of  the  backside  of  the  probe  due  to  wake  for¬ 
mation,  and  the  use  of  a  relatively  simple  ion  current  collec¬ 
tion  theory  for  a  high-velocity  plasma  in  the  presence  of  a 
magnetic  field.  The  effects  of  the  magnetic  field  on  particle 
collection  by  the  probe  can  be  measured  through  the  electron 
and  ionic  Larmor  radii.  For  the  Hall  discharge,  re***a  and 
r(>w>ay  where  re  and  rt  are  the  electronic  and  ionic  Lar¬ 
mor  radii,  a  is  the  probe  radius,  and  w  is  the  channel  width. 
While  electron  collection  will  be  affected  by  the  magnetic 
field,  the  ions  are  essentially  not  magnetized.  Therefore,  as 
pointed  out  by  Batani  et  al  [17],  the  effect  of  the  magnetic 
field  on  a  measurement  of  the  electron  density  from  the  ion 
current  will  be  negligible  compared  to  the  20%  accuracy 
expected  by  this  probe  theory.  Ideally,  the  entire  Langmuir 
probe  trace  would  be  used  for  determining  the  electron  tem¬ 
perature,  plasma  density,  and  plasma  potential.  As  men¬ 
tioned  above,  scanning  the  probe  voltage  into  the  electron- 
collection  region  is  not  desirable  in  this  discharge,  as  it  may 
strongly  perturb  the  plasma,  and  requires  an  appropriate  the¬ 
oretical  treatment  of  probe  traces  when  electrons  are  strongly 
magnetized. 

As  mentioned  above,  the  ion  current  density  is  also  deter¬ 
mined  directly  using  a  guarded  ion-flux  probe  [18].  The  ion 
(hence  electron)  density  can  be  determined  from  these  probe 
measurements  and  the  measurements  of  the  ion  velocity.  The 
ion  probe  was  constructed  of  1.5-mm  diameter  tungsten  rod 
surrounded  by  an  alumina  sleeve.  A  tantalum  guard  ring  bi¬ 
ased  to  the  same  potential  as  the  probe  was  mounted  around 
the  center  probe  to  establish  a  one-dimensional  sheath  across 
the  probe  collection  area,  and  to  minimize  outer  sheath  col¬ 
lection  area  growth  with  probe  bias.  At  all  locations  within 
the  discharge,  the  probe  was  biased  at  20  V  below  the  local 
floating  potential.  It  is  noted,  however,  that  inserting  the 
probe  into  the  discharge  channel  caused  a  significant  increase 
in  the  total  discharge  current,  up  to  20%.  The  intrusive  per¬ 
turbation  of  this  probe  precluded  measurements  upstream  of 
the  discharge  exit  plane. 

Measurements  of  plasma  fluctuations  were  also  made  to 
give  insight  into  the  role  of  oscillations  in  electron  transport. 
Low-impedance  shielded  Langmuir  probes  were  used  to 
measure  fluctuations  in  ion-saturation  current,  which  roughly 
correspond  to  fluctuations  in  plasma  density.  These  probes 
were  designed  to  detect  propagating  plasma  waves  in  the 
Hall  discharge,  experiments  described  in  great  detail  in  an¬ 
other  paper  [19].  For  this  experiment,  one  probe  was  con¬ 
nected  to  a  direct  current  (dc)  coupled  digital  oscilloscope 
through  a  50  D  terminator.  Due  to  the  high  floating  potential 
inside  the  discharge,  this  grounded  probe  collected  current  in 
the  ion-saturation  regime  of  the  Langmuir  probe  trace. 

III.  RESULTS  AND  ANALYSIS 

The  plasma  potential  profiles  as  measured  by  the  hot- 
filament  probe  are  shown  in  Fig.  4  for  the  three  operating 
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FIG.  4.  Plasma  potential  profiles  measured  by  the  emissive 
probe. 

conditions  investigated  in  this  paper.  The  floating  potential 
of  an  emitting  probe  should  be  slightly  below  the  actual 
plasma  potential,  since  the  emission  will  never  be  enough  to 
neutralize  the  sheath  [14].  It  can  be  seen  from  Fig.  4  that  as 
the  probe  approaches  the  anode,  the  measured  plasma  poten¬ 
tial  comes  quite  close  to  the  applied  anode  potential.  With 
the  probe  5  mm  away  from  the  anode,  the  measured  values 
are  within  4  V  of  the  applied  potential,  suggesting  that  the 
uncertainty  in  the  measured  potential  is  within  5  V,  and  pos¬ 
sibly  less,  depending  on  the  strength  of  the  anode  fall.  The 
electric  field  inside  the  discharge  channel,  shown  in  Fig.  5, 
was  calculated  from  the  plasma  potential  data  using  the  re¬ 
lation  Ez-  - d(f>ldz .  Smooth  curves  were  fit  to  the  potential 
profiles  before  differentiation  to  prevent  noise  amplification. 
It  is  noteworthy  that  the  peak  in  the  electric  field  coincides 
closely  with  the  peak  in  the  magnetic  field,  although  the 
distribution  is  more  narrow,  suggesting  that  the  electron  E 
XB  drift  velocity,  ud^Ez!Br,  is  also  highly  peaked  near  z 
~  —  10  mm. 

The  plasma  potential  data  was  used  in  conjunction  with 
floating  potential  data  (not  shown)  to  calculate  the  electron 
temperature: 

kTe  l  nii  \ 

<» 

Here,  we  assume  that  the  ions  enter  the  sheath  at  the  Bohm 
velocity. 
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FIG.  5.  Electric-field  profiles  calculated  from  plasma  potential. 
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FIG.  6.  Electron  temperature  profiles  calculated  from  plasma 
and  floating  potential  data. 

The  electron  temperature  is  shown  in  Fig.  6.  These  tem¬ 
perature  profiles  are  used  in  the  remaining  calculations 
needed  to  estimate  the  Hall  parameter.  The  electron  tempera¬ 
ture  is  seen  to  peak  above  16  eV.  Although  the  temperature 
could  be  in  error  by  a  few  volts  due  to  underestimation  of  the 
plasma  potential,  the  measured  values  do  agree  well  with 
those  from  other  investigators  in  similar  discharges  [20,21]. 
The  peak  in  the  electron  temperature  is  coincident  with  the 
peaks  in  the  electric  and  magnetic  fields,  as  expected  since 
the  Ohmic  dissipation  will  scale  as  JezEz .  We  see  an  unex¬ 
pected  rise  in  the  electron  temperature  near  the  anode  for  the 
higher  power  cases  studied,  most  likely  due  to  the  formation 
of  an  anode  glow  region,  as  the  production  of  electrons  is 
needed  in  this  region  to  support  the  higher  current  densities. 
These  striking  distributions  suggest  that  the  discharge  length 
is  not  optimized,  and  can  be  shortened  to  40  mm  for  more 
efficient  operation  as  a  plasma  accelerator  at  the  higher- 
voltage  conditions. 

The  electron  temperature  is  used  to  determine  the  electron 
(ion)  density  from  the  ion  saturation  current  collected  with 
the  Langmuir  probe  discussed  above, 

fkT~ 

Jis  =  0.61ene\  — (6) 
v  /n, 

The  electron-density  profiles  measured  from  the  two  probe 
methods  are  plotted  together  in  Fig.  7.  As  mentioned  above, 
the  error  in  the  electron  temperature  will  have  a  small  impact 
on  this  measurement.  For  example,  an  error  as  large  as  2  eV 
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FIG.  7.  Electron-density  profiles  from  probe  experiments. 
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FIG.  8.  Ion  velocity  profiles  measured  by  laser-induced  fluores¬ 
cence  velocimetry. 

in  a  measured  electron  temperature  of  6  eV  will  lead  to  an 
error  of  15% -20%  in  the  electron  density — within  the  ex¬ 
pected  accuracy  of  the  method. 

The  intrusive  nature  of  the  guarded  ion  probe  precluded 
measurements  with  this  probe  within  the  discharge  channel. 
Furthermore,  the  Langmuir  probe  inserted  through  the  slot  in 
the  insulator  was  unable  to  collect  data  within  a  distance  of 
10  mm  upstream  of  the  channel  exit,  because  of  the  presence 
of  the  magnetic  pole  piece.  As  a  result,  the  two  probe  mea¬ 
surements  were  indeed  complementary,  and  a  cubic-spline 
interpolation  was  used  to  estimate  the  plasma  density  in  the 
range  between  the  discharge  channel  exit  and  the  location  of 
the  base  of  the  front  magnetic  plate.  The  resultant  spline  fits 
are  also  shown  in  Fig.  7.  The  electron  number  densities  are 
seen  to  peak  in  the  same  vicinity  where  the  other  plasma 
properties  reach  their  maximum  values.  The  error  in  the 
electron-density  measurement  inside  the  channel  is  quite  sig¬ 
nificant,  but  it  is  also  bound.  The  values  have  an  upper  bound 
of  100%  propellant  utilization  (i.e.,  the  ion  flow  leaving  the 
channel  equals  the  mass  flow  entering  the  channel),  and  the 
ion-flux  measurements  effectively  provide  a  lower  bound  on 
the  density  near  the  exit  plane. 

The  ion  velocity  profiles  measured  using  LIF  are  shown 
in  Fig.  8.  Just  downstream  of  the  high-temperature  ionization 
region,  the  produced  ions  are  seen  to  accelerate  to  high  ve¬ 
locities  consistent  with  the  measured  plasma  potential  distri¬ 
butions. 

The  data  given  above  for  Br ,  Ez ,  ne ,  and  Vt  (ion  veloc¬ 
ity)  is  enough  to  calculate  the  effective  Hall  parameter  along 
the  Hall  discharge  channel  from  Eqs.  (3)  and  (4).  This  Hall 
parameter  can  be  compared  to  the  “classical”  Hall  param¬ 
eter  expected  based  on  the  electron  momentum-transfer  col¬ 
lision  frequency  vne .  For  this  calculation,  the  neutral  gas 
density  (the  dominant  electron  collision  partner)  is  deter¬ 
mined  from  the  neutral  xenon  velocity  data,  shown  in  Fig.  9, 
and  the  discharge  mass  flow  rate, 


nijA 


channel 


~neVi\y~' 


(7) 


The  axial  variation  in  the  inferred  neutral  density  is  shown  in 
Fig.  10. 
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FIG.  9.  Neutral  xenon  velocity  profiles  measured  by  laser- 
induced  fluorescence  velocimetry. 

The  dip  in  the  neutral  density  near  the  exit  plane  at  higher 
discharge  voltages  is  due  to  the  measurement  of  a  significant 
increase  in  the  ion  flux,  reflecting  the  ionization  process  tak¬ 
ing  place  between  z=~5  and  as  -10  mm.  It  is  noteworthy 
that  the  calculation  of  the  neutral  gas  density  does  not  take 
ingestion  of  background  gas  from  the  vacuum  chamber  into 
account.  Background  chamber  neutral  xenon  ingestion  would 
cause  an  underestimate  of  the  neutral  xenon  density,  as  the 
actual  mass  flow  rate  would  be  higher.  These  measurements 
of  xenon  density  should  be  interpreted  as  lower  limits  on  the 
actual  values. 

The  electron  momentum-transfer  collision  frequency  was 
calculated  using  the  cross  section  given  in  the  SIGLO©  da¬ 
tabase  [22]  assuming  a  Maxwellian  electron  energy  distribu¬ 
tion,  using  the  electron  temperature  measurements  described 
above. 


IV.  DISCUSSION 


The  deduced  Hall  parameter  is  compared  to  the  “classi¬ 
cal”  Hall  parameter  (that  determined  by  the  direct  calcula¬ 
tion  of  the  electron  momentum-transfer  collision  frequency), 
in  Fig.  11.  In  the  figure,  we  plot  the  inverse  of  the  Hall 
parameter  (1  lo)cer),  since  this  is  proportional  to  the  electron 
cross-field  mobility  [see  Eq.  (2)]. 

It  is  apparent  from  Fig.  11  that  elastic  electron-neutral 
collisions  are  the  likely  mechanisms  for  cross-field  electron 
transport  at  locations  between  the  anode  (z=  -80mm)  and 


~  io‘'  1 


10  i 


(0 
c 
a> 

O 
c 
o 
c 
to 

rrt  10  1 


3 

to 

Z 


10 


|ao 

□ 

°  0  D  □ 

|A  A  2  °  o  °  □ 

A  A  O  n  □ 


O  100V 
□  160V 
A  200V 


28 


aa  §0q 

A  A  Q 


A  0 


□  A 
A  □ 


-60 


-40 


-20 


Distance  from  Exit  (mm) 


026410-5 


FIG.  10.  Neutral  xenon  density  profiles. 
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FIG.  11.  Effective  and  classical  inverse  Hall  parameter  profiles 
inside  Hall  discharge  channel. 


FIG.  13.  Comparison  of  statistical  fluctuation  analysis  to  experi¬ 
mental  inverse  Hall  parameter  profile  at  100  V  operation. 


approximately  40  mm  upstream  of  the  discharge  channel 
exit.  It  is  noteworthy  that  this  location  is  where  the  magnetic- 
field  strength  reaches  approximately  l/5th  of  its  peak  value 
(see  Fig.  3).  We  also  see  that  beyond  this  location,  there  is  a 
striking  departure  from  this  classical  value,  although  the  un¬ 
certainty  in  the  absolute  measurements  is  still  significant,  a 
result  mainly  attributable  to  the  uncertainty  in  the  total  dis¬ 
charge  current  caused  by  the  intrusive  nature  of  the  diagnos¬ 
tic  probes.  At  locations  downstream  of  the  z=  ~40mm  po¬ 
sition,  where  the  magnetic  field  reaches  a  maximum,  the 
inverse  Hall  parameter  is  in  remarkable  agreement  with  the 
Bohm  value  of  (ocer=  16  [23].  In  this  region,  we  [24-28], 
and  others  [2-5,29-32]  have  reported  on  the  existence  of 
intense  fluctuations  in  the  plasma  properties.  These  results 
are  supportive  of  the  conjecture  that  plasma  fluctuations  are 
partially  responsible  for  the  anomalous  cross-field  transport, 
although  additional  experiments  may  be  necessary  to  obtain 
more  accurate  quantitative  measurements  of  the  resulting 
mobility.  The  results  reported  on  here  do,  however,  provide 
the  first  direct  support  of  a  model  based  on  the  use  of  a 
“mixed”  mobility  (or  collision  frequency). 
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as  has  been  originally  proposed  by  Boeuf  and  Garrigues  for 
numerical  simulations  of  Hall  discharges  [7].  The  Bohm 


FIG.  12.  Comparison  of  mixed  mobility  model  to  experimental 
and  classical  inverse  Hall  parameter  profiles  at  200  V  operation. 


value  appears  to  be  reasonable  for  the  oscillation-induced 
part  of  the  mobility  (Fig.  12),  despite  the  lack  of  a  full  un¬ 
derstanding  of  the  instabilities  in  this  region  of  the  plasma. 

Additional  support  for  the  conjecture  that  oscillations  are 
at  least  partly  responsible  for  the  enhanced  cross-field  mobil¬ 
ity  can  be  found  in  statistical  analysis  of  probe  fluctuation 
data.  For  a  homogeneous  plasma  subjected  to  homogeneous 
turbulent  fluctuations  in  the  plasma  density,  the  cross-field 
electron  diffusion  coefficient  (which  is  directly  proportional 
to  the  inverse  of  the  effective  Hall  parameter)  can  be  ex¬ 
pressed  as  [33] 


dloc7 - \ 

(WT) 


eff 


((ne~(ne))2) 


(ne)2 


(9) 


The  quantity  inside  the  square  brackets  is  the  mean-square 
deviation  of  the  observed  plasma  density  fluctuation.  As 
mentioned  above,  the  raw  fluctuation  data  collected  by  an 
ion-current  probe  is  proportional  to  the  local  plasma  density. 
Equation  (9)  was  therefore  applied  directly  to  the  raw  probe 
data  inside  the  discharge  channel.  This  “statistical  Hall  pa¬ 
rameter”  is  compared  to  the  experimental  value  in  Fig.  13. 
The  agreement  is  remarkable,  although  the  uncertainty  in  the 
measured  value  of  the  Hall  parameter  is  large.  More  signifi¬ 
cantly,  the  trend  in  the  mean-square  deviation  of  the  plasma 
density  follows  the  experimental  value  very  closely,  suggest¬ 
ing  that  the  role  of  fluctuations  may  be  important. 

The  relatively  large  error  bars  in  Fig.  11  are  a  conse¬ 
quence  of  the  number  of  different  data  used  to  calculate  the 
Hall  parameter.  The  largest  source  of  error  is  the  measure¬ 
ment  of  the  total  current  through  the  discharge.  At  a  given 
voltage,  the  current  through  the  discharge  used  in  this  study 
can  drift  up  to  10%  in  a  half  hour,  presumably  due  to  fluc¬ 
tuations  in  the  facility  background  pressure.  Furthermore,  the 
insertion  of  the  ion  and/or  Langmuir  probe  gives  rise  to  a 
perturbation  that  can  be  as  high  as  20%  on  the  total  discharge 
current.  An  alternate,  nonintrusive  measurement  of  the  elec¬ 
tron  density  is  favored,  and  is  the  subject  of  our  present 
research.  Efforts  will  also  be  made  to  reduce  the  perturbation 
caused  by  the  emissive  potential  probe. 
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V.  CONCLUSIONS 

In  this  paper,  we  have  used  a  variety  of  plasma  diagnostic 
techniques,  both  intrusive  and  nonintrusive,  to  evaluate  the 
effective  electron  mobility  in  a  Hall  discharge.  We  provide 
strong  evidence  of  the  presence  of  an  anomalous  transport 
mechanism,  in  the  region  of  the  plasma  where  there  is  a 
strong  radial  magnetic  field,  and  where  we  have  seen  intense 
plasma  fluctuations.  The  Hall  parameter  has  been  shown  to 
approach  a  constant  value,  within  the  margin  of  error,  near 
the  Bohm  value  of  16  in  this  region.  Although  these  results 
fall  short  of  identifying  fluctuations  in  plasma  density  and/or 
temperature  as  the  source  of  this  enhanced  cross-field  elec¬ 
tron  mobility,  they  do  suggest  that  fluctuations  in  plasma 
properties  can  account  for  the  observed  axial  electron  current 
flow.  The  Hall  discharge  accelerator  falls  within  a  difficult 
regime  for  plasma  diagnostics — the  device  is  small,  low- 
density,  fluctuating,  and  magnetized,  making  it  difficult  to 
interpret  probe-based  measurements  of  electron  density  us¬ 
ing  available  probe  theory.  Furthermore,  although  the 
electron-density  measurements  themselves  are  not  the  largest 
source  of  uncertainty,  we  find  that  the  probes  used  for  these 
measurements  have  a  weak  perturbing  effect  on  the  total  dis¬ 
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charge  current,  which  does  strongly  contribute  to  the  mea¬ 
surement  accuracy.  Future  studies  will  focus  on  the  develop¬ 
ment  of  nonintrusive  measurements  of  plasma  density  to 
improve  the  measurement  accuracy. 

The  most  significant  finding  is  that  the  results  are  in  re¬ 
markable  agreement  with  a  simple  model  for  the  Hall  param¬ 
eter  that  is  based  on  the  scalar  addition  of  the  electron  colli¬ 
sion  frequencies  (elastic  collision  induced  plus  fluctuation 
induced),  as  first  proposed  by  Boeuf  and  Garrigues  [7].  The 
transition  from  a  fluctuation-dominated  mobility  (Bohm-like) 
to  collision  induced  mobility  is  found  to  occur  at  the  point 
where  the  magnetic  field  is  approximately  20%  of  its  peak 
value,  and  where  our  other  studies  [19,28]  identify  the  spatial 
onset  of  strong,  quasicoherent  azimuthally  propagating 
waves  in  plasma  density  and  potential. 
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Abstract 

Experimental  results  are  presented  for  studies  of  low  (2-20  kHz)  and  intermediate-frequency  (20-100 
kHz)  oscillations  in  crossed-field  closed-electron  drift  Hall  discharges.  Conditional  sampling  using 
two  electrostatic  probes  is  used  to  identify  and  extract  properties  of  coherent  structures  associated 
with  the  propagation  of  azimuthal  and  longitudinal  instabilities  within  the  discharge  channel.  The 
azimuthal  component  phase  velocities  are  determined  for  a  wide  range  of  wave  frequencies  and  over 
characteristic  regimes  of  operation  of  these  devices.  A  variety  of  propagation  modes  are  observed 
and  analyzed,  including  the  appearance  of  an  induced  mode  due  to  the  presence  of  the  probes 
themselves.  This  later  result  is  believed  to  be  the  first  direct  evidence  of  how  fluctuations  can  be 
influenced  in  these  Hall  discharges  using  relatively  simple  actuation  methods. 

Index  Terms  -  Hall  discharges,  Langmuir  probes,  plasma  propulsion,  plasma  instabilities. 
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I.  INTRODUCTION 

Closed  electron-drift  Hall  discharges  are  low-pressure  (-0.01  -  0.1  Pa),  weakly  collisional, 
magnetized  plasma  sources  that  generate  a  relatively  high  velocity  ion  beam  suitable  for  use  in  space 
propulsion  applications.  A  particular  class  of  Hall  plasma  thrusters,  the  so-called  “Stationary  Plasma 
Thruster,”  has  been  used  in  a  number  of  space  missions  in  the  former  Soviet  Union  [1].  A  variation 
of  this  Hall  thruster  that  has  an  annular  (co-axial)  discharge  channel  100  mm  in  diameter  -  the  SPT- 
100,  has  a  high  specific  impulse  (1100  -  2000  sec),  operates  at  moderately  high  thrust  levels  (55  - 
150  mN),  and  has  an  exceptionally  high  thrust  efficiency  (40-60%).  Because  of  this  performance, 
this  plasma  source  is  now  being  aggressively  developed  for  use  in  station  keeping  applications  on 
western  satellites. 

In  a  typical  co-axial  geometry  Hall  discharge  (Fig.  1),  the  plasma  is  sustained  in  imposed 
orthogonal  electric  and  magnetic  fields.  The  discharge  electrons,  a  large  fraction  of  which  are 
emitted  by  an  external  cathode,  are  magnetized,  whereas  the  more  massive  propellant  ions,  usually 
xenon,  are  not.  Consequently,  the  electrostatic  fields  established  by  the  retarded  electron  flow 
accelerate  the  ions  to  high  velocities,  typically  50-60%  of  the  discharge  voltage  (-1 00-300 V).  The 
maximum  acceleration  occurs  in  the  region  between  the  magnetic  poles,  where  the  magnetic  field  is 
a  maximum.  In  a  co-axial  geometry,  the  electrons  are  constrained  to  move  in  the  closed,  azimuthal  E 
x  B  drift,  with  cross-field  diffusion  providing  the  necessary  current  to  sustain  the  discharge.  An 
annular  ceramic  channel  confines  the  electron  flow  towards  the  anode,  located  at  its  base. 

It  is  widely  known  that  the  Hall  discharge  plasmas  exhibit  a  rich  spectrum  of  fluctuations  in 
plasma  properties  [2].  While  it  is  not  yet  known  if  and  how  these  fluctuations  can  impact  the 
performance  of  a  Hall  thruster,  it  is  believed  that  fluctuations  in  the  bulk  plasma  properties  are  partly 
responsible  for  anomalous  electron  transport  across  the  imposed  magnetic  field  [3],  While  some 
studies  characterizing  the  presence  and  origin  of  these  fluctuations  and  their  possible  control  were 
published  in  the  mid  60’s  to  mid  70’s  [3-9],  they  have  received  increased  attention  recently  [10-15], 
as  there  is  a  growing  need  to  extend  and  enhance  the  performance  of  these  thrusters  for  a  broader 
range  of  space  missions. 

Our  research  is  motivated  by  the  possibility  of  enhancing  Hall  thruster  operation  by  the  active 
control  of  these  fluctuations,  or  by  the  passive  suppression  of  the  fluctuations  in  regions  of  the 
discharge  channel  where  a  reduction  in  electron  current  is  desired.  To  do  so,  it  is  important  that  the 
nature  of  these  fluctuations  is  adequately  understood.  In  this  paper,  we  present  an  experimental 
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characterization  of  the  propagation  of  these  disturbances  over  a  range  of  Hall  discharge  operation, 
and  demonstrate  that  it  is  possible  to  affect  the  oscillations  within  the  discharge  by  the  insertion  of 
biased  probes. 

II.  EXPERIMENT 

The  Hall  discharge  plasma  source  used  in  this  study  is  a  laboratory  version  of  a  low-power  Hall 
thruster,  and  is  described  in  more  detail  in  previous  papers  [16-18].  This  particular  source  is 
intended  to  be  used  as  a  test  bed  for  studying  the  discharge  physics  and  not  to  serve  as  an  operational 
prototype  plasma  accelerator,  although  the  principal  design  is  similar  to  that  used  in  practise.  The 
time-averaged  plasma  properties  within  the  annular  discharge  channel  have  been  extensively 
characterized  for  a  range  of  operating  conditions  [16].  The  source  consists  of  an  annular  alumina 
channel,  90  mm  in  diameter,  1 1  mm  in  width,  and  80  mm  in  length.  A  magnetic  circuit  consisting  of 
four  outer  coils,  one  inner  coil,  and  three  iron  plates  provides  a  magnetic  field  (mostly  radial  in 
direction)  peaked  near  the  exit  of  the  discharge  channel,  as  shown  in  Fig.  2.  The  mapping  of  the 
axial  and  radial  components  of  the  magnetic  field  strength  within  the  channel  is  performed  with  a 
Hall  probe  in  the  absence  of  the  plasma.  The  radial  component  of  the  magnetic  field  is  found  to  drop 
off:  by  approximately  15%  at  the  inner  and  outer  walls  of  the  acceleration  channel.  Details  of  the 
two-dimensional  magnetic  field  distribution  for  magnet  currents  used  here  can  be  found  in  [18].  A 
hollow  stainless  steel  ring  with  32  holes  of  0.5  mm  diameter  serves  both  as  the  anode  and  the 
propellant  (gas)  input  of  the  discharge.  A  commercial  hollow  cathode  (Ion  Tech  HCN-252)  is  used 
to  neutralize  the  resulting  ion  beam  and  provide  the  necessary  electron  current  to  sustain  the 
discharge.  The  cathode  is  mounted  such  that  its  exit  aperture  is  approximately  2  cm  downstream  of 
the  plasma  source  exit.  The  cathode  body  was  kept  at  the  vacuum  chamber  ground  potential. 
During  discharge  operation,  the  1  m  diameter  by  1.5  m  long  chamber  within  which  the  discharge 
was  operated  was  maintained  at  a  pressure  of  10"4  Torr  (uncorrected  for  xenon),  as  indicated  by  an 
ionization  gauge.  Flanged  elbows  (0.5  m  in  diameter)  at  each  end  of  the  vacuum  chamber  supported 
two  0.5  m  diameter  diffusion  pumps,  operated  without  baffles,  for  maximum  pumping  speed 
(approximately  6000  1/s). 

For  the  measurements  reported  on  here,  the  xenon  flow  rate  was  2  mg/s  and  a  constant  solenoid 
winding  current  of  125  mA  provided  a  peak  magnetic  induction  of  about  100  G  at  5mm  upstream  of 
the  exit  plane  (inside  the  channel).  The  characteristics  of  the  plasma  fluctuations  on  centerline  within 
the  channel  were  studied  at  various  discharge  potentials. 
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Plasma  density  fluctuations  were  detected  by  two  azimuthally-placed  low-impedance  Langmuir 
probes  biased  negatatively  with  respect  to  the  plasma  potential  to  collect  the  ion  saturation  current, 
with  the  current  density  given  as  [19]: 


4=0-61 


(1) 


Here,  Te  is  the  electron  temperature,  ne  is  the  electron  number  density,  m,-  is  the  xenon  ion  mass,  e  is 
the  electron  charge,  and  k  is  the  Boltzmann  constant.  In  Eqn.  1,  it  is  assumed  that  the  distant  plasma 
is  quasi-neutral  (ne  =  n,),  and  that  the  ions  enter  the  collisionless  sheath  at  the  Bohm  velocity.  While 
the  fluctuations  in  collected  current  can  be  a  result  of  both  fluctuations  in  electron  density  and 
temperature,  in  prior  studies,  it  has  been  shown  that  these  low  frequency  disturbances  at  low 
operating  potential  are  largely  isothermal  [3,  7]. 

A  schematic  of  the  Langmuir  probe  construction  is  shown  in  Fig.  3.  The  exposed  part  of  the 
probe  consisted  of  a  0.254  mm  diameter,  3  mm  long  tungsten  wire.  The  base  of  the  probe  was  of  a 
complex  design,  intended  to  minimize  stray  capacitance  for  extended  frequency  response.  The  probe 
base  consisted  of  an  alumina  tube  directly  surrounding  the  tungsten  wire  (0.  508  mm  ED,  1.27  mm 
OD),  followed  by  a  stainless  steel  tube  (1.48  mm  ID,  1.58  mm  OD),  connected  to  the  braided  ground 
of  a  50  Q  co-axial  cable.  The  inner  tungsten  wire  was  connected  to  the  center  pin  of  the  co-axial 
cable.  This  outer  stainless  steel  tubing  served  as  a  shield  to  isolate  the  extended  probe  base  from 
disturbances  other  than  at  the  tip  of  the  probe.  The  stainless  steel  tubing  was  then  surrounded  by 
another  alumina  tube  (1.6  mm  ED,  2.3  mm  OD)  and  the  entire  inside  air  spaces  were  potted  with  an 
alumina  paste. 

The  probe  size  was  minimized  so  as  to  reduce  its  overall  capacitance,  and  for  the  studies 
reported  on  here,  was  approximately  10  cm  in  length.  Coaxial  transmission  line  feed-throughs 
provided  the  transfer  of  the  probe  signal  through  the  vacuum  chamber  to  a  National  Instruments 
DAQScope  5102  digital  oscilloscope  card  operating  at  an  800  kHz  sample  rate.  Data  was  stored  on  a 
personal  computer  for  future  processing.  The  probes  were  terminated  with  50  Q,  at  the  input  of  the 
oscilloscope  card.  In  this  way,  the  probe  tips  were  negatively  biased  (close  to  system  ground)  with 
respect  to  the  plasma,  which  is  predominantly  at  high  positive  potential,  since  the  anode  is 
maintained  at  a  minimum  of  86V,  and  since  the  near  exit  plane  potential  is  always  greater  than  40V 
due  to  the  cathode  fall  and  finite  electron  temperature.  Occasional  scans  of  the  Langmuir  probe 
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across  the  electron  retarding  regime  to  regions  of  high  negative  bias  verified  that  system  ground 
potential  was  always  in  the  ion  current  saturation  regime. 

The  probes  were  oriented  such  that  the  exposed  center  wire  axis  was  parallel  to  the  axial 
coordinate  of  the  thruster,  and  placed  midway  between  the  inner  and  outer  insulating  walls,  as 
illustrated  in  Fig.4.  Probe  translation  along  the  axial  direction  (the  direction  of  the  annular  channel) 
was  provided  by  a  translation  stage  driven  by  a  Slo-Syn  stepper  motor,  powered  and  controlled  by  a 
Compumotor  model  SX  Microstepping  Drive/Indexer  System.  The  azimuthal  positions  of  the  two 
probes  were  manually  set  before  the  beginning  of  each  experiment,  maintaining  an  equal  angular 
distance  from  the  two  nearest  outer  magnetic  coils.  For  this  study,  probe  separations  between  30  and 
180  degrees  were  examined.  In  most  cases,  the  axial  probe  locations  varied  between  a  distance  of 
approximately  35  mm  from  the  anode  (x  =  -45  mm)  and  10  mm  beyond  the  exit  plane  (x  =  10  mm), 
with  the  x  =  0  reference  taken  to  be  the  exit  plane  of  the  discharge,  and  negative  positions  implying 
that  the  probe  locations  are  inside  the  discharge  channel. 

A  fast  Fourier  transform  of  the  signal  from  each  probe  was  performed  across  a  temporal  window 
that  was  1.25  msec  wide  and  spaced  every  0.625  msec  in  time.  This  operation  was  used  to  convert 
the  digital  signal  in  the  time  domain  into  the  frequency  domain,  for  every  temporal  interval,  over  the 
length  of  the  temporal  data  set  (in  some  cases,  0.625  sec  long).  To  study  the  fluctuating  component 
to  the  signals,  the  mean  signal  was  removed  (numerically)  from  the  transient  signal  prior  to  its 
processing.  A  comparison  of  the  mean  signals,  and  the  power  spectrum  from  each  of  the  two  probes 
verified  that  the  probes  had  equal  sensitivities  (to  within  -10%)  and  frequency  responses.  To  avoid 
aliasing,  only  frequencies  below  400  kHz  (the  Nyquist  frequency,  defined  by  the  800  kHz  sampling 
rate  used  for  all  of  the  studies  described  here)  were  considered  for  interpretation.  In  the  future,  higher 
sampling  rates  will  be  used  to  better  explore  higher  instability  frequencies,  since  most  of  the  probes 
used  here  were  found  to  have  a  20dB  cut-off  of  well  beyond  1  MHz. 

Typical  current  (I)  -  voltage  (V)  discharge  characteristics,  as  monitored  by  digital  multi-meters, 
for  the  flow  rates  and  magnetic  field  strength  used  for  these  studies  is  shown  in  Fig.  5.  In  some 
cases,  these  mean  values  were  verified  by  recording  the  discharge  current  oscillations  with  a 
differential  amplifier  placed  across  a  4-ohm  ballast  resistor  in  the  anode  circuit.  At  the  higher 
discharge  voltages  (e.g.,  -200V),  the  probes  were  seen  to  clearly  affect  the  mean  discharge  current, 
perturbing  it  by  as  much  as  20-30%.  During  the  collection  of  the  data,  the  residence  time  of  the 
probes  within  the  channel  was  several  minutes.  Damage  to  the  probes  from  ion  bombardment 
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precluded  the  collection  of  data  beyond  approximately  185V.  While  it  is  difficult  to  say  how  the 
immersion  of  the  probes  into  the  plasma  affected  the  spontaneous  instabilities  detected,  we  discuss 
later  in  the  paper  the  clear  presence  of  waves  that  are  induced  by  the  presence  of  the  probes 
themselves  within  the  discharge  channel. 

ffl.  DATA  ANALYSIS 

The  data  collected  by  a  single  probe  was  cast  to  display  three-dimensional  renderings  of  the 
plasma  density  fluctuations,  as  shown  in  Figs.  6  and  7.  The  axial  variation  in  the  mean  electron 
density  for  three  operating  conditions  has  been  reported  previously  [16],  and  is  provided  again  here 
for  reference  (Fig.  8).  In  Fig.  6,  the  amplitude  of  the  plasma  density  fluctuations  are  displayed 
verses  operating  voltage,  for  four  axial  locations  x  =  12.7  mm,  0  mm,  -12.7  mm,  and  -25.4  mm. 
Prominent  in  these  figures  is  the  presence  of  strong,  relatively  low  frequency  disturbances  in  the 
range  of  5  -  30  kHz,  and  a  broader  range  of  disturbances  towards  higher  frequencies,  at  very  low 
voltages  downstream  of  and  at  the  exit  plane,  and  at  high  voltages  upstream  of  the  exit  plane. 
Similar  spectral  maps  have  been  reported  for  the  behavior  of  discharge  current  oscillations  in  the 
Russian  Stationary  Plasma  Thruster  SPT-70  [15]. 

In  Fig.  7,  the  Fourier  component  amplitudes:  of  the  plasma  density  fluctuations  and  the  ratios  of 
the  amplitudes  to  the  mean  are  displayed  verses  axial  position  for  a  range  of  operating  voltages.  We 
find  that  the  fluctuations  in  the  plasma  density  can  be  as  high  as  the  mean  values.  While  it  is 
apparent  that  the  fluctuations  are  most  intense  near  the  exit  plane,  it  is  noteworthy  that  the  relative 
fluctuations  (amplitude  of  the  fluctuations  divided  by  the  mean  value)  actually  diminish  in  the  region 
between  x  =  -20mm  and  x  =  0  mm  over  a  wide  range  of  operating  conditions.  Note  also  that  the 
broadband  fluctuations  are  concentrated  at  x  =  -10  mm  and  beyond  the  exit  plane  at  low  voltages,  but 
are  damped  in  this  spatial  region  at  high  voltages.  This  dramatic  difference  suggests  that  the  nature 
of  these  high  frequency  disturbances  differ  depending  on  the  operating  voltage. 

A  comparison  of  the  response  of  two  probes  located  at  the  same  axial  position,  but  separated  by 
some  angle  on  the  azimuth,  provides  additional  information  on  the  nature  of  these  disturbances.  A 
cross-spectral  analysis  of  the  signals  from  the  two  probes  provides  a  measure  of  their  coherence, 
suggests  the  direction  of  propagation,  and  can  be  used  to  estimate  the  azimuthal  phase  velocity.  The 
azimuthal  phase  velocity  Vpy.  of  the  disturbances  can  be  obtained  from: 
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VPy{f)  =  2xf 


(  e^R  N 


(2) 


where  o  is  the  angular  probe  separation  (in  radians),  R  is  the  channel  radius,  /  =  ©/ 2n  is  the 
frequency,  and  ^  2  are  the  phase  shifts  of  the  individual  probe  signals,  determined  from: 


^,2(/)  =  tan_1 


(3) 


with  Re  and  Im  the  real  and  imaginary  components  of  the  complex  Fourier  transform  F(f)  =  Re(f)  + 
ilm(f),  respectively. 

The  signals  from  the  two  probes  were  used  to  derive  wave  dispersion  maps,  rendered  as  the  phase 
velocity  verses  wave  frequency.  In  all  cases,  conditional  sampling  is  employed,  to  isolate  only  the 
strongest  correlated  disturbances  detected  by  the  two  probes.  Two  examples  of  dispersion  maps  are 
shown  in  Fig.  9,  for  a  relatively  low  operating  voltage  (100V)  and  an  axial  position  of  x  =  0,  and  for 
a  relatively  high  voltage  (172V)  and  a  location  just  12  mm  beyond  the  exit  plane.  For  these 
dispersion  studies,  the  azimuthal  probe  separation  was  ep  =  30°.  Disturbances  were  considered  only 

if  the  amplitude  of  the  component  on  the  first  probe  was  at  least  three  times  its  mean,  and  that  on  the 
second  probe,  within  its  mean  and  five  times  its  mean.  The  condition  imposed  on  the  second  probe 
allows  for  the  possible  damping  of  the  wave  or  loss  of  intensity  due  to  out-of-plane  propagation. 

Also  drawn  in  the  dispersion  maps  of  Fig.  9  are  lines  corresponding  to  the  azimuthal  component 
phase  velocities  for  the  m=l,  m=4,  and  m  -  12  azimuthal  modes  of  the  cylindrical  annulus.  In 
general,  we  can  express  the  azimuthal  phase  velocity  for  a  tilted  azimuthal  mode: 

v  _  27rcos(a)  f_  iTtcosja)  f  (4) 

Here,  kx  and  ky  are  the  component  wavenumbers  along  the  axial  and  azimuthal  directions, 
respectively.  We  assume  that  ky  =  m/r,  and  m  is  the  integer  azimuthal  mode  number.  For  a  purely 
azimuthal  wave,  the  axial  wave  number  component  is  kx  =  0.  For  m  =  1  and  kx  -  0  condition,  the 
oscillation  wavelength  is  equal  to  the  channel  circumference,  2nr  where  as  it  is  one-fourth  and  one- 
twelfth  the  circumference  for  the  pure  azimuthal  m  =  4  and  m  =  12  modes,  respectively.  In  the  top 
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frame  of  Fig.  9,  we  also  depict  vertical  lines  corresponding  to  possible  m  =  0  modes.  For  m  =  0,  the 
wave  is  longitudinal  and  the  azimuthal  phase  velocity  is  infinite.  In  general,  the  observed  natural 
disturbances  can  have  both  azimuthal  and  axial  components,  and  can  propagate  with  an  angle  a, 
defined  by: 


tana=^L  (5) 

ky 

An  examination  of  the  dispersion  maps,  and  comparison  to  the  lines  associated  with  the  pure 
azimuthal  m  =  1  mode  suggests  that  some  of  these  disturbances  may  be  tilted  out  of  the  azimuthal 
plane  (ie  kx  t  0).  Indeed,  for  the  m  =  1  modes  shown  in  the  top  and  bottom  frames  of  Fig.  9,  a 
propagation  angle  of  15°  and  20°  respectively  has  been  used  to  best  capture  the  trends  seen  for  these 
disturbances.  The  strong  m  =  1  low-frequency  disturbances  at  x  =  0  also  appear,  however,  to  have  a 
range  in  the  possible  propagation  directions,  preferring  lower  propagation  angles  at  the  lower 
frequencies. 

IV.  RESULTS  DISCUSSION 

A  number  of  features  are  identifiable  in  the  dispersion  maps  of  Fig.  9,  attributed  here  to  four  or 
possibly  five  intrinsic  instabilities,  and  an  instability  induced  by  the  probes  themselves.  We  shall 
discuss  them  here,  some  within  the  context  of  what  has  been  seen  experimentally  in  previous  Hall 
discharge  studies. 

The  most  prominent  mode  seen  at  low  voltages  is  that  of  a  relatively  strong,  m  =  1  azimuthally 
propagating  disturbance  of  frequency  in  the  5-10  kHz  range,  and  of  phase  velocity  Vpy  ~  1000  -  4000 
m/s,  nearly  equal  to  the  critical  ionization  velocity  of  the  xenon  propellant.  As  mentioned  above, 
that  the  data  in  this  region  departs  slightly  from  the  superimposed  line  m  =  1  (a  =  15°)  line  suggests 
that  this  disturbance  lacks  a  unique  tilt  in  its  propagation  vector,  and  that  the  propagation  angle  tends 
to  increase  slightly  with  increasing  frequency.  These  waves  are  similar  in  their  properties  to  those 
identified  by  Janes  and  Lowder  [3],  and  others  [4,7],  as  “rotating  spokes”,  and  are  attributed  to  tilted 
ionization  waves  of  still  unknown  origin,  although  the  connection  to  the  critical  ionization  velocity 
may  be  more  than  coincidental  since  an  appropriate  scaling  is  found  when  operating  on  different 
inert  gases  [3],  A  second  prominent  mode  seen  at  low  discharge  voltage  is  that  associated  with 
nearly  purely  m  =  0  axial  modes  at  higher  frequencies,  a  few  of  which  are  identified  on  the  figure  by 
superimposing  vertical  lines  at  characteristic  frequencies.  These  modes  will  be  more  apparent  on 
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dispersion  plots  discussed  later.  While  the  behavior  of  these  modes  are  consistent  with  an  m  =  0 
azimuthal  mode,  they  could  be  higher-order  azimuthal  modes  that  are  propagating  in  the  axial 
direction  (the  azimuthally  spaced  probes  would  not  resolve  this).  These  high  frequency  modes 
correspond  to  the  nearly  constant  frequency  bursts  seen  as  intense  horizontal  striations  in  the  low 
voltage  (86V  and  100V)  spectral  maps  in  Fig.  7.  It  is  noteworthy  that  these  disturbances  are 
relatively  turbulent  and  broadband,  and  they  exhibit  a  curious  power-law  behavior  in  their  spectral 
distribution.  Fig.  10  displays  the  Fourier  transform  of  the  plasma  density  fluctuations  for  the  100V 
and  also  86V  cases,  showing  a  characteristic  -5/3  power  law  power-law  decay  in  the  amplitudes  at 
frequencies  beyond  20  kHz  (indicated  by  the  dashed  line  in  the  figures).  It  is  natural  at  first  to 
suspect  that  these  turbulent,  longitudinal  disturbances  may  correspond  to  the  so-called  “transit  time” 
oscillations  first  identified  by  Morizov  and  colleagues  [4,7],  and  conjectured  to  be  associated  with 
transit  time  of  ions  within  the  acceleration  channel.  We  dismiss  this  possibility  here,  because  we  see 
that  they  are  excited  at  relatively  low  discharge  voltage  (in  the  ionization  branch  of  the  I-V 
characteristic),  unlike  the  broadband  oscillations  described  in  Ref.  4  and  7.  Furthermore,  the 
frequencies  are  much  lower  than  the  inverse  transit  time  of  the  ions  in  the  acceleration  region  of  the 
thruster,  which,  for  our  case,  is  approximately  500  kHz. 

A  third  mode  seen  in  the  low  discharge  voltage  dispersion  plot  in  Fig.  9  appears  actually  as  a 
low- velocity  cut-off  to  the  higher  frequency  longitudinal  modes,  and  is  identified  on  the  frame  with  a 
line  corresponding  to  an  m  =  4  mode.  The  origin  of  this  clear  demarcation  is  not  presently 
understood,  although  we  suspect  that  it  is  associated  with  an  azimuthal  asymmetry  due  to  the  four 
equally  spaced  magnetic  solenoids  that  connect  to  the  front  outer  pole-piece.  The  persistence  of  this 
demarcation  is  seen  to  weaken  with  increasing  voltage,  and  is  non-existent  at  the  higher  discharge 
voltages  studied.  This  result  is  consistent  with  an  examination  of  the  axial  and  azimuthal  variation 
in  the  mean  signal  from  a  single  probe,  which  reflects  the  spatial  variation  in  mean  plasma  density, 
as  seen  in  Fig.  11  for  a  representative  low  voltage  (86V)  and  high  voltage  (150V)  condition.  It  is 
seen  that  the  four-fold  asymmetry  associated  with  the  solenoid  placement  is  more  apparent  in  the 
azimuthal  variation  in  the  plasma  density  at  lower  voltages.  The  interaction  between  the  axial 
disturbances  and  the  m  =  4  azimuthal  asymmetry  is  the  subject  of  ongoing  research  in  our 
laboratory. 

A  fourth  mode  which  appears  at  high  operating  voltages  is  that  of  an  azimuthal  m  =  1  wave, 
tilted  to  approximately  20°,  and  in  the  frequency  range  above  20  kHz  (see  Fig.  9  panel  for  172V 
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condition).  High  frequency  azimuthal  disturbances  (extending  as  far  as  ~  1  MHz)  have  been 
identified  in  the  past  [8]  and  have  been  attributed  to  drift  instabilities  of  “magnetosonic”  waves, 
excited  primarily  in  the  axial  region  where  the  magnetic  field  gradients  have  negative  values  (i.e:,  B 
decreases  with  increasing  distance  from  the  anode).  Simple  theories  have  been  suggested  for  the 
growth  rate  and  the  dispersion  of  these  waves  [8]  based  on  a  multi-fluid  description  of  the 
magnetohydrodynamic  flow,  although  a  rigorous  comparison  to  experimental  measurements  has  not 
yet  been  carried  out.  While  it  is  tempting  to  suggest  that  the  disturbances  at  low  frequencies  seen 
here  may  be  the  band  edge  of  such  drift-instabilities,  our  results  do  not  support  such  an  assignment, 
since  these  azimuthal  waves,  responsible  for  the  high  frequency  disturbances  shown  in  Fig.  7  (for  the 
high  voltage  conditions)  are  not  largely  excited  in  the  region  of  negative  magnetic  field  gradients.  In 
fact,  there  is  a  sudden  damping  of  the  waves  at  approximately  x  =  -5  mm,  beyond  which  there  is  very 
little  high  frequency  activity.  This  location,  which  defines  the  edge  of  the  activity  of  these 
disturbances,  is  about  the  location  of  the  peak  in  the  magnetic  field  (see  Fig.  2).  It  is  also  tempting  to 
suggest  that  the  azimuthal  disturbances  seen  at  higher  discharge  voltages  may  simply  be  an  extension 
of  the  low  frequency  azimuthal  waves  seen  at  lower  voltage.  However,  the  azimuthal  waves  seen  at 
the  lower  voltages  are  predominantly  excited  in  regions  where  the  magnetic  field  gradients  are 
negative,  consistent  with  what  has  been  seen  in  the  literature  [4], 

A  fifth  mode  that  is  apparent  at  high  voltage  conditions  is  that  of  an  exceptionally  strong  low- 
frequency  (-5-10  kHz)  axial  (m  =  0 )  disturbance  that  is  present  within  the  ionization  zone.  This 
disturbance,  which  is  associated  with  very  strong  overall  discharge  current  fluctuations,  is  known  to 
interact  with  the  external  circuitry  driving  the  discharge  [20].  In  the  literature,  this  instability  has 
been  referred  to  as  the  characteristic  “breathing”  mode  [20],  the  “circuit”  instability  [2],  or  the  “loop’ 
instability  [9].  The  origin  of  this  instability  is  now  reasonably  well  understood  and  its  behavior  has 
been  captured  in  many  discharge  simulations  [21-23].  The  instability  is  characteristic  of  a  predator- 
prey  behavior  in  the  neutral  xenon  and  xenon  ion  dynamics,  coupled  by  the  ionization  process  [22], 
It  is  manifested  as  a  disturbance  in  the  balance  between  the  rapid  removal  of  xenon  ions  within  the 
discharge  channel  via  electrostatic  acceleration,  and  the  replenishment  of  the  channel  with  ions 
created  from  the  slow-moving  neutrals.  The  characteristic  frequency  of  this  instability  is  expected  to 
be  close  to  the  inverse  of  the  residence  time  of  the  neutrals  within  the  ionization  zone,  and  has  been 
found  to  scale  appropriately  with  diminishing  discharge  scale  lengths  [24]. 
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A  sixth  mode  that  is  identifiable  in  both  frames  of  Fig.  9  is  that  of  an  m  =  12  azimuthal  wave, 
which  has  a  wavelength  equal  to  that  of  the  separation  between  the  two  probes,  and  is  believed  to  be 
“anchored”  by  the  placement  of  the  probes  within  the  discharge  annulus.  The  mode  number  of  these 
disturbances  is  found  to  track  the  relative  placement  of  the  two  probes  used  to  study  the  intrinsic 
plasma  fluctuations,  as  shown  in  Fig.  12,  for  exit  plane  disturbances  (x  =  0  mm),  and  a  moderate 
operating  voltage  (125  V).  An  intriguing  feature  of  the  plots  in  Fig.  11  is  that  the  placement  of  the 
probes  also  impacts  the  behavior  of  the  other  prominent  m  =  1  mode.  It  is  noteworthy  that  as  the 
mode  number  corresponding  to  the  separation  of  the  probes  diminishes,  there  seems  to  be  an 
increased  “coupling”  of  the  disturbances,  as  seen  by  the  tendency  of  the  points  on  them  =  1  line  to 
“gravitate”  towards  the  line  corresponding  to  the  dispersion  of  the  induced  modes.  The  discovery  of 
this  induced  perturbation  and  the  coupling  of  these  disturbances  to  the  m  =  1  mode  may  have 
applications  in  the  control  of  thruster  oscillations,  and,  quite  possibly,  the  control  of  electron 
transport  if  strongly  correlated  azimuthal  waves  are  largely  responsible  for  the  anomalous  axial 
electron  flow. 

Fig.  13  gives  the  dispersion  plots  obtained  for  some  of  the  operating  conditions  of  the  spectral 
maps  shown  in  Fig.  7.  A  comparison  is  made  between  the  plots  obtained  at  three  axial  locations.  At 
an  axial  location  corresponding  to  the  exit  plane  of  the  discharge  (x  =  0),  we  see  that  high  discharge 
voltage  operation  is  dominated  by  relatively  strong  high  frequency  m  =  1  azimuthal  waves  (tilted 
slightly  with  respect  to  the  azimuthal  plane),  and  low  frequency  longitudinal  waves,  most  likely 
associated  with  the  “breathing”  mode  ionization  instability.  At  lower  voltages,  the  high  frequency 
azimuthal  waves  become  apparently  more  longitudinal  in  direction,  and  there  is  a  growth  of  very 
strong  low-frequency  azimuthal  waves,  associated  with  the  ionization  process.  For  the  low  voltage 
case  shown,  the  azimuthal  asymmetry  associated  with  the  placement  and  separation  of  the  magnetic 
solenoids  plays  a  significant  role  in  the  wave  dynamics  only  beyond  x  =  0  mm,  and  appears  to  have  a 
dominating  effect  on  the  propagation  of  the  disturbances  in  this  region.  A  comparison  of  the 
dispersion  plots  at  the  three  locations  for  the  highest  operating  voltage  shown  confirms  the 
conclusions  that  were  drawn  earlier  that  the  growth  and  propagation  of  these  higher  frequency 
azimuthal  waves  do  not  depend  on  the  sign  of  the  local  gradient  in  the  magnetic  field  in  contradiction 
to  what  has  been  postulated  in  the  literature  for  azimuthal  waves  of  higher  frequency  [4,7].  On  the 
other  hand,  a  comparison  of  the  dispersion  plots  for  these  three  axial  locations  for  the  lowest 


li 


IEEE  Transactions  on  Plasma  Sciences,  Submitted,  August,  2000 


operating  voltage  shown  suggests  that  the  low-frequency  azimuthal  waves  indeed  are  seen  to  be 
active  most  notably  in  the  region  where  the  magnetic  field  gradients  are  negative. 

V.  SUMMARY 

A  detailed  examination  of  the  propagation  behavior  of  low  (2-20  kHz)  and  intermediate- 
frequency  (20-100  kHz)  oscillations  in  closed-electron  drift  Hall  a  discharge  is  presented. 
Conditional  sampling  using  two  electrostatic  probes  collecting  ion  saturation  current  is  used  to 
identify  and  extract  properties  of  correlated  structures  associated  with  the  propagation  of  azimuthal 
and  longitudinal  instabilities  within  the  discharge  channel.  The  azimuthal  component  phase 
velocities  (dispersion  plots)  are  determined  for  a  range  of  wave  frequencies  and  over  characteristic 
regimes  of  operation  of  the  discharge.  A  variety  of  propagation  modes  are  observed  and  analyzed, 
some  of  which  have  been  documented  in  the  previous  literature.  In  addition  to  the  low-frequency 
azimuthal  “spokes”  and  the  so-called  “breathing  mode”  that  is  well  known  to  be  excited  in  these 
discharges,  we  have  identified  the  appearance  longitudinal  disturbances  (that  are  prevalent  at  low 
discharge  voltages)  and  azimuthal  disturbances  (at  higher  discharge  voltages)  in  the  20-  100  kHz 
range  that  do  not  seem  to  have  been  described  previously.  Finally,  we  have  seen  that  the  presence  of 
azimuthal  non-uniformities  due  to  the  location  of  the  magnetic  solenoids,  and,  due  to  the  placement 
of  the  diagnostic  probes,  can  couple  strongly  to  the  intrinsic  natural  instabilities.  Such  methods  of 
coupling  to  the  natural  dynamical  behavior  of  these  discharges  may  be  useful  in  the  control  or 
enhancement  of  turbulence,  which  is  known  to  influence  the  axial  electron  transport. 
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Figure  Captions 

FIG.  1.  Schematic  of  a  typical  Hall  discharge. 

FIG.  2.  Axial  profile  of  radial  magnetic  field  strength  for  the  Hall  discharge  studied. 

FIG.  3.  Schematic  of  the  Langmuir  probes  used  to  detect  plasma  density  fluctuations. 

FIG.  4.  Schematic  of  the  experimental  arrangement  showing  the  location  and  orientation  of  the 
Langmuir  probes. 

FIG.  5.  Typical  current  -  voltage  Hall  discharge  characteristics. 

FIG  6.  Spectral  maps  of  the  amplitude  of  the  fluctuations  in  the  probe  signals  over  discharge 
voltage  for  axial  locations  of  12.7  mm,  0,  -12.7  mm,  and  -25.4  mm  (top  to  bottom). 

FIG.  7.  Spectral  maps  of  the  amplitude  fluctuations  in  probe  signal  (left)  and  ratio  of  the 
amplitude  fluctuations  to  the  mean  signal  (right)  for  discharge  conditions  of  86V,  100V,  128V, 
161V  and  184V  (top  to  bottom). 

FIG.  8.  Axial  variation  in  the  mean  electron  density  [16]. 

FIG.  9.  Typical  experimental  dispersion  plots  derived  from  the  azimuthal  probes  located  at  the 
discharge  exit  (x  =  0  mm)  but  separated  by  30°.  The  top  frame  is  for  a  discharge  voltage  of 
100V,  and  the  bottom  frame  is  for  a  discharge  voltage  of  172V. 

FIG.  10.  Fourier  spectrum  of  the  fluctuations  in  plasma  density  at  the  exit  of  the  Hall  discharge. 
The  top  frame  is  for  a  discharge  voltage  of  86V,  and  the  bottom  frame  is  for  a  discharge  voltage 
of  100V. 

FIG.  11.  Axial  and  azimuthal  variation  in  the  mean  plasma  density  for  86V  and  150V  operating 
conditions.  The  gray  scale  is  such  that  white  corresponds  to  peak  plasma  densities  of 
approximately  3  x  101 '  m'3. 

FIG.  12.  Dispersion  plots  obtained  using  (top  to  bottom)  probe  separations  of  180,120,  90,  and 
60,  degrees. 

FIG.  13.  Dispersion  maps  for  the  discharge  conditions  of  Fig.6,  for  just  two  voltage  conditions 
(128V  -  top  panels,  and  186V  -  bottom  panels)  with  the  probes  located  at  x  =  -12.7  mm  (left 
panels),  x  =  0  mm  (center  panels),  and  x  =  12.7  mm  (right  panels).  The  dashed  line  (— ) 
corresponds  to  the  pure  m  =  12  azimuthal  mode,  the  dotted  line  (...)  to  the  m  =  1  azimuthal 
mode,  thee  dot-dash  line  (vertical  line)  to  the  m  =  0  mode,  and  the  sold  line  to  the  m  =  4 
azimuthal  mode. 
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Using  recent  experimental  data  on  the  time-averaged,  spatially  varying  plasma  properties  within  a 
Hall  discharge  plasma,  we  present  in  this  article,  a  theoretical  study  of  the  response  of  this  plasma 
to  small  (linear)  perturbations  in  its  properties.  As  a  starting  point  for  this  analysis,  we  assume  a 
two-dimensional  fluid  description  that  includes  a  simplified  equation  for  the  electron  energy,  and 
constrain  the  azimuthal  wave  vector  such  that  we  excite  only  the  dominant  (m  =  l)  azimuthal 
modes.  The  growth  rate  and  frequencies  of  predominantly  axial  and  azimuthally  propagating  plasma 
disturbances  are  obtained  by  numerical  solution  of  the  resulting  eigenvalue  problem  under  a 
quasiuniform  plasma  condition,  along  the  entire  discharge  channel.  The  results  identify  the 
persistence  of  a  low  frequency  instability  that  is  associated  with  the  ionization  process,  concentrated 
largely  in  the  vicinity  of  the  exit  plane,  where  the  magnetic  field  is  at  its  maximum  value,  consistent 
with  experimental  observations  for  the  relatively  low  operating  voltages  (—100  V)  considered  in 
this  study.  ©  2001  American  Institute  of  Physics.  [DOI:  10.1063/1.1346656] 


I.  INTRODUCTION 

Closed  electron  drift  thrusters  (CEDTs)  are  crossed-field 
co-axial  Hall  discharge  plasma  sources  presently  considered 
for  use  as  satellite  propulsion  devices.  Their  potential  as  high 
specific  impulse  thrusters  for  space  propulsion  date  back  to 
the  early  1960s.1 2  In  a  traditional  CEDT,  the  plasma  is  sus¬ 
tained  in  an  orthogonal  electric  and  magnetic  field.  Electrons 
emitted  by  an  external  cathode  are  magnetized  by  the  mag¬ 
netic  field  and  drift  in  the  EXB  direction.  Ions,  generated 
through  volume  ionization,  having  a  Hall  parameter  much 
less  than  unity,  are  not  magnetized,  but  are  accelerated  by  the 
electrostatic  fields  established  by  the  retarded  electron  flow. 
In  a  co-axial  geometry,  the  electrons  are  constrained  to  move 
in  the  azimuthal  direction  of  the  closed  EXB  drift  (hence  the 
name  “closed  electron  drift  thruster”).  The  cross-field  drift 
of  electrons  toward  the  anode  is  known  to  be  anomalously 
high  (greater  than  that  determined  by  the  plasma  collisional- 
ity)  and  the  mechanism  for  this  anomalous  transport  is  still 
the  subject  of  ongoing  debate. 

Modem  Hall  discharges  developed  for  space  propulsion 
operate  in  a  quasisteady  mode,  i.e.,  in  some  operating  re¬ 
gimes,  the  discharge  current  can  fluctuate  by  as  much  as  50% 
or  more  of  its  mean  value.3  Furthermore,  they  exhibit  a  rich 
spectrum  of  fluctuations  in  their  plasma  properties,4-8  with 
intense  fluctuations  occurring  predominantly  in  the  region 
where  the  EXB  electron  drift  velocity  is  near  its  maximum.9 
These  plasma  fluctuations  are  believed  to  be  partly  respon¬ 
sible  for  the  anomalous  electron  transport  across  the  mag¬ 
netic  field,3,4,9  and  furthermore,  these  fluctuations  impact  the 
overall  design  of  power  supplies  for  propulsion 
applications.10 


a)  Author  to  whom  correspondence  should  be  addressed;  electronic  mail: 
cap@leland.stanford.edu 


The  characteristic  spectrum  of  fluctuations  in  Hall  dis¬ 
charges  at  lower  frequencies  (^100  kHz)  are  particularly 
sensitive  to  operating  conditions,  such  as  the  magnetic  field 
strength  and  the  applied  discharge  voltage.  In  a  recent  review 
article,  Choueiri  attempted  to  classify  these  disturbances 
when  the  magnetic  field  is  varied.3  A  similar  demarcation 
between  “operating  modes”  is  seen  when  the  magnetic  field 
is  held  constant,  and  the  applied  voltage  is  varied.10,11  Figure 
1  shows  the  variation  in  the  spectral  intensity  of  electron 
density  fluctuations  with  discharge  voltage,  at  the  exit  of  an 
annular  Hall  discharge  channel  (x=0)  where  the  magnetic 
field  is  a  maximum.  These  spectral  maps  were  generated 
with  a  Langmuir  probe  collecting  ion  saturation  current.8,11 
At  very  low  voltages,  when  the  discharge  is  operating  in  the 
so-called  “ionization  regime’  ’  (when  current  rises  sharply 
with  increased  voltage — a  regime  characterized  by  an  in¬ 
creasing  ionization  fraction  with  increased  current),  the  fluc¬ 
tuations  are  predominantly  narrow  band  and  of  relatively  low 
frequency  (5-20  kHz),  although  some  weaker  higher  fre¬ 
quencies  persist  to  beyond  100  kHz.  At  relatively  high  cur¬ 
rents  and  high  voltages,  when  the  current  saturates  with  in¬ 
creasing  voltage  (in  the  so-called  “current  saturation 
regime”),  the  fluctuations  have  similar  spectral  characteris¬ 
tics,  but  the  strong  features  at  low  frequency  are  noticeably 
shifted  from  those  at  lower  voltages.  At  moderate  voltages, 
say  120  V  (in  the  “transition  regime”),  the  discharge 
abruptly  transitions  from  the  ionization  to  current  saturation 
regimes.  Here  we  see  that  the  plasma  is  relatively  quiescent 
at  this  location,  with  the  exception  of  the  strong  disturbance 
at  15-20  kHz. 

The  disturbances  are  known  to  be  quite  dependent  on  the 
location  within  the  plasma  discharge  channel.  Figure  2 
shows  the  variation  in  the  spectral  intensity  of  the  electron 
density  disturbances  with  axial  position  within  the  annular 
discharge  channel,  for  the  same  general  conditions  of  Fig.  1 , 
but  at  a  discharge  voltage  of  100  V.  With  the  exception 
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Discharge  Voltage  (V) 


FIG.  1.  Spectral  map  of  disturbances  in  plasma  density  at  the  discharge  exit 
for  various  operating  voltages.  The  peak  magnetic  field  is  100  G,  and  the 
xenon  mass  flow  rate  is  3  mg/s. 

perhaps  of  the  narrow-band  5-15  kHz  disturbance,  we  see 
that  the  disturbances  are  fairly  broadband,  ranging  from  20 
to  100  kHz,  and  concentrated  largely  near  the  exit  plane,  near 
where  the  radial  component  of  the  magnetic  field  is  stron¬ 
gest,  and  perhaps  in  the  region  where  the  gradient  in  the 
magnetic  field  is  negative.11  These  broadband  disturbances 
are  easily  separated  from  the  lower  frequency  narrower  band 
features,  suggesting  that  these  two  waves  are  of  different 
origins  and  character.  A  further  difference  is  seen  between 
frequency-position  renderings  of  the  fluctuations  for  other 
operating  voltages11  suggesting  that  the  nature  of  the  insta¬ 
bilities  is  quite  sensitive  to  the  regime  of  operation  defined 
by  the  discharge  voltage. 

A  comparison  of  the  response  of  two  probes  located  at 
the  same  axial  position,  but  separated  by  some  angle  on  the 
azimuth,  provides  additional  information  on  the  nature  of 
these  disturbances  Details  of  these  measurements  are  pro¬ 
vided  in  Ref.  11.  A  conditional  cross-spectral  analysis  of  the 
signals  from  the  two  probes  provides  a  measure  of  their  co- 


Probe  Position  Along  Thruster  Channel  (mm) 

FIG.  2.  Spectral  map  of  disturbances  in  plasma  density  at  various  positions 
within  the  channel  (x=0  corresponds  to  the  channel  exit,  and  jc 
=  -  80  mm  corresponds  to  the  location  of  the  anode). 


FIG.  3.  Wave  dispersion  plot  of  plasma  density  disturbances  detected  at  the 
exit  plane  of  the  discharge  channel  (a=0  mm),  for  a  discharge  voltage  of 
100  V,  peak  magnetic  field  of  100  G,  and  a  mass  flow  rate  of  3  mg/s.  There 
is  an  apparent  low  frequency  m=  1  mode  (solid  line),  high  frequency  m 
=  0  modes  (vertical  dashed  line),  an  m-4  cutoff  (dashed  line),  and  an  m 
=  12  probe-induced  mode  (dotted  line). 


herence,  suggests  the  direction  of  propagation,  and  can  be 
used  to  estimate  the  azimuthal  phase  velocity.  Figure  3 
shows  a  typical  dispersion  plot,  for  discharge  conditions  of 
Fig.  2,  with  the  two  probes  separated  by  30°  along  the  azi¬ 
muthal  direction  in  a  plane  located  at  x  —  0  mm. 

Four  characteristic  modes  are  identified  in  the  dispersion 
plot.  The  most  prominent  is  that  of  an  azimuthally  propagat¬ 
ing  mode  (m  =  l  mode  number,  suggesting  an  azimuthal 
wavelength  equal  to  the  circumference  of  the  channel,  as 
indicated  by  the  overlying  solid  line)  of  frequencies  in  the 
5-20  kHz  range,  and  of  phase  velocity  Vp 
~  1000-4000  m/s,  nearly  equal  to  the  critical  ionization  ve¬ 
locity  of  the  xenon  propellant.  The  superimposed  line  is  that 
which  would  be  expected  for  an  azimuthal  m  =  1  mode  with 
an  out  of  plane  tilt  of  15°  to  the  propagation  vector.  The 
spread  of  these  disturbances  above  the  m  =  1  line  indicates 
that  the  tilt  in  the  propagation  direction  is  not  unique  for  all 
frequencies.  These  waves  are  similar  in  behavior  to  those 
identified  by  Janes  and  Lowder2  as  “rotating  spokes,”  and 
are  attributed  to  ionization  waves  of  still  unknown  origin.  A 
second  prominent  mode  is  that  associated  with  nearly  purely 
axial  modes  at  higher  frequency,  one  of  which  is  identified 
on  the  figure  by  superimposing  a  vertical  dashed  line.  While 
the  behavior  of  these  modes  is  consistent  with  an  m  =  0  azi¬ 
muthal  mode,  it  could  be  higher-order  modes  that  are  propa¬ 
gating  in  the  axial  direction  (the  azimuthally  spaced  probes 
would  not  resolve  this).  These  high  frequency  modes  corre¬ 
spond  to  the  nearly  monochromatic  bursts  in  Fig.  2.  It  is 
noteworthy  that  these  disturbances  are  relatively  broadband 
when  averaged  over  time  (i.e.,  they  occur  over  a  broad  range 
of  frequencies),  and  exhibit  a  curious  —5/3  power-law  be¬ 
havior  in  their  spectral  energy  density.  Figure  4  displays  the 
amplitude  of  the  Fourier  transform  of  the  plasma  density 
fluctuations  for  the  100  and  86  V  cases,  showing  the  power- 
law  decay  in  the  amplitudes  at  frequencies  beyond  20  kHz. 

The  third  mode,  which  is  identifiable  in  Fig.  3,  is  that  of 
an  m  =  12  azimuthal  wave  (superimposed  dotted  line),  which 
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FIG.  4.  Fourier  spectra  of  the  plasma  density  fluctuations  for  discharge 
conditions  of  86  V  (top)  and  100  V  (bottom).  The  power-law  behavior 
beyond  20  kHz  is  apparent,  and  illustrated  by  the  superimposed  line. 

has  a  wavelength  equal  to  that  of  the  separation  between  the 
two  probes,  and  is  believed  to  be  “anchored”  by  the  place¬ 
ment  of  the  probes  within  the  discharge  annulus.  The  mode 
number  of  these  disturbances  is  found  to  track  the  relative 
placement  of  the  two  probes  used  to  study  the  intrinsic 
plasma  fluctuations.11  While  the  discovery  of  this  induced 
perturbation  may  have  potential  applications  in  the  control  of 
thruster  oscillations,  we  do  not  consider  it  further  in  this 
analysis,  since  it  is  not  naturally  excited  within  the  discharge. 

Finally,  the  fourth  mode  in  Fig.  3  appears  actually  as  a 
velocity  limit  to  the  higher  frequency  modes.  It  is  identified 
in  the  figure  as  a  dashed  line  corresponding  to  an  m  =  4 
mode.  The  origin  of  this  clear  demarcation  is  not  presently 
understood,  although  we  suspect  that  it  is  associated  with  an 
azimuthal  asymmetry  due  to  the  four  equally  spaced  mag¬ 
netic  solenoids  that  connect  to  the  front  outer  pole  piece.  The 
persistence  of  this  demarcation  weakens  at  higher  discharge 
voltages.11  The  interaction  between  the  axial  disturbances 
and  the  m  —  4  azimuthal  asymmetry  is  the  subject  of  ongoing 
research  in  our  laboratory. 

II.  BACKGROUND 

In  this  article,  we  present  the  results  of  a  study  aimed  at 
developing  an  understanding  of  the  low  frequency  oscillatory 
behavior  of  closed  electron-drift  Hall  discharges.  In  prior 
theoretical  research  on  the  instabilities  of  these  devices,  most 
notably  in  earlier  Russian  literature,5,6  the  emphasis  was 
placed  largely  on  understanding  the  origins  of  higher  fre¬ 
quency  azimuthal  and  longitudinal  disturbances  that  were  ex¬ 
cited  by  magnetic-field  gradient  drifts.  In  almost  all  cases, 
simplifications  were  made  to  the  governing  equations  so  as 


to  derive  closed-form  analytical  expressions  for  the  disper¬ 
sion  relations  for  these  waves  and,  almost  exclusively,  ion¬ 
ization,  which  is  believed  to  be  predominantly  responsible 
for  the  low  frequency  disturbances  that  we  observe  at  low 
discharge  voltages  (e.g.,  80-125  V),  and  is  usually  neglected 
in  these  prior  theoretical  treatments. 

One  motivation  for  this  work  is  to  test  the  ability  of  a 
multifluid  treatment  of  a  Hall  discharge  to  capture  the  oscil¬ 
latory  behavior  seen  in  the  laboratory.  The  success  of  such  a 
test  would  lead  to  added  confidence  in  the  direct  numerical 
solution  of  the  fluid  equations,  or  perhaps  of  a  hybrid  (fluid/ 
particle)  treatment  in  two  or  three  dimensions,12,13  to  further 
our  understanding  of  these  complex  ionized  plasma  flows.  In 
this  article,  we  make  use  of  an  extensive  and  comprehensive 
experimental  mapping  of  the  background  (time-averaged) 
plasma  properties  in  a  Hall  discharge,  which  provides  us 
with  the  necessary  local  background  conditions  for  the  theo¬ 
retical  stability  analysis.14  We  believe  that  at  present,  this  is 
a  better  approach  to  understanding  fluctuation  phenomenon 
in  Hall  discharges,  since  a  satisfactory  theory  for  the  “sta¬ 
tionary"  ’  plasma  properties  do  not  exist,  in  part  because  of 
the  inherent  instabilities  in  these  devices.  We  shall  focus  first 
at  understanding  low  discharge  voltage  conditions  (e.g.,  100 
V)  for  our  particular  Hall  discharge,15-17  where  we  are  on  the 
ionization  branch  of  the  current- voltage  (7-  V)  characteris¬ 
tics.  Although  this  is  not  the  “optimum”  condition  for  op¬ 
erating  a  Hall  discharge  as  a  plasma  propulsion  source,  it 
provides  a  starting  point  for  a  more  complete  analysis  of  the 
entire  7- V  range,  and  is  within  an  operating  regime  where 
we  believe  that  the  instabilities  are  coupled  in  some  way  to 
the  ionization  process. 

III.  FORMULATION  OF  THE  PROBLEM 

The  coordinate  geometry  for  the  analysis  is  selected  so 
that  the  *  axis  is  along  the  thruster  axis  [along  the  direction 
of  the  applied  electric  field  (E0)],  the  y  axis  represents  the 
azimuthal  coordinate,  with  the  azimuthal  electron  drift  veloc¬ 
ity  in  the  negative  y  direction,  and  the  z  axis  is  along  the 
thruster  radius  corresponding  to  the  principle  direction  of  the 
applied  radial  magnetic  field  B0 .  The  plasma  is  assumed  to 
be  uniform  in  the  z  direction,  and  we  restrict  our  attention  to 
motion  in  the  x-y  plane,  rendering  the  problem  to  two  di¬ 
mensions. 

The  plasma  behavior  is  governed  by  a  set  of  13  multi¬ 
fluid  magnetohydrodynamic  equations  that  describe  the  elec¬ 
tron  (subscript  e ),  ion  z,  and  neutral  xenon  n  number  densi¬ 
ties  (ne, rt; ,«„),'  corresponding  velocities  (u,v,w),  the 
electron  temperature  Te ,  the  components  of  the  magnetic 
field  along  z,  Bz ,  and  the  components  of  the  electric  field 
along  x  and  y  ( Ex ,  and  Ey).  These  equations  are  given  below 
in  vector  form  for  completeness: 

Electron  species  conservation 

dne 

—  +  V-(«eu)  =  a,nenn,  (1) 

Ion  species  conservation 

dfii 

—  +  V-(n,v)  =  a!,n,72„,  (2) 
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Neutral  species  conservation 
Bnn 

—  +  V-(n„w)  =  -ajnenn,  (3) 

Electron  momentum  equation 

Bu 

nem  —  -\-nemu-Vu—  —eneE'~enenXB—nem  Fei(u“v) 

~nemvt  n(u-w),  (4) 

Ion  momentum  equation 
B\ 

n  —  +  n  v  *  V v=  en  ,E —  n  jin  vei( v— u) 

at 

~»im  in^Jv-w),  (5) 

Neutral  momentum  equation 
Bw 

n„M  —  +nnMw-'Vyr=  —n„mven(w—u) 

nnTHjnv-m{’W  v),  (6) 

Electron  energy  equation 
3  BT e 

2nekB-^-  =  eneu-'E-ainennei,  (7) 

Ampere’s  law 

1  o*E 

VxB=^i  —  +/M0en(\-u),  (8) 

Faraday’s  law 

<?B 

VXE--— .  (9) 

Here,  at  is  the  volumetric  rate  constant  for  ionization  and 
ven  the  averaged  electron  momentum  transfer  collision  fre¬ 
quencies,  which  are  determined  from  the  known  electron 
temperature,  electron  and, neutral  densities,  and  the  ioniza¬ 
tion  and  momentum  transfer  cross  sections  for  xenon.18  vc[  is 
the  electron-ion  momentum  transfer  collision  frequency, 
calculated  using  the  Spitzer-Harm  expression,19  and  vm  is 
the  ion-neutral  momentum  transfer  collision  frequency,  ap¬ 
proximated  using  a  constant  cross  section  determined  from 
the  molecular  diameter  of  xenon.  The  ionization  rate  and 
collision  frequencies  are  taken  to  be  time  independent,  al¬ 
lowing  a  linearization  of  the  equations  in  plasma  and  neutral 
density,  and  so  the  ionization  instabilities  captured  here  are 
driven  to  first  order  by  fluctuations  only  in  the  plasma  den¬ 
sity.  The  other  variables  in  the  governing  equations  are  the 
electron  charge  and  mass  e  and  m,  the  ion  and  neutral  xenon 
mass  M,  the  reduced  mass  m  the  speed  of  light  c,  Boltz¬ 
mann  constant  kB ,  the  permeability  of  free  space  jjlq,  and 
the  ionization  energy  of  xenon  . 

In  addition  to  the  volumetric  source  of  ionization,  we 
model  the  radial  transport  of  electrons  and  ions  as  driven 
primarily  by  the  loss  of  electrons/ions  at  the  channel  wall.  As 
a  result,  the  electron  and  ion  continuity  equations  are  rewrit¬ 
ten  to  show  an  explicit  dependence  on  the  wall  loss  rate  vwall 


FIG.  5.  Empirically  determined  ionization  rate  constant  and  wall  loss  rates 
for  the  low  voltage  (100  V)  discharge  conditions. 


dne 

Bneux 

Bneuy 

(10a) 

-p+ 

Bx 

dy  ^e^wall? 

Brii 

Bn  iV  x 

Bn-Vy 

(10b) 

Bx 

Qy  ~(Xi^enn  ^e^wall’ 

In  expressing  the  continuity  equations  in  this  way,  we  as¬ 
sume  that  there  is  a  net  loss  of  charge  at  the  wall,  controlled 
largely  by  the  electron  diffusion.  While  the  net  wall  loss  rate 
is  expected  to  depend  in  a  complex  way  on  the  secondary 
electron  emission  process,  we  determine  its  value  experi¬ 
mentally,  from  the  measured  axial  variations  in  the  back¬ 
ground  plasma  properties,  and  from  the  computed  volumetric 
ionization  rate  constant  (assuming  a  Maxwellian  electron  en¬ 
ergy  distribution  at  the  measured  electron  temperature): 

1  dJe 

•W  —  jj-.  (11) 

The  spatial  variation  in  the  ionization  rate  constant  and 
the  wall  loss  rate  are  shown  in  Fig.  5  for  the  100  V  discharge 
condition. 

IV.  STABILITY  ANALYSIS 

In  general,  for  the  case  of  an  initially  axially  nonuniform 
plasma,  a  stability  analysis  would  require  the  solution  of  a 
set  (13  in  this  case)  of  homogeneous  ordinary  differential 
equations  for  the  perturbation  amplitudes,  obtained  by  as¬ 
suming  a  perturbation  of  the  form 

f=Mx)+fi(x,y)e~‘M,  (12) 

where  /  is  the  plasma  property  of  interest,  /0  is  the  steady- 
state  background  property  value  dependent  on  position  for 
the  nonuniform  plasma  case,  and  fx  is  its  small  perturbation 
amplitude  (which  can  also  depend  on  position  and  the  mode 
frequency).  In  this  study,  we  shall  simplify  the  analysis  by 
assuming  that  the  plasma  is  initially  uniform  in  the  azimuthal 
direction  y,  and  that  the  characteristic  length  of  the  axial 
perturbations  ( — Kx ,  the  wavelength  component  along  x)  are 
smaller  than  the  characteristic  length  scale  over  which  the 
background  properties  change,  i.e., 
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(13) 


This  scale  length  is  approximately  the  length  of  the  discharge 
channel  (~8  cm)  but,  in  some  cases,  is  best  characterized  by 
the  length  of  the  ionization  zone  (~2  cm).  Under  these  con¬ 
ditions,  the  problem  reduces  to  a  homogeneous  set  of  linear, 
algebraic  equations  with 


fi(x,y,t)=hei{kxX+kyy- 


(Ot ) 


(14) 


where /^w)  is  the  mode  frequency-dependent  amplitude  of 
the  perturbation.  Although  simplified,  this  first  approach 
leads  to  a  direct  comparison  between  the  predicted  and  ob¬ 
served  azimuthal  and  axial  mode  propagation  behavior  at 
various  locations  within  the  discharge  channel.  The  analysis 
results  in  the  prediction  of  linear  perturbation  amplitudes  and 
frequencies  that  may  depend  on  position.  In  the  analysis, 
gradients  in  these  background  properties  are  evaluated 
explicitly,  based  on  their  local  values,  determined 
experimentally.14 

The  resulting  set  of  linear  homogeneous  equations  is 
(where  we  have  dropped  the  "  symbol  over  the  perturbation 
amplitudes  for  clarity): 
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The  electron  conservation  equation  is  modified  to  be  ex¬ 
plicitly  dependent  on  the  gradient  in  the  magnetic  field  by 
solving  the  electron  momentum  equations  for  the  component 
velocities  (neglecting  collision  terms),  i.e., 

m  dulx  m  duix  Elx  w0v 

“,y=  ~  7B^  IT’  eB^,U°y~dy  HtBlz ' 

(28a) 


m  diiiy  m  du{y  E\y 

Ulx=7B^~dT+7B^U°y~^r+BZ' 


(28b) 


and  by  substituting  these  back  into  the  electron  species  equa¬ 
tion.  In  this  way,  we  are  able  to  capture  any  magnetic  field- 
gradient  driven  modes  in  a  manner  similar  to  that  suggested 
by  Esipchuck  and  Tilinin.6  However,  in  our  case,  we  have 
also  kept  for  now  terms  associated  with  gradients  in  other 
plasma  properties  and  also  terms  associated  with  the  electron 
inertia,  although  these  inertia  terms  are  expected  to  be  im¬ 
portant  for  only  very  high  frequency  modes. 

Our  approach  is  to  solve  the  so-called  “temporal  formu¬ 
lation”  of  the  problem,  where  the  wave  vectors  kx  and  kY  are 
taken  as  real  quantities,  and  solve  the  system  of  equations  for 
the  complex  cd,  searching  through  the  roots  for  positive  fre¬ 
quencies  (real  component  of  o>)  and  positive  growth  rates 
(imaginary  component  of  cd).  For  purely  azimuthal  modes, 
we  constrain  ky=  —m/r  (negative  sign  since  we  are  looking 
for  disturbances  propagating  in  the  negative  azimuthal  direc¬ 
tion),  with  m=l  (to  be  consistent  with  what  is  seen  experi¬ 
mentally)  and  r=4cm,  the  radius  of  the  midpoint  of  the 
annular  discharge  channel.  For  predominantly  azimuthal 
waves  with  small  propagation  angles  a  out  of  the  azimuthal 
plane,  we  use  kx=ky  tan(a)  for  the  out  of  plane  component 
of  the  wave  vector.  For  our  particular  discharge  channel,  we 
find  that  we  satisfy  the  weak  axial  nonuniformity  condition 
above  provided  that  20°  S5  a  5^70°.  For  other  propagation 
angles,  the  axial  component  of  the  wavelength  is  comparable 
to  and  slightly  greater  than  the  length  of  the  channel,  and  so 
the  results  for  these  intermediate  propagation  angles  must  be 
interpreted  cautiously. 

To  find  the  mode  frequencies,  we  solve  a  generalized 
eigenvalue  problem  of  the  form 


1 

tA  —  a)I 

i 


?i  =  0, 


(29a) 


where  the  vector  describing  the  perturbation  amplitudes  is 

fj  \fl  ie  ,Wjj*  , Mi y  \x  9V  \y  }XW 9E \e  , 

Elx,Ely,Blz]  (29b) 

and  all  of  the  elements  of  the  coefficient  matrix  A  are  func¬ 
tions  of  the  axial  position,  known  from  experiments.  Here,  I 
is  the  identity  matrix.  For  each  eigenmode  representing  a 
nondamped  physical  instability,  we  can  determine  the  corre¬ 
sponding  set  of  perturbation  amplitudes  (and  relative  phases) 
as  complex  components  of  that  specific  eigenvector.  It  is 
noteworthy  that  to  cast  the  equation  set  in  this  eigenvalue 
form,  we  had  to  neglect  the  local  acceleration  (inertia  terms) 


in  the  electron  momentum  equation.  We  find  that  for  the  low 
frequency  modes  of  interest  here,  these  terms  do  not  affect 
the  overall  results  of  the  calculation. 


V.  RESULTS 

The  eignevalue  problem  was  applied  to  the  conditions 
for  Hall  discharge  operation  at  low  voltage  (100  V),  as  de¬ 
scribed  in  Ref.  14.  All  of  the  background  plasma  properties 
needed  to  determine  the  coefficient  matrix  A  are  specified  in 
that  article,  for  the  length  of  the  entire  discharge  channel, 
extending  to  about  10  mm  beyond  the  exit  plane.  Apart  from 
the  trivial  solution  (o)  =  0),  we  generally  find  four  modes 
very  close  to  the  electron  plasma  frequency  coe(~10nHz), 
two  modes  very  close  to  the  lower  hybrid  frequency 
(~  106- 107  Hz),  three  modes  between  20  and  100  kHz,  and 
three  other  modes  between  1  and  15  kHz. 

Of  the  lower  frequency  roots,  there  is  usually  only  one 
root  that  satisfies  the  condition  for  unstable  growth  of  the 
disturbance.  In  some  cases,  two  roots  are  found,  and  that 
which  has  the  largest  growth  rate  is  selected  for  display.  The 
frequencies  of  these  roots  are  found  to  depend  strongly  on 
both  the  axial  position  and  on  the  propagation  angle.  A  three- 
dimensional  rendering  of  the  frequency  verses  propagation 
angle  and  axial  position  is  shown  in  Fig.  6(a).  The  corre¬ 
sponding  growth  rate  and  electron  density  amplitude  are 
shown  in  Figs.  6(b)  and  6(c).  It  is  apparent  from  the  figure 
that  the  disturbances  at  higher  frequencies  favor  longitudinal 
propagation,  whereas  azimuthal  propagation  is  largely  asso¬ 
ciated  with  frequencies  below  20  kHz,  in  qualitative  agree¬ 
ment  with  what  is  seen  experimentally.  An  examination  of 
the  amplitudes  and  growth  rate  also  indicate  that  the  distur¬ 
bances  are  concentrated  between  the  axial  locations  x  — 
—  20  and  x~  10  mm,  also  in  agreement  with  the  experimen¬ 
tally  measured  plasma  density  fluctuations  as  seen  in  Fig.  2. 
A  projection  of  all  of  the  low  frequency  points  on  the 
frequency-axial  position  plane  (see  Fig.  7),  shows  that  the 
lower  frequency  azimuthal  waves  extend  throughout  a 
greater  region  of  the  channel,  whereas  the  higher  frequency 
waves  are  concentrated  mostly  in  the  region  between  x  — 
“10mm  and  x  =  0. 

We  find  that  the  gradients  in  the  background  plasma 
properties  have  a  significant  influence  on  the  resulting  mode 
behavior.  In  the  numerical  calculation,  we  can  selectively 
consider  the  gradients  of  certain  properties,  to  see  the  influ¬ 
ence  that  these  gradients  have  on  the  resulting  dynamical 
system.  The  results  of  turning  off  all  of  the  gradients  in  the 
calculation  are  shown  in  Figs.  8  and  9.  Although  the  quali¬ 
tative  appearance  of  the  plots  remain  the  same,  a  detailed 
inspection  reveals  that  the  axial  gradients  in  properties,  all  of 
which  are  positive  between  x  =  —  20  and  x=—5  mm,  pro¬ 
mote  low  and  high  frequency  activity  in  this  region,  whereas 
in  the  region  beyond  this  (between  x  =  —  5  mm  and  x 
— 10  mm),  there  seems  to  be  suppressed  activity.  In  this  re¬ 
gion,  the  plasma  density  and  magnetic  field  have  negative 
axial  gradients,  but  the  ion  velocity  is  strong  and  positive  in 
its  gradient. 

To  further  understand  how  gradients  influence  the  devel¬ 
opment  of  these  instabilities,  the  calculations  were  repeated 
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FIG.  6.  Three-dimensional  renderings  of  the  results  of  the  linear  stability 
analysis  for  the  low  voltage  (100  V)  discharge  conditions:  (a)  computed 
frequencies,  (b)  initial  linear  growth  rates,  and  (c)  electron  density  ampli¬ 
tude. 


with  only  the  plasma  density  and  magnetic  field  gradients 
turned  on.  The  frequencies  are  shown  in  Fig.  10  as  projec¬ 
tions  onto  the  frequency-axial  position  plane  for  direct  com¬ 
parison  to  Figs.  7  and  9.  It  is  quite  apparent  that  the  plasma 
density  and  magnetic  field  gradients  promote  the  growth  of 
the  instabilities  upstream  of  the  location  of  the  peak  in  the 
magnetic  field,  where  these  gradients  have  a  positive  value. 
On  the  other  hand,  they  seem  to  have  little,  if  any  effect  on 
the  frequencies  downstream  of  the  peak  in  the  magnetic 
field.  It  is  noteworthy  that  turning  on  just  the  ion  velocity 
gradient,  which  is  strong  in  the  region  between  jc 
=  —  10  mm  and  *=10  mm,  seems  to  result  in  much  reduced 
activity  overall,  and  is  most  likely  the  cause  of  the  absence 
of  high  frequencies  seen  beyond  the  exit  plane  (*  =  0)  in 
Fig.  7. 


FIG.  7.  Projection  of  the  frequencies  from  Fig.  6(a)  onto  the  frequency-axial 
position  plane. 


(a) 


(c) 


FIG.  S.  Three-dimensional  renderings  as  in  Fig.  6,  but  with  all  of  the  axial 
plasma  property  gradients  turned  off. 
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FIG.  9.  Projection  of  the  frequencies  from  Fig.  8  onto  the  frequency-axial 
position  plane.  All  of  the  gradients  in  plasma  properties  have  been  turned 
off. 

Finally,  to  identify  the  source  of  these  instabilities,  we 
examined  the  sensitivity  of  the  computed  frequencies,  ampli¬ 
tudes,  and  growth  rates  to  the  volume  ionization  and  wall 
loss  source  terms.  The  computed  results  for  the  case  where 
all  of  the  gradients  are  turned  on,  but  with  the  ionization 
process  turned  off,  are  shown  in  Fig.  11.  It  is  apparent  from 
the  figure  that  the  instabilities  nearly  all  but  disappear,  indi¬ 
cating  that  these  instabilities  are  largely  ionization  waves 
driven  in  part  by  the  axial  inhomogeneities  in  the  density  and 
magnetic  field. 

While  this  analysis  does  not  provide  an  analytical  solu¬ 
tion  for  the  dispersion  of  these  ionization  instabilities,  it  pro¬ 
vides  the  first  clear  indication  of  the  importance  of  electro¬ 
thermal  processes  such  as  ionization,  in  the  dynamical 
behavior  of  these  discharges,  at  least  for  the  low  voltage 
operation  examined  here.  It  confirms  the  conjecture  that  the 
low  frequency  azimuthally  propagating  “spokes’  *  are  ioniza¬ 
tion  waves,  as  first  suggested  by  Janes  and  Lowder.2  It  also 
supports  the  conjecture  that  the  higher  frequency  longitudi¬ 
nally  propagating  disturbances  (as  denoted  by  the  m  =  0 
mode  assignments  in  Fig.  3)  are  longitudinal  ionization 
waves.  The  coupling  between  these  longitudinal  disturbances 
and  the  azimuthal  asymmetry  in  plasma  properties  imposed 
by  the  location  of  the  magnetic  solenoids  is  yet  to  be  re¬ 
solved.  Future  theoretical  analyses  will  concentrate  on  devel- 


FIG.  10.  Projection  of  the  calculated  frequencies  onto  the  frequency-axial 
position  plane.  Only  gradients  in  the  radial  magnetic  field  and  plasma  den¬ 
sity  have  been  turned  on. 


xlO4 


<b) 

FIG.  11.  Three-dimensional  renderings  of:  (a)  computed  frequencies  and  (b) 
growth  rate,  with  all  of  the  axial  plasma  property  gradients  turned  on  but 
with  the  electron  and  ion  source  term  (ionization  and  wall  loss)  turned  off. 

oping  a  simple  analytical  understanding  of  these  ionization 
instabilities  now  that  possible  sources,  which  influence  their 
dynamical  behavior,  have  been  partially  clarified. 

VI.  SUMMARY 

The  availability  of  experimental  data  on  the  background 
(time-averaged)  plasma  properties  within  a  Hall  discharge 
has  permitted  a  study  of  the  response  of  this  plasma  to  linear 
perturbations  in  its  properties.  A  two-dimensional  multifluid 
description  accounting  for  electron  heating  is  used  as  the 
basis  for  the  unsteady  analysis.  The  disturbances  are  con¬ 
strained  such  that  only  the  m  —  1  azimuthal  mode  is  excited. 
The  growth  rate  and  frequencies  of  predominantly  axial  and 
azimuthally  propagating  plasma  disturbances  are  obtained  by 
numerical  solution  of  the  resulting  eigenvalue  problem  under 
a  quasiuniform  plasma  condition,  along  the  entire  discharge 
channel.  The  results  identify  the  persistence  of  a  low  fre¬ 
quency  instability  that  is  associated  with  the  ionization  pro¬ 
cess,  concentrated  largely  in  the  vicinity  of  the  exit  plane, 
where  the  magnetic  field  is  at  its  maximum  value,  consistent 
with  experimental  observations  for  the  relatively  low  operat¬ 
ing  voltages  (—100  V)  considered  in  this  study. 
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Abstract.  This  paper  presents  a  discussion  of  the  physics  of  modem  Hall  plasma  thrusters 
and  its  impact  on  the  design  of  new  plasma  thrusters  of  varying  geometry  and  power.  A 
particular  emphasis  is  placed  on  the  design  and  development  of  a  linear-geometry  (non¬ 
coaxial)  source  with  an  open  electron-drift  current.  The  operating  characteristics  of  a 
linear-geometry  Hall  discharge  scaled  to  operate  in  the  50  to  100  W  power  range  are 
presented.  Two  thruster  acceleration  channels  were  fabricated— one  of  alumina  and  one  of 
boron  nitride.  Differences  in  operation  with  the  two  channel  materials  are  attributable  to 
differences  in  the  secondary  electron  emission  properties.  In  either  case,  however,  operation 
is  achieved  despite  the  lack  of  a  closed  electron  current  drift  in  the  Hall  direction,  suggesting 
that  there  is  an  anomalous  axial  electron  mobility,  due  to  either  plasma  fluctuations  or 
collisions  with  the  channel  wall.  Strong  low-frequency  oscillations  in  the  discharge  current, 
associated  with  the  depletion  of  propellant  within  the  discharge,  are  seen  to  appear  and  vary 
with  changes  in  the  applied  magnetic  field  strength.  The  frequency  of  this  oscillatory  mode  is 
higher  than  that  seen  in  larger  (and  higher  power)  discharges,  due  to  the  decreased  residence 
time  of  the  propellant  within  the  channel. 


1.  Introduction 

Hall  discharge  plasma  accelerators  have  been  considered  for 
use  in  satellite  propulsion  since  the  early  1960s  [1,2].  In 
a  Hall  plasma  source,  a  low-pressure  discharge  is  sustained 
within  a  bounded  dielectric  channel  in  crossed  electric  and 
magnetic  fields.  Electrons  emitted  from  a  cathode  external  to 
the  channel,  or  created  by  the  ionization  processes,  drift  along 
the  channel  towards  the  anode  located  at  the  channel  base. 
The  anode  also  serves  as  the  source  of  neutral  propellant 
(typically  xenon).  The  radial  component  to  the  magnetic 
field  is  designed  to  be  a  maximum  near  the  channel  exit, 
and  in  this  region,  the  electrons  become  highly  magnetized, 
as  the  classical  electron  Hall  parameter  is  much  greater 
than  unity.  In  typical  Hall  discharge  plasma  sources,  the 
geometry  is  co-axial  (figure  1)  with  an  annular  channel  and 
surrounding  solenoids  generate  a  radial  magnetic  field.  In 
this  co-axial  configuration,  the  electrons  are  constrained  to 
move  in  the  azimuthal  direction  of  the  closed  E  x  B  drift, 
with  cross-field  drift  providing  the  necessary  electron  current 
to  sustain  the  discharge.  As  the  electron  Hall  parameter 
is  much  greater  than  unity,  the  Hall  current  density  can  be 
many  orders  of  magnitude  greater  than  the  axial  current 
density.  According  to  classical  electron  transport  theory, 
electrons  can  circle  the  annular  channel  in  the  Hall  direction 
many  times  before  being  captured  at  the  anode.  A  co¬ 
axial  geometry  therefore  allows  for  this  ‘closed’  electron 
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drift  in  the  Hall  direction,  and  uninterrupted  Hall  current. 
The  region  of  trapped  electrons  acts  as  a  volumetric  zone  of 
ionization  that  in  some  devices  may  occupy  a  small  fraction 
of  the  overall  channel  length.  The  ions  generated  in  this 
region,  unaffected  by  the  magnetic  field  because  of  their 
large  inertia,  are  accelerated  by  the  electric  field  resulting 
from  the  impeded  electron  flow,  producing  thrust.  Very  high 
ionization  fractions  and  ion  velocities  can  be  generated  with 
these  discharges  and,  due  to  their  high  efficiencies  and  high 
specific  impulse.  Hall  plasma  thrusters  in  the  1-5  kW  power 
range  are  being  evaluated  for  use  on  commercial,  military 
and  research  spacecraft  [3]. 

A  precise  theory  is  lacking  for  the  mechanism  of  cross¬ 
field  electron  transport  in  Hall  plasma  thrusters.  Early 
experiments  on  Hall  plasma  sources  indicated  that  classical 
electron  transport  theory  could  not  account  for  the  measured 
‘anomalous’  axial  (cross-field)  electron  current  densities. 
Janes  and  Lowder  [2]  drew  attention  to  the  presence  of  density 
and  electric  field  fluctuations  within  the  channel  of  a  Hall 
discharge,  and  first  suggested  that  these  plasma  disturbances 
enhance  the  axial  electron  current.  Indirect  measurements  of 
the  ‘effective’  Hall  parameter  as  a  result  of  these  fluctuations 
were  in  agreement  with  the  anomalous  transport  coefficient 
first  identified  by  Bohm  et  al  [4],  which  characterizes 
the  process  now  widely  recognized  as  ‘anomalous’  Bohm 
diffusion  [5].  The  Bohm  mechanism  predicts  an  electron 
mobility  that  scales  inversely  with  the  magnetic  field  strength 
(as  opposed  to  the  classical  B~ 2  scaling),  and  an  effective 
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Figure  1.  Schematic  illustration  of  a  typical  co-axial  Hall 
discharge. 

electron  Hall  parameter  of  about  16.  At  conditions  typical 
of  Hall  plasma  thrusters  near  the  region  where  the  magnetic 
field  is  strongest,  the  classical  Hall  parameter  is  about  500- 
1000.  A  value  of  16  represents  a  significant  enhancement 
in  the  cross-  field  drift,  and  indicates  that  the  ratio  of  Hall 
current  density  to  axial  current  density  may  be  much  less 
than  that  suggested  by  classical  transport  theory.  While 
an  enhanced  electron  current  due  to  fluctuations  is  one 
possible  mechanism  for  enhanced  electron  transport,  the 
operation  of  modem  Hall  plasma  thrusters  seems  to  depend 
significantly  on  the  properties  of  the  dielectric  wall  [6]. 
Previous  researchers  have  proposed  the  possibility  of  an 
enhanced  ‘near-wall  conductivity’  due  to  the  ‘wall  scattering’ 
of  electrons.  While  it  seems  the  precise  knowledge  of 
which  mechanism  is  responsible  for  transport  is  necessary 
to  properly  scale  a  Hall  discharge,  we  show  below  that  either 
of  these  mechanisms  exhibits  the  necessary  dependence  on 
discharge  parameters  to  achieve  a  desired  scaling  in  discharge 
size  or  power. 

Despite  the  progress  that  has  been  made  in  the 
development  of  co-axial  Hall  plasma  thrusters  that  operate  in 
the  kilowatt  power  range,  a  need  has  developed  for  low-thrust, 
high-efficiency  propulsion  devices  to  be  used  for  precise  orbit 
control  on  small,  power-limited  satellites.  A  low-  power 
(10-100  W)  Hall  thruster  could  fill  this  need.  The  proper 
scaling  of  a  Hall  plasma  thruster  for  efficient  operation  at  such 
low  powers  requires  a  renewed  examination  of  the  discharge 
physics  that  controls  thruster  performance.  Also,  alternative 
geometries  that  can  potentially  reduce  thruster  mass  and/or 
size  should  be  investigated.  The  scaling  of  co-axial  Hall 
thrusters  to  lower  powers  has  been  discussed  previously  in 
the  literature  [7, 8].  To  our  knowledge,  however,  no  one 
has  reported  on  the  operation  of  a  Hall  plasma  source  with 
a  linear  geometry  and  hence  an  open  electron-drift.  The 


merits  of  a  linear-  geometry  thruster  are  appealing,  although 
any  discharge  model  based  on  either  classical  or  ‘Bohm’ 
electron  transport  indicates  that  even  for  moderate  aspect 
ratios  (depth  to  channel  length  ratio,  D/L ),  such  a  geometry 
would  interrupt  the  electron  Hall  current.  A  linear  geometry 
allows  compact  packaging  of  the  thruster  in  a  limited  space, 
making  the  magnetic  circuit  amenable  to  the  use  of  permanent 
magnets.  Also,  a  multiple  array  of  linear  thrusters  could  be 
efficiently  stacked  in  order  to  extend  the  operating  envelope 
of  the  propulsion  system.  This  modular  approach  could 
be  used  to  maintain  operation  at  maximum  efficiencies  by 
simply  turning  stacked  low-power  linear  thrusters  on  and  off 
as  needed  to  change  the  thrust  level  rather  than  by  changing 
the  operating  point  of  a  single  thruster. 

In  this  paper,  we  discuss  the  design  and  operation  of 
a  low-power  linear-geometry  Hall  plasma  thruster  based  on 
scaling  arguments  that  we  have  presented  previously  [8].  The 
linear  thruster  tested  here  is  scaled  to  operate  at  a  power  level 
that  is  1CM5%  that  of  co-axial  discharges  built  and  tested  in 
our  laboratory  in  previous  years  [9, 1 0].  The  linear-geometry, 
non-coaxial  Hall  plasma  thruster  has  been  fabricated  and 
operated  at  near-design  conditions.  Operating  characteristics 
are  presented  for  both  alumina  and  boron  nitride  acceleration 
channels  for  a  range  of  peak  magnetic  field  strengths. 

2.  Review  of  Hall  thruster  physics 

Modem  co-axial  Hall  plasma  thrusters  that  operate  in  the  1- 
5  kW  power  range  have  been  shown  to  operate  with  very 
high  thrust  efficiencies — around  50%.  These  thrusters  have 
acceleration  channel  diameters  ranging  from  50  to  280  mm. 
One  feature  common  to  these  thrusters  is  that  the  channel 
width  ( W)  is  approximately  1 5%  of  the  outer  diameter,  which 
itself  is  about  twice  the  acceleration  channel  depth  ( D ).  In 
scaling  these  discharges  to  operate  at  various  power  ranges, 
it  is  often  desirable  to  preserve  the  geometrical  relationship 
between  channel  width,  diameter,  and  depth,  although  the 
physical  basis  for  the  commonly  used  geometrical  parameters 
is  not  well  understood. 

In  a  typical  Hall  thruster,  the  magnetic  field  near  the 
channel  exit  is  sufficient  to  trap  the  electrons  in  cyclotron 
motion.  The  electron  orbit  radius  (Larmor  radius)  is  generally 
smaller  than  the  electron  mean  free  path  X  and  the  acceleration 
channel  width  W.  In  this  way,  the  electrons  are  confined  to 
the  magnetized  portion  of  the  plasma  discharge.  The  Larmor 
radius,  being  dependent  on  particle  mass,  is  much  larger  for 
ions,  so  they  are  largely  unaffected  by  the  magnetic  field.  The 
electron  Larmor  radius,  re,  scales  as: 

r‘~~B r  (1) 

Here  B  is  the  magnetic  field  strength  and  Te  is  the  electron 
temperature.  In  the  design  of  a  low-power  (and  hence 
presumably  smaller)  discharge,  a  decrease  in  W  requires  a 
corresponding  decrease  in  re .  The  magnetic  field  strength 
can  be  tailored  for  proper  scaling;  however,  the  electron 
temperature  is  not  easily  adjusted,  as  it  is  a  consequence 
of  a  more  complex  relationship  between  geometry  and 
operating  conditions.  The  electron  temperature  is  established 
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Abstract.  This  paper  presents  a  discussion  of  the  physics  of  modem  Hall  plasma  thrusters 
and  its  impact  on  the  design  of  new  plasma  thrusters  of  varying  geometry  and  power.  A 
particular  emphasis  is  placed  on  the  design  and  development  of  a  linear-geometry  (non¬ 
coaxial)  source  with  an  open  electron-drift  current.  The  operating  characteristics  of  a 
linear-geometry  Hall  discharge  scaled  to  operate  in  the  50  to  100  W  power  range  are 
presented.  Two  thruster  acceleration  channels  were  fabricated — one  of  alumina  and  one  of 
boron  nitride.  Differences  in  operation  with  the  two  channel  materials  are  attributable  to 
differences  in  the  secondary  electron  emission  properties.  In  either  case,  however,  operation 
is  achieved  despite  the  lack  of  a  closed  electron  current  drift  in  the  Hall  direction,  suggesting 
that  there  is  an  anomalous  axial  electron  mobility,  due  to  either  plasma  fluctuations  or 
collisions  with  the  channel  wall.  Strong  low-frequency  oscillations  in  the  discharge  current, 
associated  with  the  depletion  of  propellant  within  the  discharge,  are  seen  to  appear  and  vary 
with  changes  in  the  applied  magnetic  field  strength.  The  frequency  of  this  oscillatory  mode  is 
higher  than  that  seen  in  larger  (and  higher  power)  discharges,  due  to  the  decreased  residence 
time  of  the  propellant  within  the  channel. 


1.  Introduction 

Hall  discharge  plasma  accelerators  have  been  considered  for 
use  in  satellite  propulsion  since  the  early  1960s  [1,2].  In 
a  Hall  plasma  source,  a  low-pressure  discharge  is  sustained 
within  a  bounded  dielectric  channel  in  crossed  electric  and 
magnetic  fields.  Electrons  emitted  from  a  cathode  external  to 
the  channel,  or  created  by  the  ionization  processes,  drift  along 
the  channel  towards  the  anode  located  at  the  channel  base. 
The  anode  also  serves  as  the  source  of  neutral  propellant 
(typically  xenon).  The  radial  component  to  the  magnetic 
field  is  designed  to  be  a  maximum  near  the  channel  exit, 
and  in  this  region,  the  electrons  become  highly  magnetized, 
as  the  classical  electron  Hall  parameter  is  much  greater 
than  unity.  In  typical  Hall  discharge  plasma  sources,  the 
geometry  is  co-axial  (figure  1)  with  an  annular  channel  and 
surrounding  solenoids  generate  a  radial  magnetic  field.  In 
this  co-axial  configuration,  the  electrons  are  constrained  to 
move  in  the  azimuthal  direction  of  the  closed  E  x  B  drift, 
with  cross-field  drift  providing  the  necessary  electron  current 
to  sustain  the  discharge.  As  the  electron  Hall  parameter 
is  much  greater  than  unity,  the  Hall  current  density  can  be 
many  orders  of  magnitude  greater  than  the  axial  current 
density.  According  to  classical  electron  transport  theory, 
electrons  can  circle  the  annular  channel  in  the  Hall  direction 
many  times  before  being  captured  at  the  anode.  A  co¬ 
axial  geometry  therefore  allows  for  this  ‘closed’  electron 
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drift  in  the  Hall  direction,  and  uninterrupted  Hall  current. 
The  region  of  trapped  electrons  acts  as  a  volumetric  zone  of 
ionization  that  in  some  devices  may  occupy  a  small  fraction 
of  the  overall  channel  length.  The  ions  generated  in  this 
region,  unaffected  by  the  magnetic  field  because  of  their 
large  inertia,  are  accelerated  by  the  electric  field  resulting 
from  the  impeded  electron  flow,  producing  thrust.  Very  high 
ionization  fractions  and  ion  velocities  can  be  generated  with 
these  discharges  and,  due  to  their  high  efficiencies  and  high 
specific  impulse.  Hall  plasma  thrusters  in  the  1-5  kW  power 
range  are  being  evaluated  for  use  on  commercial,  military 
and  research  spacecraft  [3]. 

A  precise  theory  is  lacking  for  the  mechanism  of  cross¬ 
field  electron  transport  in  Hall  plasma  thrusters.  Early 
experiments  on  Hall  plasma  sources  indicated  that  classical 
electron  transport  theory  could  not  account  for  the  measured 
‘anomalous’  axial  (cross-field)  electron  current  densities. 
Janes  and  Lowder  [2]  drew  attention  to  the  presence  of  density 
and  electric  field  fluctuations  within  the  channel  of  a  Hall 
discharge,  and  first  suggested  that  these  plasma  disturbances 
enhance  the  axial  electron  current.  Indirect  measurements  of 
the  ‘effective’  Hall  parameter  as  a  result  of  these  fluctuations 
were  in  agreement  with  the  anomalous  transport  coefficient 
first  identified  by  Bohm  et  al  [4],  which  characterizes 
the  process  now  widely  recognized  as  ‘anomalous’  Bohm 
diffusion  [5].  The  Bohm  mechanism  predicts  an  electron 
mobility  that  scales  inversely  with  the  magnetic  field  strength 
(as  opposed  to  the  classical  B~ 2  scaling),  and  an  effective 
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Figure  1 .  Schematic  illustration  of  a  typical  co-axial  Hall 
discharge. 

electron  Hall  parameter  of  about  16.  At  conditions  typical 
of  Hall  plasma  thrusters  near  the  region  where  the  magnetic 
field  is  strongest,  the  classical  Hall  parameter  is  about  500- 
1000.  A  value  of  16  represents  a  significant  enhancement 
in  the  cross-  field  drift,  and  indicates  that  the  ratio  of  Hall 
current  density  to  axial  current  density  may  be  much  less 
than  that  suggested  by  classical  transport  theory.  While 
an  enhanced  electron  current  due  to  fluctuations  is  one 
possible  mechanism  for  enhanced  electron  transport,  the 
operation  of  modem  Hall  plasma  thrusters  seems  to  depend 
significantly  on  the  properties  of  the  dielectric  wall  [6]. 
Previous  researchers  have  proposed  the  possibility  of  an 
enhanced  ‘near-wall  conductivity’  due  to  the  ‘wall  scattering’ 
of  electrons.  While  it  seems  the  precise  knowledge  of 
which  mechanism  is  responsible  for  transport  is  necessary 
to  properly  scale  a  Hall  discharge,  we  show  below  that  either 
of  these  mechanisms  exhibits  the  necessary  dependence  on 
discharge  parameters  to  achieve  a  desired  scaling  in  discharge 
size  or  power. 

Despite  the  progress  that  has  been  made  in  the 
development  of  co-axial  Hall  plasma  thrusters  that  operate  in 
the  kilowatt  power  range,  a  need  has  developed  for  low-thrust, 
high-efficiency  propulsion  devices  to  be  used  for  precise  orbit 
control  on  small,  power-limited  satellites.  A  low-  power 
(10-100  W)  Hall  thruster  could  fill  this  need.  The  proper 
scaling  of  a  Hall  plasma  thruster  for  efficient  operation  at  such 
low  powers  requires  a  renewed  examination  of  the  discharge 
physics  that  controls  thruster  performance.  Also,  alternative 
geometries  that  can  potentially  reduce  thruster  mass  and/or 
size  should  be  investigated.  The  scaling  of  co-axial  Hall 
thrusters  to  lower  powers  has  been  discussed  previously  in 
the  literature  [7,8].  To  our  knowledge,  however,  no  one 
has  reported  on  the  operation  of  a  Hall  plasma  source  with 
a  linear  geometry  and  hence  an  open  electron-drift.  The 


merits  of  a  linear-  geometry  thruster  are  appealing,  although 
any  discharge  model  based  on  either  classical  or  ‘Bohm’ 
electron  transport  indicates  that  even  for  moderate  aspect 
ratios  (depth  to  channel  length  ratio,  Z)/L),  such  a  geometry 
would  interrupt  the  electron  Hall  current.  A  linear  geometry 
allows  compact  packaging  of  the  thruster  in  a  limited  space, 
making  the  magnetic  circuit  amenable  to  the  use  of  permanent 
magnets.  Also,  a  multiple  array  of  linear  thrusters  could  be 
efficiently  stacked  in  order  to  extend  the  operating  envelope 
of  the  propulsion  system.  This  modular  approach  could 
be  used  to  maintain  operation  at  maximum  efficiencies  by 
simply  turning  stacked  low-power  linear  thrusters  on  and  off 
as  needed  to  change  the  thrust  level  rather  than  by  changing 
the  operating  point  of  a  single  thruster. 

In  this  paper,  we  discuss  the  design  and  operation  of 
a  low-power  linear-geometry  Hall  plasma  thruster  based  on 
scaling  arguments  that  we  have  presented  previously  [8].  The 
linear  thruster  tested  here  is  scaled  to  operate  at  a  power  level 
that  is  10-15%  that  of  co-axial  discharges  built  and  tested  in 
our  laboratory  in  previous  years  [9, 10].  The  linear-geometry, 
non-coaxial  Hall  plasma  thruster  has  been  fabricated  and 
operated  at  near-design  conditions.  Operating  characteristics 
are  presented  for  both  alumina  and  boron  nitride  acceleration 
channels  for  a  range  of  peak  magnetic  field  strengths. 

2.  Review  of  Hall  thruster  physics 

Modem  co-axial  Hall  plasma  thrusters  that  operate  in  the  1- 
5  kW  power  range  have  been  shown  to  operate  with  very 
high  thrust  efficiencies — around  50%.  These  thrusters  have 
acceleration  channel  diameters  ranging  from  50  to  280  mm. 
One  feature  common  to  these  thrusters  is  that  the  channel 
width  ( W)  is  approximately  1 5%  of  the  outer  diameter,  which 
itself  is  about  twice  the  acceleration  channel  depth  (D).  In 
scaling  these  discharges  to  operate  at  various  power  ranges, 
it  is  often  desirable  to  preserve  the  geometrical  relationship 
between  channel  width,  diameter,  and  depth,  although  the 
physical  basis  for  the  commonly  used  geometrical  parameters 
is  not  well  understood. 

In  a  typical  Hall  thruster,  the  magnetic  field  near  the 
channel  exit  is  sufficient  to  trap  the  electrons  in  cyclotron 
motion.  The  electron  orbit  radius  (Larmor  radius)  is  generally 
smaller  than  the  electron  mean  free  path  X  and  the  acceleration 
channel  width  W.  In  this  way,  the  electrons  are  confined  to 
the  magnetized  portion  of  the  plasma  discharge.  The  Larmor 
radius,  being  dependent  on  particle  mass,  is  much  larger  for 
ions,  so  they  are  largely  unaffected  by  the  magnetic  field.  The 
electron  Larmor  radius,  ret  scales  as: 

7*1/2 

"~1r  (1) 

Here  B  is  the  magnetic  field  strength  and  Te  is  the  electron 
temperature.  In  the  design  of  a  low-power  (and  hence 
presumably  smaller)  discharge,  a  decrease  in  W  requires  a 
corresponding  decrease  in  re.  The  magnetic  field  strength 
can  be  tailored  for  proper  scaling;  however,  the  electron 
temperature  is  not  easily  adjusted,  as  it  is  a  consequence 
of  a  more  complex  relationship  between  geometry  and 
operating  conditions.  The  electron  temperature  is  established 
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through  a  balance  between  ohmic  dissipation,  electron- 
particle  collisions  (including  ionization),  and  electron-wall 
collisions.  It  was  decided  instead  to  scale  the  magnetic  field 
strength  as  necessary  and  apply  reasonable  scaling  arguments 
to  preserve  the  mean  electron  energy  from  one  design  to 
another.  It  is  seen  from  (1)  that  if  the  electron  temperature 
is  to  be  preserved  in  the  scaling  to  lower  powers,  reducing 
the  characteristic  size  of  the  thruster  requires  a  concomitant 
increase  in  the  operating  magnetic  field  strength. 

In  a  Hall  discharge’s  use  as  a  propulsion  device,  it  is 
desirable  to  efficiently  utilize  the  propellant,  by  achieving 
as  high  an  ionization  fraction  possible.  In  scaling  a  higher 
power  Hall  discharge  to  lower  powers,  it  is  therefore  desirable 
to  preserve  the  ratio  of  the  characteristic  time  to  ionize  the 
propellant  to  the  residence  time  of  the  propellant  in  the 
discharge  channel.  The  ionization  time  can  be  found  from  the 
inverse  of  the  volumetric  rate  of  ionization  /?,,  which  scales 
linearly  as  the  electron  and  neutral  densities  ( ne  and  na): 

Ri  =  nenaai(Te).  (2) 

Here,  at(Te)  is  the  temperature-dependent  electron  impact 
ionization  rate  coefficient.  The  characteristic  time  for 
ionization  is  r7  =  ne/R\ : 


T i  -  l/naoti.  (3) 

The  residence  time  for  a  neutral  atom  can  be  found  by 
dividing  the  acceleration  channel  depth  D  by  the  velocity 
of  the  neutrals,  so  it  is  expected  to  scale  as: 

r*  ~  D/Ta'/2.  (4) 

Here,  Ta  is  the  neutral  xenon  temperature,  which  is  assumed 
to  be  relatively  uniform,  and  which  will  largely  control  the 
gas  dynamic  behaviour  of  the  neutrals  within  the  channel. 
The  ratio  of  these  two  parameters,  the  ionization  time  over 
the  residence  time,  scales  as 

—  ~  Dna  (5) 

tir 

where  we  assumed  that  the  neutral  xenon  temperature  (along 
with  the  electron  temperature)  is  invariant  to  scale.  As  we 
shall  see,  this  assumption  regarding  the  invariance  in  Ta  may 
be  tenuous,  since  the  xenon  temperature  will  depend  on  the 
anode  and  channel  wall  temperatures,  both  of  which  are  likely 
to  be  considerably  higher  for  a  low-power  device  because 
of  the  geometric  scaling  conclusions  arrived  at  below.  A 
consequence  of  (5)  is  that  a  geometric  reduction  in  the 
channel  depth  requires  a  corresponding  increase  in  the  neutral 
density  to  preserve  the  ratio  of  time  scales.  As  we  shall  see 
from  the  next  section,  this  density  increase  is  achieved  by 
properly  scaling  the  mass  flow  rate  and  the  channel  area. 

The  axial  variation  in  the  magnetic  field  is  also  known  to 
have  a  large  impact  on  discharge  performance.  In  a  modem 
co-axial  Hall  thruster,  the  radial  magnetic  field  is  sharply 
peaked  near  the  exit  of  the  acceleration  channel,  with  a 
distribution  width  that  is  much  less  then  the  channel  depth.  A 
high  magnetic  field  near  the  anode  can  lead  to  a  large  anode 
fall  loss  as  electrons  experience  resistance  to  current  flow. 
Since  magnetic  fields  are  difficult  to  shape,  especially  for  co¬ 
axial  designs,  the  depth  of  the  channel  is  often  dictated  more 


by  the  magnetic  field  distribution  than  geometric  scaling  of 
the  channel  length.  An  advantage  of  a  linear  geometry  over 
a  co-axial  one  is  the  ease  at  which  a  desired  magnetic  field 
distribution  can  be  achieved  with  a  less  complicated  magnetic 
circuit. 


3.  Thruster  scaling  implications 


Based  on  the  physics  presented  in  the  previous  section,  the 
scaling  of  the  discharge  is  relatively  straightforward.  We 
treat  the  desired  discharge  voltage  fa  as  a  design  parameter, 
as  it  directly  determines  the  ion  velocity  (and  hence  specific 
impulse  of  the  thruster),  which  is  often  dictated  by  the  satellite 
mission  objectives. 

We  continue  with  the  assumption  that  the  electron 
temperature  can  be  preserved  with  proper  scaling.  This  is 
justified  if  we  can  argue  that  for  a  reduction  of  the  total 
power  by  some  factor  £,  the  rates  of  energy  loss  and  thrust 
power  are  correspondingly  reduced  by  the  same  factor.  The 
reduction  in  discharge  power  without  a  reduction  in  the 
discharge  voltage  implies  a  reduction  in  the  overall  discharge 
current.  However,  for  proper  geometric  scaling,  the  area 
is  correspondingly  reduced  by  the  factor  £2,  so  the  current 
densities  must  be  increased  by  the  factor  l/£.  The  necessary 
scaling  in  the  ion  current  density  (and  hence  thrust  power)  is 
achieved  if  the  plasma  density  is  correspondingly  increased, 
since  the  velocity  is  unchanged.  The  necessary  scaling  in 
the  axial  electron  current  density  is  achieved  if  the  axial 
electron  drift  velocity,  Vedj  is  arguably  scale-invariant.  As 
we  have  discussed  previously  [8],  both  the  anomalous  Bohm 
transport  and  wall  collisions  will  give  rise  to  drift  velocities 
that  are  scale-  invariant.  The  axial  drift  velocity  associated 
with  Bohm  transport  is  determined  by  the  ratio  of  the  electric 
field  strength,  £,  to  the  magnetic  field  strength 

€  E 

Ved  Bohm  =  ~  -  (6) 

which  will  be  preserved  through  a  geometric  scaling.  If  the 
cross-field  transport  is  largely  controlled  by  wall  collisions, 
then,  for  highly  magnetized  electrons  ( coce  =  eB/me  » 
vwau  =  Ce/W ,  the  wall  scattering  frequency),  the  axial 
electron  drift  velocity  is  approximately 


Ved  Wall  ” 


eEvwaii 

mea?ce 


EW 


(7) 


which  will  also  be  preserved  with  the  proper  geometric 
scaling  since  the  magnetic  field  scales  as  B  1  /  W,  as  discussed 
earlier  (here,  Ce  is  the  mean  thermal  electron  speed,  which 
is  preserved  if  the  temperature  is  preserved,  and  e  and  me 
are  the  charge  and  mass  of  the  electron,  respectively).  The 
increased  electron  number  density  (by  the  factor  1/f)  is 
achieved  because  the  corresponding  decrease  in  the  mass  flow 
rate  results  in  an  increase  in  since  the  area  is  decreased  by 
the  factor  £2 .  This  relies  on  the  assumption  that  the  ionization 
fraction  is  preserved,  which  is  reasonable  if  the  ratio  of  time 
scales  presented  in  (5)  is  also  preserved. 

Finally,  in  order  to  preserve  the  electron  temperature,  we 
must  argue  that  the  electron  energy  loss  rates  will  also  scale 
in  proportion  to  the  decrease  in  power.  It  is  easily  shown  that 
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the  necessary  scaling  is  obtained  if  the  dominant  energy  loss 
mechanism  is  through  wall  collisions.  It  is  noteworthy  that 
volumetric  ionization  will  also  satisfy  the  scaling  condition, 
since  the  energy  loss  rate  through  ionization  is 

Ei  =  nenaoL\Vce  -  g.  (8) 

Here,  Vc  is  the  channel  volume  and  the  ionization  energy 
of  xenon. 

One  undesirable  consequence  of  the  geometric  scaling 
for  operation  at  reduced  power  levels  is  an  increase  in  heat 
flux  to  the  channel  walls  [7].  Since  the  power  is  reduced  by 
the  scaling  factor  f  and  the  wall  area  reduced  by  £2,  the  heat 
flux  to  the  walls  will  increase  by  a  factor  of 1/f .  This  scaling 
consequence  may  prove  problematic  for  very  low-power  (and 
consequently  reduced  size)  Hall  plasma  thrusters. 

It  is  noteworthy  that  the  decreased  residence  time  of  the 
neutral  xenon  in  the  channel  (see  (4))  should  result  in  a  shift 
to  high  frequencies  in  the  characteristic  breathing  instability 
often  seen  in  the  7-10  kHz  frequency  range  in  higher  power 
devices  [11].  While  we  predicted  this  shift  in  an  earlier 
paper  [8],  the  poor  performance  of  the  low-power  thruster 
presented  in  that  study  precluded  such  a  measurement.  We 
have  characterized  the  current  oscillations  in  the  discharge 
reported  on  here,  and  have  found  the  oscillation  frequencies 
to  be  consistent  with  the  nearly  1/10  scaling  carried  out  in 
this  study. 

In  summary,  if  it  is  desired  to  scale  the  power  of  a  Hall 
thruster  by  some  arbitrary  factor  £,  then  the  characteristic 
scale  lengths  of  the  thruster  and  mass  flow  rates  should  be 
scaled  by  the  same  factor,  £.  The  appropriate  adjustment 
to  the  magnetic  field  (preserving  its  shape)  is  to  increase  it 
by  the  factor  1  /£ .  With  these  scaling  laws,  according  to  the 
arguments  presented  above,  the  electron  temperature  should 
be  preserved,  as  well  as  the  ratio  of  electron  current  to  ion 
current. 

4.  Linear  geometry  implications 

In  a  weakly  collisional  steady- state  plasma,  where  the 
electron  Hall  parameter  satisfies  the  condition 

COceTe  »  1  (9) 

the  ratio  of  the  cross-field  (axial)  electron  current  to  the  Hall 
current  is 

(10) 

J eH  0)ceXe 

Here,  re  is  the  time  between  electron  collisions.  If  we  use 
the  classical  electron  collision  time  in  (9),  we  would  find 
that  for  most  modem  thrusters,  the  resulting  Hall  parameter 
is  typically  in  the  range  of  100-1000.  Note  that  the  current 
ratio  described  in  (10)  is  a  scale  invariant  in  that  the  scaling 
laws  introduced  here  would  increase  the  electron  cyclotron 
frequency  in  proportion  to  the  decrease  in  the  electron 
collision  time.  It  is  precisely  this  vast  inequality  between  the 
axial  and  Hall  currents  which  prompted  the  use  of  a  co-axial 
design  in  early  thrusters  since,  as  mentioned  above,  a  co-axial 
geometry  with  a  closed  electron  drift  allows  the  electrons  to 
traverse  the  annulus  many  times  prior  to  anode  capture. 


The  presence  of  an  anomalous  electron  transport 
mechanism,  whether  fluctuation  or  possibly  wall-scattering 
induced,  reduces  the  demand  placed  on  the  ratio  of  the  Hall  to 
axial  electron  current.  A  value  of  16  for  the  ‘effective'  Hall 
parameter,  as  suggested  by  the  anomalous  Bohm  mobility, 
still  implies  an  electron  drift  direction  that  is  predominantly 
in  the  direction  of  the  crossed  electric  and  magnetic  field. 
However,  we  note  that  the  value  of  1 6  for  the  Bohm  coefficient 
is  strictly  speculative,  as  the  effective  Hall  parameters  in 
modem  Hall  thrusters  have  not  been  accurately  characterized, 
and  coefficients  within  a  factor  of  two  or  three  of  this  value 
have  been  obtained  for  other  plasma  devices  [5].  It  is 
therefore  conceivable  that  the  Bohm  coefficient  can  be  less 
than  this  value.  If  so,  then  the  necessity  for  a  closed  electron 
drift  is  removed,  and  with  an  adequate  aspect  ratio  (ratio 
of  channel  length  to  channel  depth),  a  linear  Hall  thruster 
with  an  open  electron  drift  may  perform  equally  well  in 
comparison  to  closed-drift  designs.  However,  even  with 
an  effective  Hall  parameter  of  unity,  the  linear  design  does 
impose  an  asymmetry  in  the  electron  flow,  giving  rise  to 
expected  asymmetric  current  densities  within  the  channel  that 
may  impact  on  discharge  performance. 

5.  Experiment 

5.1.  Test  facility 

The  Stanford  high  vacuum  test  facility  has  been  discussed 
extensively  elsewhere  [8-10].  It  consists  of  a  non-magnetic 
stainless  steel  tank  approximately  1  m  in  diameter  and  1 .5  m 
in  length.  The  facility  is  pumped  by  two  50  cm  diffusion 
pumps,  backed  by  a  425  1  s"1  mechanical  pump.  The 
base  pressure  of  the  facility  is  approximately  10-6  Torr 
as  measured  by  an  ionization  gauge  uncorrected  for  mass 
species.  Thruster  testing  at  xenon  flow  rates  of  2-5  seem 
results  in  chamber  background  pressures  in  the  region  of 
4  x  10“5  Torr.  This  indicates  that  the  facility  has  a  xenon 
gas  pumping  speed  of  around  2000  Is-1.  Propellant  flow 
to  the  thruster  anode  and  cathode  is  controlled  by  two 
Unit  Instruments  1200  series  mass  flow  controllers  factory 
calibrated  for  xenon.  The  propellant  used  in  this  study  was 
research  grade  (99.99%)  xenon. 

5.2.  Linear  Hall  thruster 

The  design  of  the  linear  Hall  thruster  studied  here  is  based  on 
the  scaling  of  a  co-axial  reference  thruster  recently  built  by 
our  laboratory  and  operated  at  a  nominal  power  of 400-700  W 
[12].  A  scaling  factor  of  £  =  0.1  was  used  in  accordance  with 
the  scaling  laws  presented  in  the  previous  sections,  although 
the  performance  of  the  magnetic  circuit  precluded  the  use  of 
a  channel  depth  that  was  one-tenth  the  depth  of  the  reference 
coaxial  discharge.  The  channel  depth  deviated  from  strict 
scaling  laws  in  order  to  reduce  the  magnetic  field  strength  at 
the  anode,  and  hence  the  anode  fall  losses.  A  comparison  of 
the  coaxial  thruster  used  for  scaling  and  the  linear  thruster  is 
shown  in  table  1 .  A  schematic  of  the  linear  Hall  thruster  is 
shown  in  figure  2. 

The  magnetic  circuit  includes  four  90  mm  long 
electromagnet  windings  consisting  of  a  9.5  mm  diameter 
core  of  commercially  pure  iron  with  6  layers  of  22  gauge 
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through  a  balance  between  ohmic  dissipation,  electron- 
particle  collisions  (including  ionization),  and  electron-wall 
collisions.  It  was  decided  instead  to  scale  the  magnetic  field 
strength  as  necessary  and  apply  reasonable  scaling  arguments 
to  preserve  the  mean  electron  energy  from  one  design  to 
another.  It  is  seen  from  (1)  that  if  the  electron  temperature 
is  to  be  preserved  in  the  scaling  to  lower  powers,  reducing 
the  characteristic  size  of  the  thruster  requires  a  concomitant 
increase  in  the  operating  magnetic  field  strength. 

In  a  Hall  discharge’s  use  as  a  propulsion  device,  it  is 
desirable  to  efficiently  utilize  the  propellant,  by  achieving 
as  high  an  ionization  fraction  possible.  In  scaling  a  higher 
power  Hall  discharge  to  lower  powers,  it  is  therefore  desirable 
to  preserve  the  ratio  of  the  characteristic  time  to  ionize  the 
propellant  to  the  residence  time  of  the  propellant  in  the 
discharge  channel.  The  ionization  time  can  be  found  from  the 
inverse  of  the  volumetric  rate  of  ionization  which  scales 
linearly  as  the  electron  and  neutral  densities  ( ne  and  na ): 

Ri  = nenaoii(Te ).  (2) 

Here,  at(Te)  is  the  temperature-dependent  electron  impact 
ionization  rate  coefficient.  The  characteristic  time  for 
ionization  is  t/  =  //?,•: 

r i  =  l/nacci-  (3) 

The  residence  time  for  a  neutral  atom  can  be  found  by 
dividing  the  acceleration  channel  depth  D  by  the  velocity 
of  the  neutrals,  so  it  is  expected  to  scale  as: 

t  r-D/T'J2.  (4) 

Here,  Ta  is  the  neutral  xenon  temperature,  which  is  assumed 
to  be  relatively  uniform,  and  which  will  largely  control  the 
gas  dynamic  behaviour  of  the  neutrals  within  the  channel. 
The  ratio  of  these  two  parameters,  the  ionization  time  over 
the  residence  time,  scales  as 

^  ~  Dna  (5) 

t / 

where  we  assumed  that  the  neutral  xenon  temperature  (along 
with  the  electron  temperature)  is  invariant  to  scale.  As  we 
shall  see,  this  assumption  regarding  the  invariance  in  Ta  may 
be  tenuous,  since  the  xenon  temperature  will  depend  on  the 
anode  and  channel  wall  temperatures,  both  of  which  are  likely 
to  be  considerably  higher  for  a  low-power  device  because 
of  the  geometric  scaling  conclusions  arrived  at  below.  A 
consequence  of  (5)  is  that  a  geometric  reduction  in  the 
channel  depth  requires  a  corresponding  increase  in  the  neutral 
density  to  preserve  the  ratio  of  time  scales.  As  we  shall  see 
from  the  next  section,  this  density  increase  is  achieved  by 
properly  scaling  the  mass  flow  rate  and  the  channel  area. 

The  axial  variation  in  the  magnetic  field  is  also  known  to 
have  a  large  impact  on  discharge  performance.  In  a  modem 
co-axial  Hall  thruster,  the  radial  magnetic  field  is  sharply 
peaked  near  the  exit  of  the  acceleration  channel,  with  a 
distribution  width  that  is  much  less  then  the  channel  depth.  A 
high  magnetic  field  near  the  anode  can  lead  to  a  large  anode 
fall  loss  as  electrons  experience  resistance  to  current  flow. 
Since  magnetic  fields  are  difficult  to  shape,  especially  for  co¬ 
axial  designs,  the  depth  of  the  channel  is  often  dictated  more 
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by  the  magnetic  field  distribution  than  geometric  scaling  of 
the  channel  length.  An  advantage  of  a  linear  geometry  over 
a  co-axial  one  is  the  ease  at  which  a  desired  magnetic  field 
distribution  can  be  achieved  with  a  less  complicated  magnetic 
circuit. 


3.  Thruster  scaling  implications 


Based  on  the  physics  presented  in  the  previous  section,  the 
scaling  of  the  discharge  is  relatively  straightforward.  We 
treat  the  desired  discharge  voltage  <pd  as  a  design  parameter, 
as  it  directly  determines  the  ion  velocity  (and  hence  specific 
impulse  of  the  thruster),  which  is  often  dictated  by  the  satellite 
mission  objectives. 

We  continue  with  the  assumption  that  the  electron 
temperature  can  be  preserved  with  proper  scaling.  This  is 
justified  if  we  can  argue  that  for  a  reduction  of  the  total 
power  by  some  factor  £,  the  rates  of  energy  loss  and  thrust 
power  are  correspondingly  reduced  by  the  same  factor.  The 
reduction  in  discharge  power  without  a  reduction  in  the 
discharge  voltage  implies  a  reduction  in  the  overall  discharge 
current.  However,  for  proper  geometric  scaling,  the  area 
is  correspondingly  reduced  by  the  factor  £2,  so  the  current 
densities  must  be  increased  by  the  factor  l/£.  The  necessary 
scaling  in  the  ion  current  density  (and  hence  thrust  power)  is 
achieved  if  the  plasma  density  is  correspondingly  increased, 
since  the  velocity  is  unchanged.  The  necessary  scaling  in 
the  axial  electron  current  density  is  achieved  if  the  axial 
electron  drift  velocity,  Ved,  is  arguably  scale-invariant.  As 
we  have  discussed  previously  [8],  both  the  anomalous  Bohm 
transport  and  wall  collisions  will  give  rise  to  drift  velocities 
that  are  scale-  invariant.  The  axial  drift  velocity  associated 
with  Bohm  transport  is  determined  by  the  ratio  of  the  electric 
field  strength,  E ,  to  the  magnetic  field  strength 

Ved  Bohm  =  E  ~  -  (6) 

which  will  be  preserved  through  a  geometric  scaling.  If  the 
cross-field  transport  is  largely  controlled  by  wall  collisions, 
then,  for  highly  magnetized  electrons  ( o>ce  =  eB/me  ;$> 
Vwaii  =  Ce/w,  the  wall  scattering  frequency),  the  axial 
electron  drift  velocity  is  approximately 


Ved  Wall  = 


eEvwaii 
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(7) 


which  will  also  be  preserved  with  the  proper  geometric 
scaling  since  the  magnetic  field  scales  as  B  1/  W,  as  discussed 
earlier  (here,  Ce  is  the  mean  thermal  electron  speed,  which 
is  preserved  if  the  temperature  is  preserved,  and  e  and  me 
are  the  charge  and  mass  of  the  electron,  respectively).  The 
increased  electron  number  density  (by  the  factor  l/£)  is 
achieved  because  the  corresponding  decrease  in  the  mass  flow 
rate  results  in  an  increase  in  na,  since  the  area  is  decreased  by 
the  factor  £2 .  This  relies  on  the  assumption  that  the  ionization 
fraction  is  preserved,  which  is  reasonable  if  the  ratio  of  time 
scales  presented  in  (5)  is  also  preserved. 

Finally,  in  order  to  preserve  the  electron  temperature,  we 
must  argue  that  the  electron  energy  loss  rates  will  also  scale 
in  proportion  to  the  decrease  in  power.  It  is  easily  shown  that 
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the  necessary  scaling  is  obtained  if  the  dominant  energy  loss 
mechanism  is  through  wall  collisions.  It  is  noteworthy  that 
volumetric  ionization  will  also  satisfy  the  scaling  condition, 
since  the  energy  loss  rate  through  ionization  is 

Ej  =  nenadi  Vcs  ~  g.  (8) 

Here,  Vc  is  the  channel  volume  and  the  ionization  energy 
of  xenon. 

One  undesirable  consequence  of  the  geometric  scaling 
for  operation  at  reduced  power  levels  is  an  increase  in  heat 
flux  to  the  channel  walls  [7].  Since  the  power  is  reduced  by 
the  scaling  factor  ?  and  the  wall  area  reduced  by  £2,  the  heat 
flux  to  the  walls  will  increase  by  a  factor  of  l/£.  This  scaling 
consequence  may  prove  problematic  for  very  low-power  (and 
consequently  reduced  size)  Hall  plasma  thrusters. 

It  is  noteworthy  that  the  decreased  residence  time  of  the 
neutral  xenon  in  the  channel  (see  (4))  should  result  in  a  shift 
to  high  frequencies  in  the  characteristic  breathing  instability 
often  seen  in  the  7-10  kHz  frequency  range  in  higher  power 
devices  [11].  While  we  predicted  this  shift  in  an  earlier 
paper  [8],  the  poor  performance  of  the  low-power  thruster 
presented  in  that  study  precluded  such  a  measurement.  We 
have  characterized  the  current  oscillations  in  the  discharge 
reported  on  here,  and  have  found  the  oscillation  frequencies 
to  be  consistent  with  the  nearly  1/10  scaling  carried  out  in 
this  study. 

In  summary,  if  it  is  desired  to  scale  the  power  of  a  Hall 
thruster  by  some  arbitrary  factor  £,  then  the  characteristic 
scale  lengths  of  the  thruster  and  mass  flow  rates  should  be 
scaled  by  the  same  factor,  £.  The  appropriate  adjustment 
to  the  magnetic  field  (preserving  its  shape)  is  to  increase  it 
by  the  factor  l/£.  With  these  scaling  laws,  according  to  the 
arguments  presented  above,  the  electron  temperature  should 
be  preserved,  as  well  as  the  ratio  of  electron  current  to  ion 
current. 

4.  Linear  geometry  implications 

In  a  weakly  collisional  steady-state  plasma,  where  the 
electron  Hall  parameter  satisfies  the  condition 

CDceTe  »  1  (9) 

the  ratio  of  the  cross-field  (axial)  electron  current  to  the  Hall 
current  is 


Here,  xe  is  the  time  between  electron  collisions.  If  we  use 
the  classical  electron  collision  time  in  (9),  we  would  find 
that  for  most  modem  thrusters,  the  resulting  Hall  parameter 
is  typically  in  the  range  of  100-1000.  Note  that  the  current 
ratio  described  in  (10)  is  a  scale  invariant  in  that  the  scaling 
laws  introduced  here  would  increase  the  electron  cyclotron 
frequency  in  proportion  to  the  decrease  in  the  electron 
collision  time.  It  is  precisely  this  vast  inequality  between  the 
axial  and  Hall  currents  which  prompted  the  use  of  a  co-axial 
design  in  early  thrusters  since,  as  mentioned  above,  a  co-axial 
geometry  with  a  closed  electron  drift  allows  the  electrons  to 
traverse  the  annulus  many  times  prior  to  anode  capture. 


The  presence  of  an  anomalous  electron  transport 
mechanism,  whether  fluctuation  or  possibly  wall-scattering 
induced,  reduces  the  demand  placed  on  the  ratio  of  the  Hall  to 
axial  electron  current.  A  value  of  16  for  the  ‘effective’  Hall 
parameter,  as  suggested  by  the  anomalous  Bohm  mobility, 
still  implies  an  electron  drift  direction  that  is  predominantly 
in  the  direction  of  the  crossed  electric  and  magnetic  field. 
However,  we  note  that  the  value  of  1 6  for  the  Bohm  coefficient 
is  strictly  speculative,  as  the  effective  Hall  parameters  in 
modem  Hall  thrusters  have  not  been  accurately  characterized, 
and  coefficients  within  a  factor  of  two  or  three  of  this  value 
have  been  obtained  for  other  plasma  devices  [5].  It  is 
therefore  conceivable  that  the  Bohm  coefficient  can  be  less 
than  this  value.  If  so,  then  the  necessity  for  a  closed  electron 
drift  is  removed,  and  with  an  adequate  aspect  ratio  (ratio 
of  channel  length  to  channel  depth),  a  linear  Hall  thruster 
with  an  open  electron  drift  may  perform  equally  well  in 
comparison  to  closed-drift  designs.  However,  even  with 
an  effective  Hall  parameter  of  unity,  the  linear  design  does 
impose  an  asymmetry  in  the  electron  flow,  giving  rise  to 
expected  asymmetric  current  densities  within  the  channel  that 
may  impact  on  discharge  performance. 

5.  Experiment 

5.1.  Test  facility 

The  Stanford  high  vacuum  test  facility  has  been  discussed 
extensively  elsewhere  [8-10].  It  consists  of  a  non-magnetic 
stainless  steel  tank  approximately  1  m  in  diameter  and  1 .5  m 
in  length.  The  facility  is  pumped  by  two  50  cm  diffusion 
pumps,  backed  by  a  425  1  s-1  mechanical  pump.  The 
base  pressure  of  the  facility  is  approximately  10-6  Torr 
as  measured  by  an  ionization  gauge  uncorrected  for  mass 
species.  Thruster  testing  at  xenon  flow  rates  of  2-5  seem 
results  in  chamber  background  pressures  in  the  region  of 
4  x  10”5  Torr.  This  indicates  that  the  facility  has  a  xenon 
gas  pumping  speed  of  around  2000  1  s-1.  Propellant  flow 
to  the  thruster  anode  and  cathode  is  controlled  by  two 
Unit  Instruments  1200  series  mass  flow  controllers  factory 
calibrated  for  xenon.  The  propellant  used  in  this  study  was 
research  grade  (99.99%)  xenon. 

5.2.  Linear  Hall  thruster 

The  design  of  the  linear  Hall  thruster  studied  here  is  based  on 
the  scaling  of  a  co-axial  reference  thruster  recently  built  by 
our  laboratory  and  operated  at  a  nominal  power  of 400-700  W 
[12].  A  scaling  factor  of  £  =0.1  was  used  in  accordance  with 
the  scaling  laws  presented  in  the  previous  sections,  although 
the  performance  of  the  magnetic  circuit  precluded  the  use  of 
a  channel  depth  that  was  one-tenth  the  depth  of  the  reference 
coaxial  discharge.  The  channel  depth  deviated  from  strict 
scaling  laws  in  order  to  reduce  the  magnetic  field  strength  at 
the  anode,  and  hence  the  anode  fall  losses.  A  comparison  of 
the  coaxial  thruster  used  for  scaling  and  the  linear  thruster  is 
shown  in  table  1.  A  schematic  of  the  linear  Hall  thruster  is 
shown  in  figure  2. 

The  magnetic  circuit  includes  four  90  mm  long 
electromagnet  windings  consisting  of  a  9.5  mm  diameter 
core  of  commercially  pure  iron  with  6  layers  of  22  gauge 
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Linear-geometry  Hall  plasma  source  with  an  open  electron-drift 


Figure  2.  Schematic  of  the  linear-geometry  open  electron-drift 
Hall  discharge. 


Figure  3.  Photograph  of  the  linear-geometry  open  electron-drift 
Hall  discharge. 


insulated  copper  magnet  wire.  The  magnetic  bottom  plate 
is  3  mm  thick  silicon  steel,  whereas  the  magnetic  top  plate 
is  1.5  mm  thick  silicon  steel.  The  discharge  channel  was 
fabricated  in  two  versions,  one  constructed  of  high-purity 
alumina  ceramic  and  the  other  of  boron  nitride.  The  anode 
is  a  1.6  mm  diameter  stainless  steel  tube  with  14  propellant 
holes,  0.2  mm  in  diameter  spaced  by  1.6  mm.  A  photograph 
of  the  thruster  is  shown  in  figure  3. 


Figure  4.  Transverse  magnetic  field  strength  variation  with  axial 
position.  In  both  cases,  the  winding  current  is  1.25  A. 


Measurements  of  the  transverse  component  of  the 
magnetic  field  show  that  the  magnetic  field  near  the  anode 
is  23%  of  the  peak  value,  which  is  located  about  2  mm 
upstream  of  the  channel  exit.  The  measured  field  distribution 
for  a  winding  current  of  1 .25  A  is  shown  in  figure  4  for  two 
different  temperatures.  These  measurements  were  obtained 
ex  situ,  by  heating  the  entire  thruster  unit  in  an  oven  at  ambient 
conditions  while  measuring  the  magnetic  field  strength.  It  can 
be  seen  that  the  peak  value  of  1500  G  at  room  temperature 
drops  to  less  then  1400  G  with  a  100  °C  temperature  rise. 
This  is  significant  in  that  the  temperature  of  the  acceleration 
channel  has  been  measured  to  be  as  high  as  440  °C  by 
embedded  thermocouples  during  operation  in  the  thruster 
fabricated  out  of  boron  nitride. 

The  cathode  used  to  neutralize  the  ion  beam  and  support 
the  necessary  electric  field  is  an  Ion  Tech.  Inc.  HCN-252 
hollow  cathode.  It  is  capable  of  supplying  a  maximum  current 
of  5  A  at  xenon  flow  rates  of  0.1  to  0.5  mg  s_1 .  It  is  mounted 
in  front  of  the  thruster  such  that  the  hollow  cathode  exit  is  1 
cm  above  the  exit  of  the  channel.  The  cathode  is  the  exact 
same  unit  used  in  the  higher  power  thrusters,  and  the  flow 
rate  used  here  is  comparable  to  the  flow  rate  through  the 
thruster  itself  (2  seem).  We  expect  that  the  near  exit  plane 
xenon  density  due  to  the  cathode  flow  will  have  a  negative 
effect  on  the  discharge  performance.  However,  because  the 
neutral  gas  density  in  this  low  power  discharge  is  about  a 
factor  of  5-10  times  that  in  the  higher  power  prototype,  the 
effect  is  expected  to  be  no  greater  here  than  in  the  higher 
power  version.  No  attempt  at  designing  and  fabricating  an 
appropriately  scaled  cathode  has  been  made,  although  the 
future  development  of  low-power  (<  50  W)  Hall  thrusters 
will  rely  on  such  a  development. 

The  anode  of  the  discharge  is  powered  by  a  Sorensen 
SCR600-1.7  laboratory  power  supply  capable  of  providing 
600  V  and  1.7  A.  The  anode  also  has  a  4  Q  resistor  in 
the  power  line  to  serve  as  ballast  during  initial  start  of  the 
discharge.  The  cathode  heating  element  is  powered  by  a 
low  voltage  direct-current  (DC)  power  supply  capable  of 
providing  the  8.5  A  required  to  heat  the  cathode  for  startup 


5 


D  P  Schmidt  et  al 


Figure  5.  Photograph  of  the  linear-geometry  open  electron-drift 
Hall  discharge  while  operating. 

and  4.0  A  after  start.  The  cathode  flow  rate  of  xenon 
was  2  seem,  a  typical  value  used  during  operation  of  our 
higher  power  co-axial  discharges.  The  cathode  keeper  uses  a 
Sorensen  SCR300-6  laboratory  power  supply  to  provide  250 
V  for  initial  cathode  start  and  approximately  10  V  and  250 
mA  during  thruster  operation.  The  power  required  for  the 
magnetic  circuit  solenoids  is  provided  by  a  Tektronix  PS281 
DC  power  supply  operating  in  current  limited  mode. 

The  voltage  and  current  of  the  thruster  were  recorded  by 
acquiring  data  through  a  National  Instruments  PCI-5102  data 
acquisition  card  plugged  into  a  desktop  computer.  A  voltage 
divider  was  used  so  that  the  5  V  maximum  voltage  limit  to 
the  card  was  not  exceeded  while  testing  the  thruster  up  to  250 
V.  The  current  was  monitored  by  measuring  the  voltage  drop 
across  the  4  £2  ballast  resistor. 

6.  Results  and  analysis 

The  thruster  described  was  run  at  near-design  conditions.  To 
start  the  thruster,  a  glow  discharge  was  initiated  with  the 
magnetic  field  turned  off.  With  the  power  supply  under 
current  limit  control  and  an  upper  limit  set  on  the  voltage, 
the  magnetic  field  was  increased.  The  discharge  intensity 
greatly  increased  and  changed  colour  from  dull  pink  to  a 
bright  bluish  emission  as  the  magnetic  field  was  increased, 
which  is  also  typical  of  coaxial  designs.  A  picture  of  the 
thruster  running  at  design  conditions  is  shown  in  figure  5. 
The  voltage  gradually  increased  until  it  reached  the  voltage 
limit  setting.  The  power  supply  subsequently  switched  into 
voltage  control,  where  most  of  the  data  reported  was  taken. 

During  operation  at  magnetic  fields  above  600  G,  it  was 
noted  that  the  discharge  was  slightly  asymmetric,  being  more 
intense  near  the  side  end  of  the  channel  in  direction  of  the  E  x 
B  electron  drift.  It  was  especially  apparent  during  operation 
of  the  alumina  thruster  that  this  end  wall  became  extremely 
hot  and  glowed  intensely.  This  glow  was  contrasted  to  the 
bulk  of  the  thruster  channel,  which  did  not  show  this  intense 


heating.  On  one  occasion,  after  a  few  minutes  of  running  the 
alumina  thruster,  the  acceleration  channel  cracked  along  the 
edge  of  the  glowing  area.  This  failure  was  probably  due  to  a 
high  thermal  stress  caused  by  extreme  temperature  gradients 
in  this  region  of  the  channel  wall.  Operation  with  the 
boron  nitride  thruster  did  not  result  in  such  an  non-uniform 
temperature  field  on  the  insulating  wall.  This  difference 
between  the  boron  nitride  and  alumina  insulators  is  attributed 
to  the  difference  in  the  thermal  conductivity  values  of  the 
materials.  The  boron  nitride  channel  end  wall  in  the  direction 
of  the  electron  drift  was  instrumented  with  four  embedded 
J-type  thermocouples  distributed  along  its  length.  During 
nominal  operation  at  1500  G  and  0.7  A,  the  thermocouples 
registered  temperatures  in  excess  of  400  °C  across  the  entire 
channel,  with  the  side  wall  at  around  440  °C. 

Figure  6  shows  the  voltage  ( V)-eurrent  (/)  characteris¬ 
tics  recorded  for  both  the  boron  nitride  and  alumina  thruster 
channel  for  a  range  of  magnetic  field  strengths  and  at  a  xenon 
flow  rate  of  2  seem.  The  V-I  characteristics  for  the  boron 
nitride  channel  are  somewhat  typical  of  Hall  discharges,  with 
an  ‘ionization  branch’  at  low  currents  and  low  magnetic  field 
strengths,  and  a  relatively  steep  ‘current  saturation  branch’ 
at  high  operating  magnetic  field  strengths  and  high  discharge 
current.  However,  these  features  are  less  distinct  in  the  case 
of  the  alumina  channel.  In  fact,  the  V-I  characteristics  do  not 
show  an  obvious  current  saturation  regime  in  this  latter  case. 
At  relatively  low  magnetic  field  strengths,  the  V-I  charac¬ 
teristics  of  the  thruster  with  the  alumina  channel  wall  are 
nearly  indistinguishable  from  that  of  the  same  thruster  with  a 
boron  nitride  wall.  This  is  also  the  case  for  all  magnetic  field 
strengths  investigated,  at  currents  below  about  0.5  A.  Above 
these  current  levels,  the  V-I  characteristics  for  the  alumina 
channel  thruster  flatten  out  (note  the  apparent  ‘knee’  in  the 
figure),  whereas  the  boron  nitride  channel  thruster  voltage 
rises  sharply. 

We  speculate  that  secondary  electron  emission  from 
the  alumina  channel  is  partly  responsible  for  the  interesting 
shape  of  the  V-I  characteristics.  The  possible  influence 
of  secondary  electron  emission  in  establishing  the  electron 
transport  in  Hall  thrusters  has  been  discussed  in  the  prior 
literature  [6, 11].  Alumina  has  a  higher  secondary  electron 
emission  coefficient  than  boron  nitride  [13]  and,  in  addition, 
we  expect  that  the  secondary  electron  emission  is  sensitive 
to  wall  temperature.  It  was  apparent  that  during  thruster 
operation  with  the  alumina  walls,  the  wall  temperature 
may  have  been  sufficiently  high  to  significantly  enhance 
secondary  electron  emission.  The  higher  secondary  electron 
emission  for  the  case  of  the  alumina  insulator  wall  would 
aid  electron  transport  across  the  magnetic  field.  As  a  result, 
the  thruster  could  not  support  as  high  a  voltage  as  that 
supported  by  the  thruster  with  the  boron  nitride  wall.  This 
conjecture  would  imply  that  there  is  a  high  electron  flux 
(current)  along  the  side  wall  opposite  the  direction  of  the 
electron  Hall  current,  an  argument  that  is  consistent  with  the 
observation  that  the  side  wall  of  the  boron  nitride  thruster 
is  found  to  experience  significant  erosion.  In  a  recent  study 
by  Raitses  et  al  [6],  it  was  noted  that  there  were  significant 
differences  in  the  V-I  characteristics  of  a  thruster  operating 
with  a  machinable  glass  channel  and  a  boron  nitride  channel. 
In  that  study,  enhanced  axial  transport  in  the  thruster  with 
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Figure  2.  Schematic  of  the  linear-geometry  open  electron-drift 
Half  discharge. 


Figure  3.  Photograph  of  the  linear-geometry  open  electron-drift 
Hall  discharge. 


insulated  copper  magnet  wire.  The  magnetic  bottom  plate 
is  3  mm  thick  silicon  steel,  whereas  the  magnetic  top  plate 
is  1.5  mm  thick  silicon  steel.  The  discharge  channel  was 
fabricated  in  two  versions,  one  constructed  of  high-purity 
alumina  ceramic  and  the  other  of  boron  nitride.  The  anode 
is  a  1.6  mm  diameter  stainless  steel  tube  with  14  propellant 
holes,  0.2  mm  in  diameter  spaced  by  1 .6  mm.  A  photograph 
of  the  thruster  is  shown  in  figure  3. 


Figure  4.  Transverse  magnetic  field  strength  variation  with  axial 
position.  In  both  cases,  the  winding  current  is  1 .25  A. 


Measurements  of  the  transverse  component  of  the 
magnetic  field  show  that  the  magnetic  field  near  the  anode 
is  23%  of  the  peak  value,  which  is  located  about  2  mm 
upstream  of  the  channel  exit.  The  measured  field  distribution 
for  a  winding  current  of  1.25  A  is  shown  in  figure  4  for  two 
different  temperatures.  These  measurements  were  obtained 
ex  situ ,  by  heating  the  entire  thruster  unit  in  an  oven  at  ambient 
conditions  while  measuring  the  magnetic  field  strength.  It  can 
be  seen  that  the  peak  value  of  1500  G  at  room  temperature 
drops  to  less  then  1400  G  with  a  100  °C  temperature  rise. 
This  is  significant  in  that  the  temperature  of  the  acceleration 
channel  has  been  measured  to  be  as  high  as  440  °C  by 
embedded  thermocouples  during  operation  in  the  thruster 
fabricated  out  of  boron  nitride. 

The  cathode  used  to  neutralize  the  ion  beam  and  support 
the  necessary  electric  field  is  an  Ion  Tech.  Inc.  HCN-252 
hollow  cathode.  It  is  capable  of  supplying  a  maximum  current 
of  5  A  at  xenon  flow  rates  of  0.1  to  0.5  mg  s'"1 .  It  is  mounted 
in  front  of  the  thruster  such  that  the  hollow  cathode  exit  is  1 
cm  above  the  exit  of  the  channel.  The  cathode  is  the  exact 
same  unit  used  in  the  higher  power  thrusters,  and  the  flow 
rate  used  here  is  comparable  to  the  flow  rate  through  the 
thruster  itself  (2  seem).  We  expect  that  the  near  exit  plane 
xenon  density  due  to  the  cathode  flow  will  have  a  negative 
effect  on  the  discharge  performance.  However,  because  the 
neutral  gas  density  in  this  low  power  discharge  is  about  a 
factor  of  5-10  times  that  in  the  higher  power  prototype,  the 
effect  is  expected  to  be  no  greater  here  than  in  the  higher 
power  version.  No  attempt  at  designing  and  fabricating  an 
appropriately  scaled  cathode  has  been  made,  although  the 
future  development  of  low-power  (<  50  W)  Hall  thrusters 
will  rely  on  such  a  development. 

The  anode  of  the  discharge  is  powered  by  a  Sorensen 
SCR600-1.7  laboratory  power  supply  capable  of  providing 
600  V  and  1.7  A.  The  anode  also  has  a  4  ft  resistor  in 
the  power  line  to  serve  as  ballast  during  initial  start  of  the 
discharge.  The  cathode  heating  element  is  powered  by  a 
low  voltage  direct-current  (DC)  power  supply  capable  of 
providing  the  8.5  A  required  to  heat  the  cathode  for  startup 
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Figure  5.  Photograph  of  the  linear-geometry  open  electron-drift 
Hall  discharge  while  operating. 


and  4.0  A  after  start.  The  cathode  flow  rate  of  xenon 
was  2  seem,  a  typical  value  used  during  operation  of  our 
higher  power  co-axial  discharges.  The  cathode  keeper  uses  a 
Sorensen  SCR300-6  laboratory  power  supply  to  provide  250 
V  for  initial  cathode  start  and  approximately  10  V  and  250 
mA  during  thruster  operation.  The  power  required  for  the 
magnetic  circuit  solenoids  is  provided  by  a  Tektronix  PS281 
DC  power  supply  operating  in  current  limited  mode. 

The  voltage  and  current  of  the  thruster  were  recorded  by 
acquiring  data  through  a  National  Instruments  PCI-5 102  data 
acquisition  card  plugged  into  a  desktop  computer.  A  voltage 
divider  was  used  so  that  the  5  V  maximum  voltage  limit  to 
the  card  was  not  exceeded  while  testing  the  thruster  up  to  250 
V.  The  current  was  monitored  by  measuring  the  voltage  drop 
across  the  4  Q  ballast  resistor. 

6.  Results  and  analysis 

The  thruster  described  was  run  at  near-design  conditions.  To 
start  the  thruster,  a  glow  discharge  was  initiated  with  the 
magnetic  field  turned  off.  With  the  power  supply  under 
current  limit  control  and  an  upper  limit  set  on  the  voltage, 
the  magnetic  field  was  increased.  The  discharge  intensity 
greatly  increased  and  changed  colour  from  dull  pink  to  a 
bright  bluish  emission  as  the  magnetic  field  was  increased, 
which  is  also  typical  of  coaxial  designs.  A  picture  of  the 
thruster  running  at  design  conditions  is  shown  in  figure  5. 
The  voltage  gradually  increased  until  it  reached  the  voltage 
limit  setting.  The  power  supply  subsequently  switched  into 
voltage  control,  where  most  of  the  data  reported  was  taken. 

During  operation  at  magnetic  fields  above  600  G,  it  was 
noted  that  the  discharge  was  slightly  asymmetric,  being  more 
intense  near  the  side  end  of  the  channel  in  direction  of  the  E  x 
B  electron  drift.  It  was  especially  apparent  during  operation 
of  the  alumina  thruster  that  this  end  wall  became  extremely 
hot  and  glowed  intensely.  This  glow  was  contrasted  to  the 
bulk  of  the  thruster  channel,  which  did  not  show  this  intense 


heating.  On  one  occasion,  after  a  few  minutes  of  running  the 
alumina  thruster,  the  acceleration  channel  cracked  along  the 
edge  of  the  glowing  area.  This  failure  was  probably  due  to  a 
high  thermal  stress  caused  by  extreme  temperature  gradients 
in  this  region  of  the  channel  wall.  Operation  with  the 
boron  nitride  thruster  did  not  result  in  such  an  non-uniform 
temperature  field  on  the  insulating  wall.  This  difference 
between  the  boron  nitride  and  alumina  insulators  is  attributed 
to  the  difference  in  the  thermal  conductivity  values  of  the 
materials.  The  boron  nitride  channel  end  wall  in  the  direction 
of  the  electron  drift  was  instrumented  with  four  embedded 
J-type  thermocouples  distributed  along  its  length.  During 
nominal  operation  at  1500  G  and  0.7  A,  the  thermocouples 
registered  temperatures  in  excess  of  400  °C  across  the  entire 
channel,  with  the  side  wall  at  around  440  °C. 

Figure  6  shows  the  voltage  (V)-current  (/)  characteris¬ 
tics  recorded  for  both  the  boron  nitride  and  alumina  thruster 
channel  for  a  range  of  magnetic  field  strengths  and  at  a  xenon 
flow  rate  of  2  seem.  The  V-I  characteristics  for  the  boron 
nitride  channel  are  somewhat  typical  of  Hall  discharges,  with 
an  ‘ionization  branch’  at  low  currents  and  low  magnetic  field 
strengths,  and  a  relatively  steep  ‘current  saturation  branch’ 
at  high  operating  magnetic  field  strengths  and  high  discharge 
current.  However,  these  features  are  less  distinct  in  the  case 
of  the  alumina  channel.  In  fact,  the  V-7  characteristics  do  not 
show  an  obvious  current  saturation  regime  in  this  latter  case. 
At  relatively  low  magnetic  field  strengths,  the  V-I  charac¬ 
teristics  of  the  thruster  with  the  alumina  channel  wall  are 
nearly  indistinguishable  from  that  of  the  same  thruster  with  a 
boron  nitride  wall.  This  is  also  the  case  for  all  magnetic  field 
strengths  investigated,  at  currents  below  about  0.5  A.  Above 
these  current  levels,  the  V-I  characteristics  for  the  alumina 
channel  thruster  flatten  out  (note  the  apparent  ‘knee’  in  the 
figure),  whereas  the  boron  nitride  channel  thruster  voltage 
rises  sharply. 

We  speculate  that  secondary  electron  emission  from 
the  alumina  channel  is  partly  responsible  for  the  interesting 
shape  of  the  V-I  characteristics.  The  possible  influence 
of  secondary  electron  emission  in  establishing  the  electron 
transport  in  Hall  thrusters  has  been  discussed  in  the  prior 
literature  [6, 11].  Alumina  has  a  higher  secondary  electron 
emission  coefficient  than  boron  nitride  [13]  and,  in  addition, 
we  expect  that  the  secondary  electron  emission  is  sensitive 
to  wall  temperature.  It  was  apparent  that  during  thruster 
operation  with  the  alumina  walls,  the  wall  temperature 
may  have  been  sufficiently  high  to  significantly  enhance 
secondary  electron  emission.  The  higher  secondary  electron 
emission  for  the  case  of  the  alumina  insulator  wall  would 
aid  electron  transport  across  the  magnetic  field.  As  a  result, 
the  thruster  could  not  support  as  high  a  voltage  as  that 
supported  by  the  thruster  with  the  boron  nitride  wall.  This 
conjecture  would  imply  that  there  is  a  high  electron  flux 
(current)  along  the  side  wall  opposite  the  direction  of  the 
electron  Hall  current,  an  argument  that  is  consistent  with  the 
observation  that  the  side  wall  of  the  boron  nitride  thruster 
is  found  to  experience  significant  erosion.  In  a  recent  study 
by  Raitses  et  al  [6],  it  was  noted  that  there  were  significant 
differences  in  the  V-I  characteristics  of  a  thruster  operating 
with  a  machinable  glass  channel  and  a  boron  nitride  channel. 
In  that  study,  enhanced  axial  transport  in  the  thruster  with 
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Figure  6.  Discharge  voltage-current  characteristics  of  the 
linear-geometry  Hall  discharge  (top)  boron  nitride  channel, 
(bottom)  alumina  channel. 

the  machinable  glass  channel  was  attributed  to  wall  effects, 
decreasing  the  thruster  efficiency  at  high  operating  voltages. 
The  qualitative  findings  reported  on  in  this  study  agree 
with  these  past  observations.  It  is  also  noteworthy  that  the 
conjecture  that  there  is  an  enhanced  electron  current  due  to 
wall  collisions  is  supported  by  the  relatively  poor  efficiency 
of  this  discharge,  as  an  upper  limit  of  the  ratio  of  the  ion 
current  to  electron  current  is  no  more  than  30%,  based  on 
full  utilization  (ionization)  of  the  propellant. 

In  all  cases,  the  thrusters  exhibited  a  region  of  unstable 
operation.  At  all  magnetic  field  strengths  except  the  lowest 
value  shown,  the  discharge  had  a  low  current  limit  in  its 
operating  envelope.  Attempts  to  operate  or  start  the  discharge 
in  this  region  would  fail.  We  speculate  that  this  region  of 
instability  is  closely  tied  to  the  requirement  for  enhanced 
electron  transport  and  that,  at  low  currents  and  high  magnetic 
fields,  the  anomalous  transport  process  cannot  provide  the 
necessary  current  to  maintain  the  discharge. 

The  fact  that  the  thruster  ran  without  a  mechanism 
for  closing  the  Hall  current  confirms  the  importance  of  an 
electron  transport  process  that  is  due  to  plasma  fluctuations 
and/or  wall  effects  associated  with  secondary  electron 
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Figure  7.  Time-variation  in  discharge  current  for  the 
alumina-walled  plasma  source  at  a  range  of  peak  magnetic  field 
conditions. 


emission.  As  in  the  higher  power  co-axial  discharges, 
the  linear  discharge  studied  here  also  exhibited  plasma 
fluctuations,  which  were  detected  as  fluctuations  in  the 
external  circuit  discharge  current. 

Figure  7  shows  the  oscillations  in  the  discharge  current 
of  the  linear  thruster  operating  with  the  alumina  channel  wall 
for  a  range  of  peak  magnetic  field  strengths  and  at  a  discharge 
current  of  0.7  A.  The  fluctuations  in  the  discharge  current  for 
the  boron  nitride  wall  were  qualitatively  similar.  It  is  apparent 
that  at  low  magnetic  fields,  there  is  a  relatively  low  frequency 
oscillation,  on  which  higher  frequencies  are  superimposed. 
The  low  frequency  oscillation  increases  in  amplitude  and  in 
frequency  as  the  magnetic  field  is  increased.  These  intense 
low  frequencies  observed  in  this  linear  device  are  similar  to 
those  seen  in  co-axial  devices,  and  are  believed  to  be  the 
so-called  ‘breathing’  mode  of  oscillation  associated  with  the 
neutral  xenon  transit  through  the  ionization  zone  [14].  This 
instability  is  associated  with  the  disturbance  in  the  balance 
established  between  the  depletion  of  neutrals  in  the  channel 
as  a  result  of  ionization,  and  their  replenishment.  Since  the 
length  of  the  ionization  zone  in  this  low -power  Hall  discharge 
is  scaled  to  be  some  5-10  times  shorter  than  that  of  our 
reference  thruster,  the  frequencies  of  these  disturbances  are 
expected  to  be  at  least  a  factor  of  five  higher  than  those  seen 
in  our  higher  power  co-axial  devices. 

Figure  8  compares  the  Fourier  analysis  of  the  temporal 
fluctuations  in  discharge  current  of  the  reference  400  W  (200 
V,  20  seem,  160  G,  2  A)  co-axial  thruster  and  the  low- 
power  linear  thruster  operating  with  the  alumina  channel 
wall  (data  shown  in  figure  7).  It  is  apparent  that  the  low- 
power  thruster  has  a  strong  low-frequency  mode,  similar  to 
that  seen  in  the  co-axial  high-  power  devices,  at  frequencies 
that  are  approximately  a  factor  of  four  times  that  of  the 
higher  power  thruster.  As  the  magnetic  field  is  increased 
the  frequency  of  the  fluctuations  in  the  anode  current 
also  increase,  until  a  magnetic  field  strength  of  900  G, 
beyond  which  it  remains  constant.  This  result  is  seemingly 
inconsistent  with  the  theoretical  predictions  of  Boeuf  and 
Garrigues  [14].  However,  a  direct  comparison  to  the  results 
in  [14]  is  difficult  to  make,  since  in  our  studies,  the  current  is 
held  constant  while  the  magnetic  field  is  increased  (resulting 
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Figure  8.  Low-frequency  spectral  analysis  of  the  temporal 
fluctuations  in  the  discharge  current.  Included  in  the  figure  is  the 
spectral  distribution  of  the  fluctuations  seen  in  the  higher  power 
co-axial  discharge  for  comparison. 

in  increased  discharge  voltages).  In  the  calculations  of  Boeuf 
and  Garrigues,  the  voltage  was  varied  at  constant  magnetic 
field  (giving  rise  to  varying  current)  and/or  the  magnetic  field 
was  varied  at  constant  voltage.  As  discussed  in  [14],  the 
frequency  of  this  mode  is  seen  to  increase  dramatically  with 
voltage  (at  constant  magnetic  field).  The  response  that  we 
see  is  therefore  likely  to  be  a  result  of  the  response  in  the 
frequency  to  changes  in  both  the  voltage  and  the  magnetic 
field. 

It  is  interesting  to  note  that  at  the  highest  magnetic  field 
studied,  there  is  a  superimposed  high-frequency  oscillation 
at  about  80  kHz.  Through  particle  simulations,  Boeuf 
and  Garrigues  [14]  also  discovered  the  presence  of  strong 
disturbances  in  the  plasma  density  and  electric  field  upstream 
near  the  anode  when  the  magnetic  field  continues  to  persist 
near  the  anode.  The  weak  intensities  seen  here  at  the 
highest  magnetic  field  studied  may  be  a  consequence  of  these 
near-  anode  instabilities,  although  a  precise  characterization 
of  these  instabilities  in  a  linear  geometry  must  still  be 
performed.  In  a  related  study,  we  have  identified  the  presence 
of  near-anode  instabilities  in  this  frequency  range  (40-100 
kHz)  in  co-axial  Hall  discharges  [15]. 

7.  Summary 

An  analysis  and  arguments  for  the  scaling  of  a  Hall  plasma 
thruster  to  low  powers  was  presented.  A  linear-geometry 
version  of  a  thruster  operating  in  the  50-100  W  power  range 
has  been  fabricated  and  operated  at  near-design  conditions. 
Preliminary  results  obtained  so  far  indicate  that  at  the  scaled 
power  levels,  these  low  power  plasma  discharges  operate  at 
much  higher  channel  wall  temperatures.  This  increased  heat 
flux  may  be  a  major  impediment  to  the  extended  operation 
of  very  low  power  devices. 

The  linear  Hall  plasma  thruster  reported  on  here  is  found 
to  have  the  characteristic  discharge  instabilities  seen  in  higher 
power  co-axial  versions.  Although  we  have  not  performed  an 


extensive  analysis  of  the  thruster  performance,  it  does  appear 
that  the  linear  device  behaves  in  many  ways  similarly  to  those 
of  a  co-axial  design  with  a  closed  electron  drift.  Since  the 
linear  discharge  operates  without  a  closed  Hall  current,  it 
suggests  that  there  must  be  an  anomalous  mechanism  for 
cross-field  electron  transport. 

This  device  might  prove  most  useful  for  investigating 
various  materials  for  use  as  acceleration  channels  in  Hall 
thrusters.  The  linear-geometry  design  allows  easy  fabrication 
from  a  variety  of  materials.  A  study  of  the  operating 
characteristics  of  thrusters  constructed  with  insulating  walls 
fabricated  from  a  wide  variety  of  materials  would  be  useful 
for  understanding  the  effect  of  secondary  electron  emission 
on  electron  transport. 
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Figure  6.  Discharge  voltage-current  characteristics  of  the 
linear-geometry  Hall  discharge  (top)  boron  nitride  channel, 
(bottom)  alumina  channel. 

the  machinable  glass  channel  was  attributed  to  wall  effects, 
decreasing  the  thruster  efficiency  at  high  operating  voltages. 
The  qualitative  findings  reported  on  in  this  study  agree 
with  these  past  observations.  It  is  also  noteworthy  that  the 
conjecture  that  there  is  an  enhanced  electron  current  due  to 
wall  collisions  is  supported  by  the  relatively  poor  efficiency 
of  this  discharge,  as  an  upper  limit  of  the  ratio  of  the  ion 
current  to  electron  current  is  no  more  than  30%,  based  on 
full  utilization  (ionization)  of  the  propellant. 

In  all  cases,  the  thrusters  exhibited  a  region  of  unstable 
operation.  At  all  magnetic  field  strengths  except  the  lowest 
value  shown,  the  discharge  had  a  low  current  limit  in  its 
operating  envelope.  Attempts  to  operate  or  start  the  discharge 
in  this  region  would  fail.  We  speculate  that  this  region  of 
instability  is  closely  tied  to  the  requirement  for  enhanced 
electron  transport  and  that,  at  low  currents  and  high  magnetic 
fields,  the  anomalous  transport  process  cannot  provide  the 
necessary  current  to  maintain  the  discharge. 

The  fact  that  the  thruster  ran  without  a  mechanism 
for  closing  the  Hall  current  confirms  the  importance  of  an 
electron  transport  process  that  is  due  to  plasma  fluctuations 
and/or  wall  effects  associated  with  secondary  electron 
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Figure  7.  Time-variation  in  discharge  current  for  the 
alumina-walled  plasma  source  at  a  range  of  peak  magnetic  field 
conditions. 


emission.  As  in  the  higher  power  co-axial  discharges, 
the  linear  discharge  studied  here  also  exhibited  plasma 
fluctuations,  which  were  detected  as  fluctuations  in  the 
external  circuit  discharge  current. 

Figure  7  shows  the  oscillations  in  the  discharge  current 
of  the  linear  thruster  operating  with  the  alumina  channel  wall 
for  a  range  of  peak  magnetic  field  strengths  and  at  a  discharge 
current  of  0.7  A.  The  fluctuations  in  the  discharge  current  for 
the  boron  nitride  wall  were  qualitatively  similar.  It  is  apparent 
that  at  low  magnetic  fields,  there  is  a  relatively  low  frequency 
oscillation,  on  which  higher  frequencies  are  superimposed. 
The  low  frequency  oscillation  increases  in  amplitude  and  in 
frequency  as  the  magnetic  field  is  increased.  These  intense 
low  frequencies  observed  in  this  linear  device  are  similar  to 
those  seen  in  co-axial  devices,  and  are  believed  to  be  the 
so-called  ‘breathing’  mode  of  oscillation  associated  with  the 
neutral  xenon  transit  through  the  ionization  zone  [14].  This 
instability  is  associated  with  the  disturbance  in  the  balance 
established  between  the  depletion  of  neutrals  in  the  channel 
as  a  result  of  ionization,  and  their  replenishment.  Since  the 
length  of  the  ionization  zone  in  this  low-power  Hall  discharge 
is  scaled  to  be  some  5-10  times  shorter  than  that  of  our 
reference  thruster,  the  frequencies  of  these  disturbances  are 
expected  to  be  at  least  a  factor  of  five  higher  than  those  seen 
in  our  higher  power  co-axial  devices. 

Figure  8  compares  the  Fourier  analysis  of  the  temporal 
fluctuations  in  discharge  current  of  the  reference  400  W  (200 
V,  20  seem,  160  G,  2  A)  co-axial  thruster  and  the  low- 
power  linear  thruster  operating  with  the  alumina  channel 
wall  (data  shown  in  figure  7).  It  is  apparent  that  the  low- 
power  thruster  has  a  strong  low-frequency  mode,  similar  to 
that  seen  in  the  co-axial  high-  power  devices,  at  frequencies 
that  are  approximately  a  factor  of  four  times  that  of  the 
higher  power  thruster.  As  the  magnetic  field  is  increased 
the  frequency  of  the  fluctuations  in  the  anode  current 
also  increase,  until  a  magnetic  field  strength  of  900  G, 
beyond  which  it  remains  constant.  This  result  is  seemingly 
inconsistent  with  the  theoretical  predictions  of  Boeuf  and 
Garrigues  [14].  However,  a  direct  comparison  to  the  results 
in  [14]  is  difficult  to  make,  since  in  our  studies,  the  current  is 
held  constant  while  the  magnetic  field  is  increased  (resulting 
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Figure  8.  Low-frequency  spectral  analysis  of  the  temporal 
fluctuations  in  the  discharge  current.  Included  in  the  figure  is  the 
spectral  distribution  of  the  fluctuations  seen  in  the  higher  power 
co-axial  discharge  for  comparison. 

in  increased  discharge  voltages).  In  the  calculations  of  Boeuf 
and  Garrigues,  the  voltage  was  varied  at  constant  magnetic 
field  (giving  rise  to  varying  current)  and/or  the  magnetic  field 
was  varied  at  constant  voltage.  As  discussed  in  [14],  the 
frequency  of  this  mode  is  seen  to  increase  dramatically  with 
voltage  (at  constant  magnetic  field).  The  response  that  we 
see  is  therefore  likely  to  be  a  result  of  the  response  in  the 
frequency  to  changes  in  both  the  voltage  and  the  magnetic 
field. 

It  is  interesting  to  note  that  at  the  highest  magnetic  field 
studied,  there  is  a  superimposed  high-frequency  oscillation 
at  about  80  kHz.  Through  particle  simulations,  Boeuf 
and  Garrigues  [14]  also  discovered  the  presence  of  strong 
disturbances  in  the  plasma  density  and  electric  field  upstream 
near  the  anode  when  the  magnetic  field  continues  to  persist 
near  the  anode.  The  weak  intensities  seen  here  at  the 
highest  magnetic  field  studied  may  be  a  consequence  of  these 
near-  anode  instabilities,  although  a  precise  characterization 
of  these  instabilities  in  a  linear  geometry  must  still  be 
performed.  In  a  related  study,  we  have  identified  the  presence 
of  near- anode  instabilities  in  this  frequency  range  (40-100 
kHz)  in  co-axial  Hall  discharges  [15]. 

7.  Summary 

An  analysis  and  arguments  for  the  scaling  of  a  Hall  plasma 
thruster  to  low  powers  was  presented.  A  linear-geometry 
version  of  a  thruster  operating  in  the  50-100  W  power  range 
has  been  fabricated  and  operated  at  near-design  conditions. 
Preliminary  results  obtained  so  far  indicate  that  at  the  scaled 
power  levels,  these  low  power  plasma  discharges  operate  at 
much  higher  channel  wall  temperatures.  This  increased  heat 
flux  may  be  a  major  impediment  to  the  extended  operation 
of  very  low  power  devices. 

The  linear  Hall  plasma  thruster  reported  on  here  is  found 
to  have  the  characteristic  discharge  instabilities  seen  in  higher 
power  co-axial  versions.  Although  we  have  not  performed  an 


extensive  analysis  of  the  thruster  performance,  it  does  appear 
that  the  linear  device  behaves  in  many  ways  similarly  to  those 
of  a  co-axial  design  with  a  closed  electron  drift.  Since  the 
linear  discharge  operates  without  a  closed  Hall  current,  it 
suggests  that  there  must  be  an  anomalous  mechanism  for 
cross-field  electron  transport. 

This  device  might  prove  most  useful  for  investigating 
various  materials  for  use  as  acceleration  channels  in  Hall 
thrusters.  The  linear-geometry  design  allows  easy  fabrication 
from  a  variety  of  materials.  A  study  of  the  operating 
characteristics  of  thrusters  constructed  with  insulating  walls 
fabricated  from  a  wide  variety  of  materials  would  be  useful 
for  understanding  the  effect  of  secondary  electron  emission 
on  electron  transport. 
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